Ga-doping induced carrier tuning and multiphase engineering in n

Publication Date (Web): June 12, 2018. Copyright © 2018 American Chemical Society. Cite this:ACS Appl. Mater. Interfaces XXXX, XXX, XXX-XXX ...
8 downloads 0 Views 1MB Size
Subscriber access provided by University of Winnipeg Library

Functional Inorganic Materials and Devices

Ga-doping induced carrier tuning and multiphase engineering in n-type PbTe with enhanced thermoelectric performance Zhengshang Wang, Guoyu Wang, Ruifeng Wang, Xiaoyuan Zhou, Zhiyu Chen, Cong Yin, Mingjing Tang, Qing Hu, Jun Tang, and Ran Ang ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.8b05117 • Publication Date (Web): 12 Jun 2018 Downloaded from http://pubs.acs.org on June 12, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Ga-doping induced carrier tuning and multiphase engineering in n-type PbTe with enhanced thermoelectric performance Zhengshang Wanga, Guoyu Wangb,c, Ruifeng Wangb,c, Xiaoyuan Zhoud, Zhiyu Chena, Cong Yina, Mingjing Tanga, Qing Hua, Jun Tanga,*, Ran Anga,e,*

a

Key Laboratory of Radiation Physics and Technology, Ministry of Education, Institute of Nuclear Science and Technology, Sichuan University, Chengdu 610064, China

b

Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences, Chongqing 400714, China c

d

e

University of Chinese Academy of Sciences, Beijing 100190, China

College of Physics, Chongqing University, Chongqing 401331, China

Institute of New Energy and Low-Carbon Technology, Sichuan University, Chengdu 610065, China

ACS Paragon Plus Environment

1

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 26

KEYWORDS: Thermoelectric materials; n-type PbTe; Ga doping; Carrier tuning; Multiphase engineering

ABSTRACT: P-type lead telluride (PbTe) emerged as a promising thermoelectric material for intermediate temperature waste heat energy harvesting. However, n-type PbTe still confronted with a considerable challenge owing to its relatively low figure of merit ZT and conversion efficiency η, limiting widespread thermoelectric applications. Here, we report that Ga-doping in n-type PbTe can optimize carrier concentration and thus improve power factor. Moreover, further experimental and theoretical evidences reveal that Ga-doping induced multiphase structures with nano- to micrometer size can simultaneously modulate phonon transport, leading to dramatic reduction of lattice thermal conductivity. As a consequence, a tremendous enhancement of ZT value at 823 K reaches ~1.3 for n-type Pb0.97Ga0.03Te. In particular, in a wide temperature range from 323 to 823 K, the average ZTave value of ~0.9 and calculated conversion efficiency η of ~13% are achieved by Ga doping. The present findings demonstrate the great potential in Ga-doped PbTe thermoelectric materials through a synergetic carrier tuning and multiphase engineering strategy.

1. Introduction Thermoelectric materials have been recognized as an ideal candidate in power generation and solid-state cooling, which can directly convert waste heat into electricity and vice versa.1-4 The efficiency of thermoelectric material is governed by the dimensionless figure of merit ZT =

ACS Paragon Plus Environment

2

Page 3 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

S2σT/κtot, where S is the Seebeck coefficient, σ is the electrical conductivity, and κtot is the total thermal conductivity which comprises the lattice and electronic components κlat and κele.5,6 It is apparent that a high ZT requires a high S, high σ, and low κtot at a certain temperature. Unfortunately, these parameters S, σ, and κtot are closely interrelated with each other.7-9 Thus, how to modulate harmoniously the electrical and thermal properties of high performance thermoelectric materials is a crucial issue to be resolved.10-12 As one of the state-of-the-art thermoelectric materials for power generation at the intermediate temperature (500-900 K),13 lead telluride (PbTe) has an inherent ZT value of 0.8.10 For p-type PbTe, the electrical properties (S and σ) have been governed by carrier concentration tuning,14 band convergence,15 and band energy alignment.16 On the other hand, the thermal properties (κlat) of p-type PbTe have also been determined by introducing point defects,17 nanostructuring,18,19 high density of dislocations,20 interfaces,21 and all-scale hierarchical architecturing.22 Consequently, the thermoelectric performance exhibits ultrahigh ZT values for p-type PbTe materials. However, the ZT values for n-type PbTe materials are relatively lower than that of ptype PbTe materials.23-29 From the viewpoint of commercial applications, the implementation of thermoelectric energy conversion requires an identical thermoelectric performance between nand p-type PbTe materials.30 Therefore, it is urgent for n-type PbTe materials to enhance ZT merit to break through the limitation of large-scale application. Current strategies to improve thermoelectric performance of n-type PbTe, which is focused on enhancing power factor PF (PF = S2σ) and suppressing κlat, can be mainly summarized in the following aspects: (1) self-tuning carrier concentration with temperature-dependent solubility of impurity in the matrix,31 (2) increasing effective mass m* by producing impurity levels which locally distort density of states (DOS) of conduction band and form so-called resonant states,32

ACS Paragon Plus Environment

3

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 26

(3) scattering low-energy carriers with multiphase potential barriers resulting in a record-high ZT for PbTe-4%InSb,29 (4) decreasing κlat by synergistic alloying and nanostructuring.33-37 However, few novel methods for n-type PbTe have been developed by multiphase engineering to suppress κlat. Particularly, it is not yet known that individual contribution of various phases to reduce κlat and to what extent the multiphase can play a role in manipulating transport properties of carriers and phonons. This motivates us to seek for a possible avenue by multiphase engineering to enhancing ZT values, as well as relatively high average ZT (ZTave) and conversion efficiency (η) over a wide temperature range. In this work, we report that Ga-doping-driven n-type PbTe not only promotes a noticeable enhancement of power factor by carrier tuning, but also induces multiphase structures with nanoto micrometer size which simultaneously manipulates phonon transport. Ultimately, the maximum of ZT value for n-type Pb0.97Ga0.03Te reaches ~1.3 at 823 K, as well as the average ZTave value of ~0.9 and calculated conversion efficiency η of ~13% in a wide temperature region from 323 to 823 K. 2. Experimental Polycrystal Pb1-xGaxTe (x = 0.01, 0.02, 0.03, 0.04) samples were fabricated by convention solid-state reaction technique. Stoichiometric quantities of high purity elements [lead (99.99%), tellurium (99.999%), and gallium (99.999%)] were loaded into quartz tubes. The tubes then were evacuated, flame-sealed, slowly heated to 773 K, soaked at this temperature for 6 h, followed by melting at 1273 K for 6 h, finally cooled to room temperature rapidly. The obtained ingots were crushed into powders and densified by spark plasma sintering (SPS) method (LABOX-325, Japan) at 823 K for 10 min with a 60 MPa uniaxial pressure. All the disk shaped samples

ACS Paragon Plus Environment

4

Page 5 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

obtained were 10 mm in diameter with density no less than 97% of theoretical density (Table S1). Powder X-ray diffraction (PXRD) measurements with Cu-Kα radiation of powder samples after SPS process were performed on an X-ray diffractometer (DX-2700). Structural refinements were carried out by Rietveld analysis. The surface morphology was characterized by a FEI inspect F50 scanning electron microscope (SEM). The temperature dependence of electrical conductivity σ and Seebeck coefficient S were carried out by using apparatus (LSR-3). The temperature dependence of thermal diffusivity D and thermal conductivity κ were performed by the laser flash method (LFA 457, Netzsch). The heat capacity Cp was determined from the measurements

of

Blachnik

by Cp(kB/atom) = 3.07

+

0.00047(T/K-300) for

lead

chalcogenides.14,15 To verify the accuracy of Cp from the calculation, the experimental Cp was also measured using a TA Q20 differential scanning calorimeter (DSC). The uncertainties in ZT are nearly ±10%. The Hall coefficient RH and carrier concentration nH were determined by Hall measurements using a standard five-probe method and conducted by the Quantum Design Physical Property Measurement System (PPMS). 3. Results and discussion Figure 1a shows the PXRD patterns for Pb1-xGaxTe (x = 0.01, 0.02, 0.03, 0.04) samples, in which the main diffraction peaks in all samples are well indexed to face-centre-cubic PbTe structure. Interestingly, for the samples with x = 0.03 and 0.04, the diffraction peaks of impurity marked by the asterisks are observed at 26°, 29° and 45°, verifying the appearance of Ga6PbTe10 phase (inset of Figure 1a). In addition, we also carried out structural Rietveld refinements to reveal the variation of impurity phase (Figure S1), and the detailed phase fraction of Ga6PbTe10

ACS Paragon Plus Environment

5

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 26

is plotted in Figure S2. It is noted that no diffraction peak of Ga6PbTe10 appears in XRD patterns when x < 0.03, since the phase fraction of Ga6PbTe10 is much lower than the limited resolution of XRD technology (about 1%), which can also be verified by the following SEM results. While for x ≥ 0.03, the phase fraction of Ga6PbTe10 increases from 1.1% to 2.8%, triggering the emergence of impurity peaks as shown in the inset of Figure 1a. To probe thermoelectric performance of pristine and Ga-doped PbTe, we first measured electrical properties as a function of temperature. Evidently, Ga doping improves the electrical conductivity σ due to the increase of carrier density nH (Figure 1b and Figure 1c). Therefore, the optimal carrier concentration can be achieved by manipulating Ga content. With rising temperature, σ for each Ga-doped sample decreases (Figure 1c), exhibiting a typical characteristic of degenerate semiconductor. One might have noticed that temperature-dependent σ increases with increasing of Ga content from x = 0 to 0.03, while it drops for x = 0.04, taking into account the lower mobility µH (Figure 1b). In particular, at low temperatures, Seebeck coefficient S for Ga-doped PbTe presents a rising feature in Figure 2d, further supporting the degenerate behavior. The slow-down characteristic of S at high temperatures originates from the thermal excitation of minority carriers. It is worth noting that with increasing of Ga content, the absolute value of S declines because of the increased nH (Figure 1b). In addition, the S for pristine PbTe shows an apparent transition from p- to n-type occurring around 550 K (inset of Figure 3d), which is accordant with previous results,38 indicating the comparable carrier concentrations between electrons and holes. Especially for the sample with x = 0.03, carrier tuning by Ga doping induces a remarkable improvement of power factor PF (Figure 1e). While for the sample with x = 0.04, it has a lower carrier mobility µH (Figure 1b), resulting in the depression of PF in the entire temperature range. Theoretical Pisarenko relation based on single

ACS Paragon Plus Environment

6

Page 7 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Kane band model (SKB) and acoustic phonon scattering can explain experimental data well.39,40 The S at room temperature as a function of nH (Figure 1f) demonstrates that the experimental data matches well with the theoretical predictions of n-type PbTe (effective mass m* = 0.25 me). Such results reveal that Ga doping do not significantly affect the conduction band structure. Note that the magnitudes of nH for all samples reach the order of ~1019 cm-3. However, the low values of nH do not follow with the theoretical assumption, i.e. each Ga atom occupying the Pb site would offer one electron. Hence, a large number of Ga atoms can be favorable to the formation of precipitates when Ga point defects exceed the limit of solubility, which is also supported by the diffraction peaks of impurity marked by the asterisks (Figure 1a). Thus, it is necessary to examine the evolution of microstructures by Ga doping. The characterization of SEM and energy dispersive spectroscopy (EDS) for all of samples are shown in Figure 2. As for x = 0.01, the matrix can be clearly distinguished due to the distinct grains and boundaries (Figure 2a). The black particles marked by white circles labeled as phase 1 (P1), are different from that of pristine PbTe. Additionally, another phase distributes in the grain boundaries of PbTe with nano-scale size denoted as phase 2 (P2). As for x = 0.02, the content of P2 precipitates visibly increase (Figure 2b). In contrast, the black particles P1 start to assemble (marked by the white arrowheads in Figure 2b). Nevertheless, the atomic ratio is very close between P1 and matrix (i.e. Pb : Ga : Te = 41.15 : 52.57 : 6.28 for P1 and Pb : Ga : Te = 43.25 : 52.40 : 4.35 for matrix). Therefore, it is reasonable that both P1 and matrix should have the same contribution to scattering of phonons and carriers, as further verified in Figure S3. As for x = 0.03 and 0.04, a new kind of phase with micrometer scale emerges, denoted as P3, which is surrounded by grains and P1 (Figure 2c and Figure 2d). More details about P1, P2, and P3 are described in Figure S3. In comparison with matrix, Ga element is excessive in P2, while

ACS Paragon Plus Environment

7

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 26

both Ga and Te elements are excessive in P3. These multiple phases with scales from nanometer to micrometer provide a new route to modulate thermoelectric properties, especially for phonon transport. In order to understand the intrinsic relationship between phonon scattering and multiphase, we measured the thermal transport properties. The detailed thermal diffusivity and heat capacity are listed in Figure S4. Figure 3a shows the temperature dependence of total thermal conductivity κtot. Due to intense Umklapp scattering, the calculated lattice thermal conductivity κlat rapidly decreases in a wide temperature range (Figure 3b). To gain more insight into the contribution of κlat, we firstly simulated the values of κlat at room temperature (Figure 3c) in the framework of modified Klemen’s model,41,42 which only includes Umklapp process and point defects as the scattering sources in thermal transport. At relatively low Ga content (x ≤ 0.01), the experimental values of κlat follow fairly well with the theoretical prediction, revealing that point defects dominate phonon scattering for the reduction of κlat. While for x ≥ 0.02, the experimental values of κlat gradually deviate from the calculated lines (Figure 3c), demonstrating the strong phonon scattering arising from various multiple phases. Herein, we utilized a well-accepted model based on Debye approximation43 to simulate the experimental values of κlat. The theoretical κlat can be expressed as:

κ lat

k k T  = B2  B  2π ν  h 

3



θ D /T

0

τc ( x)

x 4e x

(e

x

− 1)

2

dx

(1)

where kB is the Boltzmann constant, ћ is Plank’s reduced constant, ν is the average sound velocity, θD is the Debye temperature, τc is combined relaxation time, and x is defined as x = ћω/kBT. As mentioned above, both matrix and P1 can be regarded as an identical scattering

ACS Paragon Plus Environment

8

Page 9 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

source. More theoretical details, including the equations for different relaxation times and the relevant parameters for calculation are described in SI. To validate the accuracy of simulation, we compare the calculated and experimental values of κlat for the sample with x = 0.04 (Figure 3b). It is apparent that the simulated κlat agrees quite well with the experimental ones, which reinforces our assumption, i.e. the multiphase plays an important role in scattering heat-carrying phonon. In addition, we also developed the Debye model to examine the individual contribution of these different phases to κlat reduction. Since all the Pb1-xGaxTe samples reach the maximum value of ZT at 823 K, it is crucial to simulate κlat at this temperature (Figure 3d). As for x = 0.01, the experimental value of κlat falls exactly on the simulated red dot line, which takes into account Umklapp process, normal process and point defects (UN + PD) as the dominant phononscattering source, confirming that point defects are the most effective scatters. As for x = 0.02, the deviation between the experimental value of κlat and the calculated line (UN + PD) emerges, demonstrating the intense phonon scattering arising from nanoprecipitates of P2. With the rising of P2 content (x = 0.03), the simulated values (UN + PD + P2) of κlat decrease and match the experimental value of κlat, indicating that the importance of P2 for phonon scattering. Since abundant P3 only appear for the sample with x = 0.04, we just subtracted a fixed value (κlat reduction by P3 at 823 K) from the green dash dot line (UN + PD + P2). Based on the consistency between experimental results and theoretical line (UN + PD + P2 + P3), we reached an unexpected conclusion that the multiple phases (P2 + P3) play a key role in scattering heatcarrying phonon. All of the input parameters based on above theoretical simulation of κlat are listed in Table S2. Such calculated simulations provide a comprehensive account for the intrinsic contribution of each phase and the fundamental role of multiphase.

ACS Paragon Plus Environment

9

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 26

As a result, we achieved a noticeable enhancement of ZT merit through manipulating Ga content. Figure 4a plots the figure of merit ZT as a function of temperature. Because of the increasing contents of P2 and P3, multiphase leads to the lowest values of κlat for the sample with x = 0.04. On the other hand, the carrier mobility µH decreases, leading to the suppression of PF. Therefore, it is not the highest value of ZT for x = 0.04. While for the sample with x = 0.03, it has the highest value of ZT ~1.3 at 823 K. Such enhanced value of ZT driven by carrier tuning and multiphase engineering is much larger than most of previous outstanding results (Figure S5a), especially compared with PbSe-0.5%Ga.37 Additionally, the high performance of Pb0.97Ga0.03Te indicates negligible variation of thermoelectric properties during the experimental heating and cooling cycles, revealing an excellent thermal stability (Figure S6). And also, it is found a good experimental repeatability for high value of ZT, as evidenced by three different samples fabricated independently for the sample with x = 0.03 (Figure S6f). Meanwhile, we should not ignore the importance of average ZTave from the viewpoint of thermoelectric application. The calculated ZTave in the temperature range from 323 to 823 K are shown in Figure 4b, as well as the thermoelectric conversion efficiency η (Figure S5b). It is emphasized that the enhanced PF and reduced κlat for x = 0.03 triggers a large ZTave of ~0.9 in a wide temperature range. Moreover, it is well known that the value of η can be calculated from ZTave, hot temperature Th, and cold side temperature Tc: 44-46

η=

1 + ZTave − 1 Th − Tc Th 1 + ZTave + Tc Th

(2)

Distinctly, from the viewpoint of practical application for thermoelectric, a great η is determined by large ZTave and broad temperature range. However, the previous results for the ntype PbTe, to the best of our knowledge, cannot reach the above requirement due to narrow

ACS Paragon Plus Environment

10

Page 11 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

temperature range (e.g. PbTe-5%InSb29), or relatively low values of ZTave (e.g. PbTe-0.5%Cr,36 Pb0.99Cd0.01Te0.25Se0.7538) as shown in Figure 4c. In comparison, our present work of Pb0.97Ga0.03Te exhibits an appropriately high value of ZTave in a wide temperature region (Th - Tc = 500 K), reflecting its potential application as an ideal candidate of thermoelectric materials. Moreover, the conversion efficiency η for the sample with x = 0.03 for Tc = 323 K and Th = 823 K is ~13%, which exceeds that of many other n-type PbT (T = Te, Se, S) materials (Figure 4d). 4. Conclusions In summary, the introduction of Ga into PbTe matrix enhances the PF through the modulation of carrier concentration. Furthermore, the Ga-doping triggers the emergence of multiple phases (i.e., P1, P2, and P3) responsible for the scattering of heat-carrying phonon, which have been demonstrated by a series of theoretical simulations. In particular, both P2 and P3 dominate the phonon transport, leading to the reduction of κlat. As a result, the ZT value at 823 K reaches ~1.3 for n-type Pb0.97Ga0.03Te, as well as the average ZTave value of ~0.9 and conversion efficiency η of ~13% from 323 to 823 K. These findings open up a novel avenue for manipulating thermoelectric properties in Ga-doped PbTe via synergetic carrier tuning and multiphase engineering.

ASSOCIATED CONTENT Supporting Information. Calculation of lattice thermal conductivity; Rietveld refinement for Ga-doped samples; phase fraction of Ga6PbTe10 in Ga-doped samples; SEM and EDX elemental mapping; thermal

ACS Paragon Plus Environment

11

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 26

diffusivity, heat capacity, Lorenz number, and electronic thermal conductivity for the Pb1-xGaxTe samples; temperature dependence of ZT and calculated conversion efficiency η; stability and reproducibility evaluation for Pb0.97Ga0.03Te; tables with density of all samples for Pb1-xGaxTe and input parameters based on the theoretical simulation of lattice thermal conductivity. AUTHOR INFORMATION Corresponding Author [email protected] (R. Ang), [email protected] (J. Tang) ACKNOWLEDGMENT R. Ang thanks the financial supports from National Natural Science Foundation of China (NSFC) (grant no. 51771126), Youth Foundation of Science & Technology Department of Sichuan Province (grant no. 2016JQ0051), Sichuan University Outstanding Young Scholars Research Funding (grant no. 2015SCU04A20), Sichuan University Talent Introduction Research Funding (grant no. YJ201537), World First-Class University Construction Funding. J. Tang is grateful to the financial support from NSFC (grant no. 11274234). G. Y. Wang thanks the financial supports from Key Research Program of Frontier Sciences, Chinese Academy of Sciences (grant no. QYZDB-SSW-SLH016).

ACS Paragon Plus Environment

12

Page 13 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

REFERENCES (1)

Bell, L. E. Cooling, Heating, Generating Power, and Recovering Waste Heat with Thermoelectric Systems. Science 2008, 321, 1457-1461.

(2)

He, J.; Kanatzidis, M. G.; Dravid, V. P. High Performance Bulk Thermoelectrics via a Panoscopic Approach. Mater. Today 2013, 16, 166-176.

(3)

Zhao, L. D.; Lo, S. H.; Zhang, Y.; Sun, H.; Tan, G.; Uher, C.; Wolverton, C.; Dravid, V. P.; Kanatzidis, M. G. Ultralow Thermal Conductivity and High Thermoelectric Figure of Merit in SnSe Crystals. Nature 2014, 508, 373-377.

(4)

Fu, C.; Bai, S.; Liu, Y.; Tang, Y.; Chen, L.; Zhao, X.; Zhu, T. Realizing High Figure of Merit in Heavy-Band P-Type Half-Heusler Thermoelectric Materials. Nat. Commun. 2015, 6, 8144.

(5)

Snyder, G. J.; Toberer, E. S. Complex Thermoelectric Materials. Nat. Mater. 2008, 7, 105114.

(6)

Wan, S.; Qiu, P.; Huang, X.; Song, Q.; Bai, S.; Shi, X.; Chen, L. Synthesis and Thermoelectric Properties of Charge-Compensated SyPdxCo4–xSb12 Skutterudites. ACS Appl. Mater. Interfaces 2018, 10, 625-634.

(7)

Shuai, J.; Mao, J.; Song, S.; Zhu, Q.; Sun, J.; Wang, Y.; He, R.; Zhou, J.; Chen, G.; Singh, D. J.; Ren, Z. Tuning the Carrier Scattering Mechanism to Effectively Improve the Thermoelectric Properties. Energy Environ. Sci. 2017, 10, 799-807.

ACS Paragon Plus Environment

13

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(8)

Page 14 of 26

Li, Z.; Xiao, C.; Fan, S.; Deng, Y.; Zhang, W.; Ye, B.; Xie, Y. Dual Vacancies: An Effective Strategy Realizing Synergistic Optimization of Thermoelectric Property in BiCuSeO. J. Am. Chem. Soc. 2015, 137, 6587-6593.

(9)

Sun, X.; Guo, Y.; Wu, C.; Xie, Y. The Hydric Effect in Inorganic Nanomaterials for Nanoelectronics and Energy Applications. Adv. Mater. 2015, 27, 3850-3867.

(10) Chen, Z.-G.; Han, G.; Yang, L.; Cheng, L.; Zou, J. Nanostructured Thermoelectric Materials: Current Research and Future Challenge. Prog. Nat. Sci. 2012, 22, 535-549. (11) Tan, G.; Liu, W.; Wang, S.; Yan, Y.; Li, H.; Tang, X.; Uher, C. Rapid Preparation of CeFe4Sb12 Skutterudite by Melt Spinning: Rich Nanostructures and High Thermoelectric Performance. J. Mater. Chem. A 2013, 1, 12657-12668. (12) Liu, H.; Shi, X.; Xu, F.; Zhang, L.; Zhang, W.; Chen, L.; Li, Q.; Uher, C.; Day, T.; Snyder, G. J. Copper Ion Liquid-Like Thermoelectrics. Nat. Mater. 2012, 11, 422-425. (13) Zhao, L.-D.; He, J.; Wu, C. I.; Hogan, T. P.; Zhou, X.; Uher, C.; Dravid, V. P.; Kanatzidis, M. G. Thermoelectrics with Earth Abundant Elements: High Performance P-Type PbS Nanostructured with SrS and CaS. J. Am. Chem. Soc. 2012, 134, 7902-7912. (14) Pei, Y.; LaLonde, A. D.; Heinz, N. A.; Shi, X.; Iwanaga, S.; Wang, H.; Chen, L.; Snyder, G. J. Stabilizing the Optimal Carrier Concentration for High Thermoelectric Efficiency. Adv. Mater. 2011, 23, 5674-5678. (15) Pei, Y.; Shi, X.; LaLonde, A.; Wang, H.; Chen, L.; Snyder, G. J. Convergence of Electronic Bands for High Performance Bulk Thermoelectrics. Nature 2011, 473, 66-69.

ACS Paragon Plus Environment

14

Page 15 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(16) Zhao, L.-D.; Hao, S.; Lo, S. H.; Wu, C. I.; Zhou, X.; Lee, Y.; Li, H.; Biswas, K.; Hogan, T. P.; Uher, C.; Wolverton, C.; Dravid, V. P.; Kanatzidis, M. G. High Thermoelectric Performance via Hierarchical Compositionally Alloyed Nanostructures. J. Am. Chem. Soc. 2013, 135, 7364-7370. (17) He, J.; Zhao, L.-D.; Zheng, J. C.; Doak, J. W.; Wu, H.; Wang, H. Q.; Lee, Y.; Wolverton, C.; Kanatzidis, M. G.; Dravid, V. P. Role of Sodium Doping in Lead Chalcogenide Thermoelectrics. J. Am. Chem. Soc. 2013, 135, 4624-4627. (18) Biswas, K.; He, J.; Zhang, Q.; Wang, G.; Uher, C.; Dravid, V. P.; Kanatzidis, M. G. Strained Endotaxial Nanostructures with High Thermoelectric Figure of Merit. Nat. Chem. 2011, 3, 160-166. (19) Wu, D.; Zhao, L.-D.; Tong, X.; Li, W.; Wu, L.; Tan, Q.; Pei, Y.; Huang, L.; Li, J.-F.; Zhu, Y.; Kanatzidis, M. G.; He, J. Superior Thermoelectric Performance in PbTe–PbS PseudoBinary: Extremely Low Thermal Conductivity and Modulated Carrier Concentration. Energy Environ. Sci. 2015, 8, 2056-2068. (20) Chen, Z.; Jian, Z.; Li, W.; Chang, Y.; Ge, B.; Hanus, R.; Yang, J.; Chen, Y.; Huang, M.; Snyder, G. J.; Pei, Y. Lattice Dislocations Enhancing Thermoelectric PbTe in Addition to Band Convergence. Adv. Mater. 2017, 29, 1606768. (21) Popescu, A.; Haney, P. M. Interface Scattering in Polycrystalline Thermoelectrics. J. Appl. Phys. 2014, 115, 123701.

ACS Paragon Plus Environment

15

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 26

(22) Biswas, K.; He, J.; Blum, I. D.; Wu, C. I.; Hogan, T. P.; Seidman, D. N.; Dravid, V. P.; Kanatzidis, M. G. High-Performance Bulk Thermoelectrics with All-Scale Hierarchical Architectures. Nature 2012, 489, 414-418. (23) Girard, S. N.; He, J.; Zhou, X.; Shoemaker, D.; Jaworski, C. M.; Uher, C.; Dravid, V. P.; Heremans, J. P.; Kanatzidis, M. G. High Performance Na-Doped PbTe-PbS Thermoelectric Materials: Electronic Density of States Modification and Shape-Controlled Nanostructures. J. Am. Chem. Soc. 2011, 133, 16588-16597. (24) Wu, H. J.; Zhao, L. D.; Zheng, F. S.; Wu, D.; Pei, Y. L.; Tong, X.; Kanatzidis, M. G.; He, J. Q. Broad Temperature Plateau for Thermoelectric Figure of Merit ZT > 2 in PhaseSeparated PbTe0.7S0.3. Nat. Commun. 2014, 5, 4515. (25) Tan, G.; Shi, F.; Hao, S.; Zhao, L.-D.; Chi, H.; Zhang, X.; Uher, C.; Wolverton, C.; Dravid, V. P.; Kanatzidis, M. G. Non-Equilibrium Processing Leads to Record High Thermoelectric Figure of Merit in PbTe–SrTe. Nat. Commun. 2016, 7, 12167. (26) Jood, P.; Ohta, M.; Kunii, M.; Hu, X.; Nishiate, H.; Yamamoto, A.; Kanatzidis, M. G. Enhanced Average Thermoelectric Figure of Merit of N-Type PbTe1−xIx-MgTe. J. Mater. Chem. C 2015, 3, 10401-10408. (27) Girard, S. N.; Chasapis, T. C.; He, J.; Zhou, X.; Hatzikraniotis, E.; Uher, C.; Paraskevopoulos, K. M.; Dravid, V. P.; Kanatzidis, M. G. PbTe–PbSnS2 Thermoelectric Composites: Low Lattice Thermal Conductivity from Large Microstructures. Energy Environ. Sci. 2012, 5, 8716-8725.

ACS Paragon Plus Environment

16

Page 17 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(28) Tan, G.; Stoumpos, C. C.; Wang, S.; Bailey, T. P.; Zhao, L.-D.; Uher, C.; Kanatzidis, M. G. Subtle Roles of Sb and S in Regulating the Thermoelectric Properties of N-Type PbTe to High Performance. Adv. Energy Mater. 2017, 7, 1700099. (29) Zhang, J.; Wu, D.; He, D.; Feng, D.; Yin, M.; Qin, X.; He, J. Extraordinary Thermoelectric Performance Realized in N-Type PbTe through Multiphase Nanostructure Engineering. Adv. Mater. 2017, 29, 1703148. (30) Xuan, X. C. Optimum Design of a Thermoelectric Device. Semicond. Sci. Technol. 2002, 17, 114. (31) Pei, Y.; May, A. F.; Snyder, G. J. Self-Tuning the Carrier Concentration of PbTe/Ag2Te Composites with Excess Ag for High Thermoelectric Performance. Adv. Energy Mater. 2011, 1, 291-296. (32) Heremans, J. P.; Jovovic, V.; Toberer, E. S.; Saramat, A.; Kurosaki, K.; Charoenphakdee, A.; Yamanaka, S.; Snyder, G. J. Enhancement of Thermoelectric Efficiency in PbTe by Distortion of the Electronic Density of States. Science 2008, 321, 554-557. (33) Fu, L.; Yin, M.; Wu, D.; Li, W.; Feng, D.; Huang, L.; He, J. Large Enhancement of Thermoelectric Properties in N-Type PbTe via Dual-Site Point Defects. Energy Environ. Sci. 2017, 10, 2030-2040. (34) Xiao, Y.; Wu, H.; Li, W.; Yin, M.; Pei, Y.; Zhang, Y.; Fu, L.; Chen, Y.; Pennycook, S. J.; Huang, L.; He, J.; Zhao, L.-D. Remarkable Roles of Cu to Synergistically Optimize Phonon and Carrier Transport in N-Type PbTe-Cu2Te. J. Am. Chem. Soc. 2017, 139, 18732-18738.

ACS Paragon Plus Environment

17

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 26

(35) Pei, Y.; Lensch-Falk, J.; Toberer, E. S.; Medlin, D. L.; Snyder, G. J. High Thermoelectric Performance in PbTe due to Large Nanoscale Ag2Te Precipitates and La Doping. Adv. Funct. Mater. 2011, 21, 241-249. (36) Zhang, Q.; Chere, E. K.; McEnaney, K.; Yao, M.; Cao, F.; Ni, Y.; Chen, S.; Opeil, C.; Chen, G.; Ren, Z. Enhancement of Thermoelectric Performance of N-Type PbSe by Cr Doping with Optimized Carrier Concentration. Adv. Energy Mater. 2015, 5, 1401977. (37) Zhang, Q.; Cao, F.; Lukas, K.; Liu, W.; Esfarjani, K.; Opeil, C.; Broido, D.; Parker, D.; Singh, D. J.; Chen, G.; Ren, Z. Study of the Thermoelectric Properties of Lead Selenide Doped with Boron, Gallium, Indium, or Thallium. J. Am. Chem. Soc. 2012, 134, 1773117738. (38) Zhang, Q.; Chere, E. K.; Wang, Y.; Kim, H. S.; He, R.; Cao, F.; Dahal, K.; Broido, D.; Chen, G.; Ren, Z. High Thermoelectric Performance of N-Type PbTe1−ySy due to Deep Lying States Induced by Indium Doping and Spinodal Decomposition. Nano Energy 2016, 22, 572-582. (39) Pei, Y.; Gibbs, Z. M.; Gloskovskii, A.; Balke, B.; Zeier, W. G.; Snyder, G. J. Optimum Carrier Concentration in N-Type PbTe Thermoelectrics. Adv. Energy Mater. 2014, 4, 1400486. (40) Ravich, Y. I.; Efimova, B. A.; Smirnov, I. A. Semiconducting Lead Chalcogenides; Plenum Press: New York, 1970. (41) Klemens, P. G. The Scattering of Low-Frequency Lattice Waves by Static Imperfections. Proc. Phys. Soc. A 1955, 68, 1113-1128.

ACS Paragon Plus Environment

18

Page 19 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(42) Klemens, P. G. Thermal Resistance due to Point Defects at High Temperatures. Phys. Rev. 1960, 119, 507-509. (43) Callaway, J.; von Baeyer, H. C. Effect of Point Imperfections on Lattice Thermal Conductivity. Phys. Rev. 1960, 120, 1149-1154. (44) Goldsmid, H. J. Introduction to Thermoelectricity; Springer: Heidelberg, 2009. (45) Yang, Y.; Pradel, K. C.; Jing, Q.; Wu, J. M.; Zhang, F.; Zhou, Y.; Zhang, Y.; Wang, Z. L. Thermoelectric Nanogenerators Based on Single Sb-Doped ZnO Micro/Nanobelts. ACS Nano 2012, 6, 6984-6989. (46) Kraemer, D.; Jie, Q.; McEnaney, K.; Cao, F.; Liu, W.; Weinstein, L. A.; Loomis, J.; Ren, Z.; Chen, G. Concentrating Solar Thermoelectric Generators with a Peak Efficiency of 7.4%. Nat. Energy 2016, 1, 16153. (47) Ibanez, M.; Korkosz, R. J.; Luo, Z.; Riba, P.; Cadavid, D.; Ortega, S.; Cabot, A.; Kanatzidis, M. G. Electron Doping in Bottom-Up Engineered Thermoelectric Nanomaterials through HCl-Mediated Ligand Displacement. J. Am. Chem. Soc. 2015, 137, 4046-4049. (48) Chere, E. K.; Zhang, Q.; McEnaney, K.; Yao, M.; Cao, F.; Sun, J.; Chen, S.; Opeil, C.; Chen, G.; Ren, Z. Enhancement of Thermoelectric Performance in N-Type PbTe1−ySey by Doping Cr and Tuning Te:Se Ratio. Nano Energy 2015, 13, 355-367.

ACS Paragon Plus Environment

19

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 26

Figure captions:

Figure 1. (a) Powder XRD patterns for Pb1-xGaxTe (x = 0.01, 0.02, 0.03, 0.04); (b) carrier concentration nH and mobility µH at 300 K. Temperature dependence of (c) electrical conductivity σ; (d) Seebeck coefficient S; (e) power factor PF; (f) Seebeck coefficient at room temperature as a function of carrier concentration for this work and refs 35, 38. The solid curve is the theoretical Pisarenko plot for n-type PbTe with effective mass of electrons of m* = 0.25 me.

Figure 2. (a) Surface morphology and corresponding EDS mapping for x = 0.01; (b) surface morphology for x = 0.02, along with the EDS spectrum and the corresponding quantitative results of matrix and P1, respectively; (c) surface morphology and corresponding EDS mapping for x = 0.03; (d) surface morphology for x = 0.04, along with the EDS spectrum and the corresponding quantitative results of P2 and P3, respectively.

Figure 3. Thermal transport properties and simulation of Pb1-xGaxTe. Temperature dependence of (a) total thermal conductivity κtot; and (b) lattice thermal conductivity κlat; the orange dot line represents the simulated κlat for x = 0.04 based on Debye approximation; (c) κlat at room temperature as a function of Ga content, and the solid line stands for κlat in the framework of Klemens-Debye model; (d) Ga content dependence of κlat at 823 K, where the dashed curves express the theoretical predictions based on a Debye approximation with various types of phonon scattering, the solid dots stand for the experimental results.

ACS Paragon Plus Environment

20

Page 21 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Figure 4. Temperature dependence of (a) ZT; and (b) average ZTave for Pb1-xGaxTe, compared with ref 38; (c) comparison of ZTave with previous n-type lead chalcogenides refs 14, 26, 28, 29, 36-38, 47, 48; (d) corresponding ZTave values and calculated conversion efficiency η in a wide temperature range from 323 to 823 K.

ACS Paragon Plus Environment

21

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 26

Figure 1

ACS Paragon Plus Environment

22

Page 23 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Figure 2

ACS Paragon Plus Environment

23

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 26

Figure 3

ACS Paragon Plus Environment

24

Page 25 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Figure 4

ACS Paragon Plus Environment

25

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 26

Table of Contents

ACS Paragon Plus Environment

26