Rationally Designing High-Performance Bulk ... - ACS Publications

Aug 31, 2016 - ABSTRACT: There has been a renaissance of interest in exploring highly efficient ... high-performance bulk thermoelectric materials. We...
2 downloads 0 Views 8MB Size
Review pubs.acs.org/CR

Rationally Designing High-Performance Bulk Thermoelectric Materials Gangjian Tan,† Li-Dong Zhao,*,‡ and Mercouri G. Kanatzidis*,† †

Department of Chemistry, Northwestern University, Evanston, Illinois 60208, United States School of Materials Science and Engineering, Beihang University, Beijing 100191, China



ABSTRACT: There has been a renaissance of interest in exploring highly efficient thermoelectric materials as a possible route to address the worldwide energy generation, utilization, and management. This review describes the recent advances in designing high-performance bulk thermoelectric materials. We begin with the fundamental stratagem of achieving the greatest thermoelectric figure of merit ZT of a given material by carrier concentration engineering, including Fermi level regulation and optimum carrier density stabilization. We proceed to discuss ways of maximizing ZT at a constant doping level, such as increase of band degeneracy (crystal structure symmetry, band convergence), enhancement of band effective mass (resonant levels, band flattening), improvement of carrier mobility (modulation doping, texturing), and decrease of lattice thermal conductivity (synergistic alloying, second-phase nanostructuring, mesostructuring, and all-length-scale hierarchical architectures). We then highlight the decoupling of the electron and phonon transport through coherent interface, matrix/precipitate electronic bands alignment, and compositionally alloyed nanostructures. Finally, recent discoveries of new compounds with intrinsically low thermal conductivity are summarized, where SnSe, BiCuSeO, MgAgSb, complex copper and bismuth chalcogenides, pnicogen-group chalcogenides with lone-pair electrons, and tetrahedrites are given particular emphasis. Future possible strategies for further enhancing ZT are considered at the end of this review.

CONTENTS 1. Introduction 2. Basic Stratagem of Optimizing ZT: Carrier Concentration Engineering 3. Means of Increasing Maximum ZT 3.1. Enhancement of Carrier Effective Mass m* 3.1.1. Increasing the Number of Band Extrema NV 3.1.2. Increase of Carrier Effective Mass mb* 3.2. Modulation Doping and Carrier Mobility Improvement 3.3. Reducing Lattice Thermal Conductivity 3.3.1. Atomic Scale: Synergistic Alloying 3.2.2. Nanoscale: Nanostructuring with Second Phases 3.3.3. Mesoscale Structuring 3.3.4. All-Scale Hierarchical Architectures 4. Decoupling of Electron and Phonon Transport 4.1. Strained Endotaxial Nanostructuring 4.2. Matrix/Precipitate Band Alignment 4.3. Compositionally Alloyed Nanostructures 5. Discoveries of New Thermoelectric Materials with Intrinsically Low Thermal Conductivity 5.1. Layered SnSe 5.2. BiCuSeO Oxyselenides 5.3. Half-Heusler MgAgSb 5.4. Copper Chalcogenides 5.5. Complex Bismuth Chalcogenides © 2016 American Chemical Society

5.6. Chalcoantimonates with Lone-Pair Electrons 5.7. Tetrahedrites 6. Summary and Outlook Author Information Corresponding Authors Notes Biographies Acknowledgments References

12123 12124 12126 12126 12126 12127 12128 12129 12129

12139 12141 12141 12142 12142 12142 12142 12142 12142

1. INTRODUCTION With more than two-thirds of the worldwide utilized energy being wasted as heat releasing to the atmosphere in vain, it is of great economic and environmental benefit to capture this largely untapped waste heat to generate emission-free renewable power. Thermoelectric materials that are capable of directly and reversibly converting heat into electricity have thus drawn increasing attention over the past several decades. The key of promoting this promising power generation technology into massive application is to increase the low conversion efficiency of current thermoelectric materials, which is technically evaluated by the dimensionless figure of merit ZT = S2σT/κ. A good thermoelectric material should have a large Seebeck coefficient (S) that is usually present in semi-

12130 12130 12131 12131 12131 12132 12133 12133 12134 12136 12137 12138 12139

Received: April 21, 2016 Published: August 31, 2016 12123

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

crystalline SnSe, a simple binary compound containing no toxic or scarce elements, displays an excitingly record high ZT of ∼2.6 at 923 K along the b axis of the room-temperature orthorhombic unit cell.40 This can be achieved because of the giant anharmonic and anisotropic bonding of SnSe which gives rise to an extremely low thermal conductivity.40 Thus, the search for new compounds with intriguing physical properties is a necessary activity to achieve high thermoelectric performance in addition to optimizing the known systems. Although in the past decade breakthroughs have been constantly made in the thermoelectric field, systematic and indepth analysis and new perspectives on the applied approaches and concepts are valuable. There were previously a few nice introductory reviews on thermoelectric materials, which could be a nice starting point for the readers who are interested in thermoelectrics. Typical examples are summarized by Wood,45 Disalvo,46 Snyder,47 and Kanatzidis.48 Interested readers are also encouraged to refer to some outstanding review articles that focus on some topics, e.g., nanostructuring,20,49,50 band structure engineering,26,51 or nanopolycrystalline composites,52,53 and some others that focus on several specific materials, such as clathrates,54 Zintl phase compounds,55 oxides,56,57 organic materials,58,59 and others.47,60−62 A recent review discusses issues of charge and phonon transport in thermoelectric compounds from the chemical bonding point of view.63 This review aims to summarize the most recent state-ofthe-art approaches of designing high-performance bulk thermoelectric materials which are most likely to enable the devices for large-scale power generation application. It covers all aspects, from improving the electrical properties in the manner of band structure and microstructure manipulations to decreasing the thermal conductivity by composition and microstructure designs and to decoupling the electron and phonon transport through matrix/precipitates band alignment and compositionally alloyed nanostructures. The advances in finding new compounds with intrinsically low thermal conductivity are also highlighted. Future potential strategies universally applicable to all compounds aimed at enhancing the performance are discussed. They can be a comprehensive guide to designing and exploring advanced bulk thermoelectric materials.

conductors, high electrical conductivity (σ) just like metals, and in the meantime poor thermal conductivity (κ) as in glasses. It is knotty to combine all of these features in a single material, enabling its high performance. Experimentally, to improve ZT, one can increase the numerator S2σ (also known as power factor) and/or decrease the denominator κ (thermal conductivity) by rational band structure and microstructure designs.1−7 Since the historic discoveries of thermoelectric phenomena by Seebeck and Peltier in the first half of the 19th century, many materials have been explored and considered useful to generate thermoelectricity, including metals,8−11 ceramics,12−14 and ultimately semiconductors.15−17 It needs to be emphasized that almost all these materials were initially obtained empirically through countless attempts based on personal experience. Therefore, the growth of the thermoelectric community was very slow in the beginning until the 1950s, the decade when the basic science of thermoelectric effects was established and heavily doped semiconductors were widely recognized as excellent thermoelectric materials.15−18 Since then the field of thermoelectrics advanced rapidly, especially in the past two decades, relying on the development of new concepts, better intuitive thinking and theoretical ideas relating to size effects,1−5,19−22 and band structuring engineering.3,23−28 Figure 1 shows the primary milestones achieved for ZT values of bulk materials over the last two decades as a function

2. BASIC STRATAGEM OF OPTIMIZING ZT: CARRIER CONCENTRATION ENGINEERING The fundamental challenge of designing high-ZT thermoelectric materials stems from the strong correlation of S, σ, and κ through carrier concentration n which can be adjusted by controlling the doping level. For a degenerated semiconductor with parabolic band dispersion, assuming that dopant does not change the scattering or band structure significantly, S is given by64

Figure 1. Current state-of-the-art bulk thermoelectric materials: the thermoelectric figure-of-merit ZT as a function of temperature and year illustrating important milestones. Green cylinders represent the ptype materials, while red cylinders represent the n-type ones. CsBi4Te6,29 BiSbTe,2 AgPbmSbTe2+m,1 Ba8Ga16Ge30,30 Tl−PbTe,23 In4Se3,31 CuxBiTeSe,32 (BaLaYb)xCo4Sb12,33 MgAgSb,34 BiSbTe+Te arrays,22 PbTe−SrTe,4 DDyFe3CoSb12,35 Mg2Si0.3Sn0.7,24 Cu2−xSe,36 Na2+xGa2+xSn4−x,37 GeTe−Bi2Te3,38 Pb(Te,Se,S),39 PbS−CdS,27 SnSe,40 Cu2−xS,41 PbTe−PbS,42 SnTe−CdS,43 and Pb1−xSbxSe.44

of both year and temperature. Great achievements have been achieved with both n- and p-type materials with the greatest ZTs around or above 2, especially since the year 2010.4,39,40,42 Traditional thermoelectric materials, namely, bismuth telluride2,22 and lead chalcogenides,1,3,4,23,25,27,39,42,44 remain the top performing ones in their respective working temperature region and attract strong attention. On the other hand, because of concerns of the scarceness of tellurium and the perceived toxicity of lead in lead chalcogenides, researchers also explore the possibility of Te- and Pb-free compounds, among which MgAgSb,34 skutterudites,33,35 and copper and tin chalcogenides36,40,41,43 appear as most promising. Particularly, single-

S=

⎛ π ⎞2/3 m*T ⎜ ⎟ ⎝ 3n ⎠ 3eh

8π 2kB2 2

(1)

where kB is the Boltzmann constant, e the electron charge, h the Plank constant, and m* the density of states effective mass of carriers. Apparently, high S is usually found in low-n semiconductors or insulators, and with increasing n, S drops rapidly, Figure 2a. On the contrary, to ensure a large σ, n should be as high as possible, see eq 2 and Figure 2a σ = neμ 12124

(2) DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

orange line in Figure 2b and 2c. The consequence of this is that the theoretically maximum performance cannot be fully realized at every working temperature. One way to deal with this issue is to use functionally graded doping67−69 by integrating two or multiple segments with dissimilar n, see green lines in Figures 2b and 3c. One common approach developed in the late 2000s

Figure 2. (a) Schematic diagram showing how ZT and its related parameters (electrical conductivity σ, Seebeck coefficient S, power factor S2σ, electrical thermal conductivity κele, lattice thermal conductivity κlat, and total thermal conductivity κ) change as a function of carrier concentration n. (b) Strategies for stabilizing the optimal carrier concentration (n*, denoted by the red line, usually shows a (T)3/2 dependence66). For most conventional dopants, the resultant carrier concentration (denoted by orange dotted line) is almost temperature independent, functionally graded doping67−69 (green lines) by use of samples with dissimilar n, and use of temperature-dependent-solubility dopant,71,72 namely, n has a strong temperature dependence (purple line). (c) Comparison of different doping methods: (1) conventional doping, (2) graded doping, and (3) T-dependent doping. (d) Enhancement of ZT values over a broad temperature range through stabilizing n* in comparison to the conventional doping approach.69,71,75

Figure 3. (a) Band convergence effect by forming solid solution A1−xBxC between AC and BC: evolution of band structure with increasing doping fraction; Eg and ΔEL−Σ represent the band gap and energy separation between L and Σ bands, respectively. (b) Seebeck coefficient, (c) power factor, and (d) ZT values as a function of both temperature and doping fraction of Sn1−xMnxTe.89 The temperaturedependent Seebeck coefficient of SnTe-3%CdTe92 is included in b for comparison. (Inset of c) Reduction of lattice thermal conductivity with increasing Mn doping concentration.89

where μ is the carrier mobility. Only within the proper carrier concentration range (e.g., 1019−1020 cm−3 for most semiconductors47) can the power factors be maximized by balancing S and σ. The thermal conductivity κ comprises two major components: the electronic contribution κele and the lattice contribution κlat (among all ZT-related physical parameters, κlat is the only one with little dependence on n (Figure 2a) and potentially can be minimized to the amorphous limit through independent crystal structure and/or microstructure design4,55) κ = κele + κlat (3)

of preparing such a graded material is the so-called spark plasma sintering or hot pressing a compacted stack of powder layers, each of which has different carrier concentrations.67,69 However, after extended duration in service the initial carrier concentration gradient in the graded material may fade or vanish due to the diffusion-induced homogenization effect, thereby deteriorating the conversion efficiency.69,70 Alternatively, by making use of the temperature-dependent solubility71,72 of some specific dopants (Figure 2c), one can create a gradient of n within a single material controlled only by temperature and its gradient (see purple line in Figure 2b and 2c). A well-known example is that of Cu, Ag, and excess Pb, which have negligible solubility within PbTe around room temperature but are highly soluble at elevated temperatures.73,74 This temperature-dependent doping is reversible in the heating−cooling process and gets rid of the diffusion problem with graded doping, making it a better candidate for actual application.4,71 Figure 2d compares the temperature-dependent ZT values of n-type PbTe using different doping methods. The green and orange dotted lines represent A and B materials of conventional n-type PbTe with n = 5 × 1018 and 3 × 1018 cm−3, respectively, and the dotted purple line is for the graded doping by integrating the two components.69 A has small n and thus high ZT at low temperatures, while B with larger n has significantly higher ZT at high temperatures. Graded doping makes the ZT fall between A and B but can achieve a much larger average ZT in the entire working temperature range. The solid black line

and κele is influenced by n through the Wiedemann−Franz relation65

κele = LσT = LneμT

(4)

where L is the Lorenz number. Equation 4 links the electrical to thermal transport and makes the optimization of ZT even tougher: there is a trade-off between high σ and low κ (κele), Figure 2a. Nonetheless, as shown in Figure 2a, one can always achieve a maximum ZT at a given temperature by finely tuning n to an optimal level using an aliovalent substitution approach. However, the optimal carrier concentration (n*) of a thermoelectric semiconductor is not constant at all temperatures and usually increases rapidly with rising temperature, roughly obeying the power law of (T)3/2,66 as shown by the red line in Figure 2b. In the case of conventional doping, the carrier concentrations are generally irrelevant to temperature, see the 12125

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

shows the theoretically optimal ZT for n-type PbTe.71 Conventional doping approaches, such as by heavy La doping75 in PbTe (solid navy line), can make the maximum ZT approach the optimal value around 750 K but cannot achieve high performance at all temperatures. The temperature-dependent solubility limit of Ag in the PbTe+Ag2Te composite (Tdependent doping, solid green line) enables an increase in n with increasing T and pushes ZT values close to the optimal levels over the entire temperature range. The T-dependent doping approach has also been successfully applied in many other systems including p-type Na-doped PbTe.4,42,72,76 It was demonstrated that the Na-rich precipitates residing at the grain boundaries at low temperature are redissolved into the PbTe matrix at T > 600 K.4,42 This gives rise to the enhancement of hole density, contributing to the increase of electrical conductivity and power factor and suppression of bipolar conduction at elevated temperature achieving superior performance. The modifications to the conventional doping method to include temperature-dependent doping behavior are helpful to acquire larger average ZT values which are especially important in technological applications.

noncubic lattice of either of the chalcopyrite compounds can be distorted to be cubic-like (c = 2a) on average. This reconstructed pseudocubic structure was perceived to be responsible for the lower energy-splitting parameter (or larger NV) and higher performance obtained in the solid solutions. More work is needed, however, in order to better understand just how much of the net ZT enhancement comes from the increase of NV since lattice thermal conductivity reduction through forming the solid solution is also contributing.82 Another attractive way to increase the effective NV is to converge different bands in the Brillouin zone within a few kBT in energy of each other. For example, both experimental and theoretical investigations suggest the presence of two valence bands in lead chalcogenides: one at the L point with NV = 4 and the other at the Σ point of the Brillouin zone with a much larger NV of 12.81 Normally, the two bands are separated by an energy ΔEL−Σ that is sufficiently large to make the Σ band irreverent to charge transport, Figure 3a. However, upon proper alloying with specific elements (for example, Mg,83 Cd,84 Sr,85 or Mn86 in PbTe; Sr in PbSe87), ΔEL−Σ can be significantly lowered, and the two valence bands become converged gradually. In that case, the effective NV becomes larger than 12, and higher performance can be achieved.25,26,87 Likewise, notable conduction band convergence was demonstrated to be effective in enhancing the n-type thermoelectric performance of Mg2Si by alloying with Sn at the Si site.24 We note that the mobility of carriers is nominally unaffected by NV, but there may be some reduction because of intervalley scattering.25 Recently, such a band convergence concept has also been extended to p-type SnTe,43,88−92 which is an analogue of lead chalcogenides but is supposed to be more environmental friendly. Pristine SnTe is not a good thermoelectric material because of high lattice thermal conductivity and the unavailability of its Σ valence band. Specifically, ΔEL−Σ is ∼0.35 eV for SnTe at 300 K,93 much larger than in PbTe (∼0.17 eV)94 or PbSe (∼0.24 eV).95 This means that a considerably larger amount of dopant is needed to produce the same band convergence effect as in PbTe or PbSe, because the doping-induced band convergence is dependent on the amount of dopant,96 Figure 3a. Cd43,92 and Hg90 are good alloying elements from the perspective of decreasing ΔEL−Σ in SnTe. However, because of their low solubility limit (∼3 mol %),90,92 the resultant ΔEL−Σ is still too high to reach the level of PbTe. For example, 3 mol % CdTe-alloyed SnTe92 has the highest Seebeck coefficient of ∼200 μV/K at ∼900 K, which is significantly larger than that of pristine SnTe but is far below that of p-type PbTe83 at the same hole concentration, Figure 3b. Most recently, Mn was explored as a valence band convergence producer in SnTe with a high solubility of >9 mol %.89,97,98 This gives rise to a much lower ΔEL−Σ that is comparable to that of PbTe.89 Correspondingly, the Seebeck coefficients and power factors of Sn1−xMnxTe are remarkably improved with increasing Mn doping fraction,89 Figure 3b and 3c. Notably, the 12 mol % Mn-doped SnTe89 displays the highest Seebeck coefficients (∼80 and ∼230 μV/K at 300 and 900 K, respectively), similar to PbTe83 with comparable hole densities, Figure 3b. Moreover, the higher doping fraction of Mn also facilitates the decrease of lattice thermal conductivity as a result of stronger point defect scattering,89 inset of Figure 3c. Concurrently, ZTs of SnTe are steadily improved upon Mn doping with a highest value approaching 1.3 at 900 K,89 a 120% increase over pristine SnTe. The above examples demonstrate

3. MEANS OF INCREASING MAXIMUM ZT Maximum performance of a given thermoelectric material is generally achieved by carrier concentration optimization, as we elaborated above. However, to increase the maximum ZT to much higher levels, one needs to carefully tailor the electronic structures and microstructures. Previous studies indicate that the dimensionless merit factor (β) of thermoelectric materials can be characterized by77,78 β=9

5/2 U ⎛ T ⎞ ⎜ ⎟ κlat ⎝ 300 ⎠

(5) 2

−1 −1

where U is the weighted mobility (in cm V s , a measure of the electrical properties) that can be written as79

U = (m*)3/2 μ

(6)

In eq 6, m* is the effective mass of carriers (in me, me is the free electron mass) and relates to the band effective mass mb* through the number of equivalent degenerated valleys of the band structure (NV)26,80 m* = (NV )2/3 mb*

(7)

Thus, to increase the maximum ZT, one needs to maximize NV, mb*, and μ and minimize κlat concurrently to achieve higher S and σ but lower κ at a constant n (according to eqs 1−4). Below we discuss the possible approaches. 3.1. Enhancement of Carrier Effective Mass m*

3.1.1. Increasing the Number of Band Extrema NV. Nv is closely related to the crystal structure symmetry and can be large when the crystal structures are highly symmetric. For example, PbTe with its cubic rocksalt structure has NV of 4 and 12 at the L and Σ points of the valence band, respectively.46,81 Other good thermoelectric materials, such as SiGe alloys, CoSb3, Mg2Si, and half-Heusler alloy, also feature NV > 1.26 Currently, strategies on how to increase the asymmetryrelated Nv are in their infancy. One illuminating attempt was recently carried out on tetragonal chalcopyrite compounds.82 By mixing two types of chalcopyrite compounds (c/2a > 1 for one but a) in a proper ratio to form solid solutions, the 12126

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

1.8 × 1019 cm−3), the La-doped PbTe displays lower ZT value than that of I-doped samples.99 This is because despite the former’s larger mb* and higher Seebeck coefficient, there is a decline of μ leading to a net decrease of power factor. Nonetheless, at a higher doping level of n, La-doped PbTe samples can have comparable or even better power factors than I-doped PbTe samples, as demonstrated by both experimental results and theoretical predictions.99 Indeed, to make best use of the benefits of larger mb* through band flattening for power factor enhancement it generally requires a much higher optimal n to reach the same Fermi level as the one with smaller mb*. Figure 4b plots the power factor as a function of carrier concentration at 800 K for several state-of-the-art thermoelectric materials with different effective masses.107 It is clearly evident that the maximum power factor values of these heavy bands compounds (e.g., ntype CoSb3) are significantly higher than those with small mb* (for example, n-type PbTe). Moreover, the optimal n increases systematically with increasing mb*. Band flattening can increase not only mb* but also the band gap Eg, according to the Kane-band model.108 This model is widely accepted to describe the band dispersion of most excellent thermoelectric materials (for example, CoSb3109 and PbTe71,75)

the feasibility of turning poor thermoelectric materials into excellent ones by energetically converging multiple bands. 3.1.2. Increase of Carrier Effective Mass mb*. Physically, the effective mass mb* in a band is related to the curvature of the bands and can be modified by distorting the band. This can increase the density of states in the vicinity of Fermi level.23,51 Two commonly adopted approaches of enhancing mb* include band flattening99,100 and resonant levels23,91,101,102 through chemical doping, Figure 4a. It should be emphasized that in

⎛ ℏ2k 2 E⎞ = E⎜⎜1 + ⎟⎟ Eg ⎠ 2mb* ⎝

(9)

In eq 9, ℏ is the reduced Plank’s constant, k the crystal momentum, and E the energy of electron states. Intuitively, this additional band gap enlargement effect by band flattening can help to suppress the detrimental effect of bipolar diffusion on the thermoelectric properties.83,92 The concept of resonant levels was first proposed in metals in 1950s.110,111 Resonance levels originate from the coupling between electrons of a dilute impurity with those of the conduction or the valence band of the host solid. This creates excess density of states near the resonant energy and can give rise to improved mb* and enhanced Seebeck coefficient if the Fermi level can be brought close in energy.51 However, resonant levels are generally effective at low temperature, and their effect tends to diminish at high temperatures where the relaxation time of acoustic phonon scattering is much shorter than that of resonant impurity scattering.43,51 The most well-known examples of resonant levels are associated with the group-III impurities in IV−VI semiconductors. For example, Al was reported to form resonant levels in the conduction band of PbSe,102 while In43,91 and Tl23 create resonant levels in the valence bands of SnTe and PbTe, respectively. In addition to IV−VI semiconductors, resonant levels are also observed in Sn-doped Bi2Te3.101 Recently, the synergy of resonant levels and band convergence to strongly increase performance has been demonstrated in In- and Cd-codoped SnTe,43 Figure 4c and 4d. In forms resonant levels in the valence band of SnTe and therefore strongly improved the Seebeck coefficient and power factors of SnTe below 500 K. However, with increasing temperature, this enhancement becomes less significant because of the diminished resonant scattering at elevated temperatures. On the other hand, Cd substitution92 creates valence band convergence in SnTe that works most efficiently at higher temperatures where two valence bands approach each other, enabling a clear enhancement of the Seebeck coefficient, Figure

Figure 4. (a) Schematic representation of the density of states of a single valence band (blue line) contrasted to that of band flattening (purple line) and the introduction of resonant states (red line). (b) Power factors as a function of carrier concentrations for various of compounds (n-type PbTe, ZrNiSn, CoSb3 and p-type NbFeSb) with different effective masses (m*) at 800 K.107 (c) Seebeck coefficient, and (d) power factor as a function of temperature for SnTe with resonance levels by In doping and/or valence band convergence by alloying with Cd.43

most cases the increase of mb* will result in the decrease of carrier mobility and through it the electrical conductivity, according to the relationship46 eτ μ= ml* (8) where ml* is the inertial mass along the transport direction (and equal to mb* for an isotropic band) and τ is the scattering time which decreases with increasing mb* so long as the carriers are scattered predominately by phonons.26 Band flattening effects can be achieved through the adoption of dopants containing highly localized orbitals which can decrease orbital overlap. For instance, both theoretical and experimental studies suggest that La103−106 or Nd103,106 substitution for Sr significantly flattens the conduction band of SrTiO3 perovskites for largely improved density of states effective mass of electrons, thereby leading to higher power factors and thermoelectric figure of merits. Likewise, in Ladoped PbTe,99,100 because of the hybridization between La f states (or d states if f states are fully localized) and Pb p states, the conduction band at the L point can be affected, contributing to the increase of mb* by band flattening.100 In comparison, iodine-doped n-type PbTe shows negligible changes in the conduction band structure, consistent with the rather broad I d states.100 For equal carrier concentrations (n = 12127

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

4c. The introduction of Cd does not extinguish the favorable resonance effect of In below 500 K but improves the Seebeck coefficients and power factors over the entire temperature range. We believe this codoping concept can improve the average ZT over a broad temperature range for many thermoelectric materials. In summary, as indicated by eq 7, the effective mass (m*) can be enhanced by increasing the valley number (NV) or by increasing local DOS effective mass (mb*) through distorting DOS. It is well known that the carrier mobility is inversely proportional to the effective mass, but the carrier mobility will be less deteriorated by increasing valley number (NV) than by increasing the local effective mass (mb*).86 3.2. Modulation Doping and Carrier Mobility Improvement

Most state-of-the-art thermoelectric materials are heavily doped semiconductors with the carrier concentrations on the order of 1019−1021 cm−3, Figure 2a. The high population of free carriers in these heavily doped semiconductors leads to a decrease of carrier mobility with respect to the lightly doped or undoped ones, because of the increased ionized impurity scattering. For example, Hilsum112 found that at room temperature, the mobility is inversely proportional to [1 + (nD/1017)1/2] for a number of important semiconductors, where nD is the donor concentration (in cm−3). Three-dimensional (3D) modulation doping has been recently demonstrated to be effective in increasing the ZT values of several important thermoelectric materials by carrier mobility enhancement.113−116 Modulation-doped samples (Figure 5b) actually are two-phase composites made of undoped (Figure 5a) and heavily doped (Figure 5c) counterparts. The undoped pristine sample has low carrier concentration but high carrier mobilities, while the uniformly and heavily doped sample has high carrier concentrations but low carrier mobilities. The Fermi level of the undoped sample usually locates at the middle of the energy gap and would move deep into the conduction band (n-type doping, Figure 5d) or the valence band (p-type doping, Figure 5e) in the case of heavily doped sample. Because of the Fermi level imbalance between the undoped and the heavily doped phases, the carriers in the modulation-doped sample spill over from the heavily doped region to the undoped region (Figure 5d and 5e). This process results in carrier mobility enhancements relative to uniform doping because of avoidance of ionized impurity scattering. The spatial separation of carriers and impurity scattering centers in 3D modulation doping enables higher carrier mobility without deteriorating the Seebeck coefficient, Figure 5b. For example, under the same carrier concentrations, the 3D modulation-doped p-type SiGe alloys113 and BiCuSeO compounds116 possess considerably larger carrier mobilities and power factors than the uniformly doped materials, Figure 5f. Recently, this 3D modulation doping approach has also been successfully applied to several other bulk thermoelectric systems (e.g., p-type ZnSb117 and n-type FeSb2118) with largely improved thermoelectric performance. One concern of modulation-doped structures is that they may functionally fail when operating at high temperature for long periods of time, due to element diffusion and thermodynamic equilibration. Therefore, modulation doping could be more promising at relatively lower temperatures. Aside from modulation doping, texturing could be another feasible approach to increase the carrier mobility, especially for

Figure 5. Schematic representation of (a) nondoped semiconductor, (b) heavily doped semiconductor where dopants are uniformly dispersed in the host matrix, and (c) modulation doping where two types of grains consisting of nondoped and heavily doped semiconductors are spatially separated. Red arrows show how carriers are scattered when mobilizing across the material. Fermi level positions in undoped, uniformly doped, and modulation-doped semiconductors: (d) n-type doping and (e) p-type doping. The Fermi level imbalance between the undoped and the heavily doped grains prompts the carriers to diffuse from the latter to the former. (f) Comparison of power factors of p-type SiGe alloys113 and p-type BiCuSeO116 in the cases of uniform doping and modulation doping, respectively. At comparable carrier concentration, materials with modulation doping clearly show much higher power factors than those by uniform doping, due to the much larger carrier mobilities.

materials with anisotropic structures. In these anisotropic compounds, carrier mobility may be higher in certain crystallographic directions. In polycrystalline sample, the random orientation of grains weakens the anisotropy of carrier transport and the carrier mobility is a statistical mean value along different crystallographic directions. To take advantage of the anisotropy and maximize carrier mobilities, it is desirable to produce grains that orient along a certain direction, so as to approach the character of the single crystals. For instance, single crystals of SnSe40 feature a high carrier mobility of ∼250 cm2 V−1 s−1 at 300 K, along the b axis of the orthorhombic unit cell, which is about 5 times larger than in its polycrystalline119 form with equal hole concentrations. The significantly higher mobility of SnSe single crystal is the main reason for its much superior thermoelectric performance as compared with its polycrystalline samples. However, as single-crystal growth is highly dependent on the equipment and is time consuming, simpler methods to realize highly oriented polycrystalline samples by hot deformation (thermoforging) technique are being developed.120 Such methods have been extensively utilized to improve the thermoelectric performance of many 12128

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

layer structured materials; for example, Bi 2 Te 3 , 121,122 Ca3Co4O9,123,124 and BiCuSeO.116 Texturing is not sensitive to temperature, so it can be applied to high-temperature thermoelectric materials. 3.3. Reducing Lattice Thermal Conductivity

In solids, the interaction between atoms creates displacements from their equilibrium positions leading to a set of vibrational waves with various wavelengths, so-called phonons. These phonons are heat carriers that propagate through the lattice, contributing to the lattice thermal conductivity κlat. At crystal defects (e.g., point defects, dislocations, interfaces, precipitates, etc.), phonon waves can be scattered, giving rise to additional thermal resistance and a reduction of κlat. For example, atomic point defects formed by elemental substitution and nanoscaled precipitates created by second-phase nucleation and growth can significantly scatter the short- and medium-wavelength phonons, respectively. Mesoscaled grains can also scatter phonons but are more effective for long-wavelength phonons.4,125 Below we will elaborate on these mechanisms. 3.3.1. Atomic Scale: Synergistic Alloying. Introduction of lattice imperfections (point defects, at the length scale of angstroms) in the host lattice (by doping or alloying) is a wellestablished approach of decreasing the lattice thermal conductivity. According to the thermal conductivity model developed by Klemens126 and Callaway,127,128 this reduction of κlat is a result of combination of mass contrast and strain field fluctuations. The degree of reduction can be evaluated by the scattering parameter (Γ) that is written as129 ⎡ ⎛ adisorder − a pure ⎞2 ⎤ ⎛ ΔM ⎞2 ⎢ ⎟⎟ ⎥ ⎟ + ε⎜ Γ = x(1 − x) ⎜ ⎜ ⎢⎝ M ⎠ a ⎝ ⎠ ⎥⎦ pure ⎣

Figure 6. Schematic representation of various types of point defects: (a) single doping, (b) cross-substitution, and (c) lattice vacancy formation. (d) Comparison of lattice thermal conductivities as a function of temperature for SnTe with different point defects.43,140,141

one Fe2+ and one Te2− concurrently, thus forming an isoelectronic compound FeSb2Te.130−132 Cross-substituted compound can have a significantly larger solubility limit of dopants than by single doping, and one could expect a much lower κlat, according to eq 10. It should be noted that crosssubstitution not only considerably impacts the phonon transport but may modify the electronic structure (band gap mostly),131,133 providing other avenues to tune the thermoelectric properties of the host material. Because only when the two aliovalent elements are present together in a specific ratio do they enter the matrix in larger amounts (than each alone can) we refer to this effect as synergistic alloying. In addition to a high doping fraction, increasing the atomic mass contrast between the guest and the host atoms can lead to low lattice thermal conductivity. However, the maximum mass contrast (ΔM/M = 1) seems unattainable by conventional doping but is possible when the guest or host atom represents a disordered vacancy, Figure 6c. For instance, in partially filled skutterudites FyCo4Sb12 (F and y represent the filler atom and filling fraction, respectively; y < 1), the lattice thermal conductivity is unexpectedly low, which has been rationally explained by the strong point defect scattering between the filler F in fully filled FCo4Sb12 and the structural vacancy □ in unfilled □Co4Sb12, a 100% mass difference.134 Other examples include clathrate compounds135 with void cages, La3−δTe4 materials136,137 with intrinsically high concentration of cationic vacancies, and InSb/GeTe/SnTe−In2Te3 solid solutions138−140 (note: one-third of the cationic sites are vacant in In2Te3) where very low lattice thermal conductivity is observed. Figure 6d compares the temperature-dependent lattice thermal conductivities of SnTe with different types of point defects: (A) single In doping,43 (B) cross-substitution by Ag and Bi,141 and (C) formation of lattice vacancies by introduction of In2Te3.140 It is clear to see that crosssubstitution and vacancy formation are superior ways of frustrating the phonon propagation in SnTe compared to single doping, as we discussed above. This synergistic alloying approach to thermal conductivity reduction is applicable to most kinds of thermoelectric materials.

(10)

where x is the doping fraction, ΔM/M is the rate of change of atomic mass, adisorder and apure represent the lattice constants of disordered and pure alloys, respectively, and ε is an elastic properties related adjusting parameter.129 From eq 10, one can conclude that to maximize Γ and achieve the lowest κlat it is necessary to have (1) a high doping fraction x, (2) a large mass difference between the dopant and the host element (ΔM/M) creating disorder in the lattice, and (3) a significant lattice mismatch between the disordered phase and the host phase (adisorder/apure). Requirements 1 and 3 can be conflicting in many cases, since large lattice and size mismatch generally renders the dopant insoluble. Figures 6a−c schematically represent the typical doping approaches frequently used in thermoelectric materials research, namely, single-element doping, cross-substitution, and formation of lattice vacancies. In the case of single doping, the dopant could be either isovalent with the host element to generate lattice disorder or aliovalent to control the carrier concentration, Figure 6a. However, due to charge imbalance, singly aliovalent doping is limited by the low doping fraction possible which prevents the further decrease of κlat. To promote a large solubility limit of aliovalent elements, cross-substitution is used. Cross substitution means replacing one or more of the host elements by pairs from other groups of the Periodic Table while keeping the overall valence electron counts constant, Figure 6b. For example, in PbTe, when the divalent Pb is simultaneously substituted by half monovalent Ag and half trivalent Sb, one could obtain AgPbmSbTe2+m, also known as LAST.1 In CoSb3, one Co3+ and one Sb1− can be replaced by 12129

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

3.2.2. Nanoscale: Nanostructuring with Second Phases. Nanoscaled defects in crystals help scatter the phonons with wavelengths on the nanometer scale. To make this scattering effective, the geometric length of the nanoscaled defects should be small enough, typically on the order of several to dozens of nanometers, and the nanostructures should be uniformly distributed for maximal interface density over a large material volume.3,4,142 There are several approaches to achieving nanoscale inhomogeneity, including ex-situ additions of guest phase by mechanical or chemical mixing143,144 and in-situ second-phase precipitation through thermodynamically/kinetically driven processes.1,3,4,145,146 While both have been reported to create nanostructures in thermoelectric materials, the latter is more favorable for charge transport because they can afford better dispersion and yield coherent interfaces between precipitate and matrix.3,4 Nucleation and growth of a second phase is the most common in-situ approach of producing nanostructures in bulk crystals,147 and its mechanism is schematically presented in Figure 7a. This works well when the second phase is completely

This nanostructuring approach has been successfully applied to numerous thermoelectric materials with a great reduction of lattice thermal conductivity. Typical examples include PbTe− MTe (M = Ag, Ca, Sr, Ba, Cd, and Hg),3,4,71,75,148,149 PbSe− MSe (M = Ca, Sr, Ba, Zn, Cd, Sb, and Bi),28,44,150 PbS−MS (M = Bi, Sb, Ca, Sr, Zn, and Cd),27,151,152 SnTe−MS (M = Cd and Zn),43,92 CoSb3−MSb (M = In, Ga, and Zn),153−156 and Bi2Te3−M (M = Cu, Zn, and Te).22,32 A specific example is shown in Figure 7c. When only ∼4% (mole fraction) of strontium chalcogenides is added as second phase to the parent lead chalcogenides, the lattice thermal conductivities are strongly decreased. Unlike point defect scattering where high alloying fraction is always pursued to achieve the lowest thermal conductivity, in the case of nanostructuring, a relatively low volume fraction of the second phase (usually a few moles %) is optimal to achieve the strongest phonon scattering. In fact, at higher fractions, the nanostructures can aggregate or coalesce into larger structures with lower interface area, weakening their scattering efficiency.130,151 Despite the numerous successes described above, we must note that nanostructuring is not always successful in scattering the heat-carrying phonons. As an example, Pb and Bi nanostructures appear to exert little influence on the lattice thermal conductivity of PbTe.157 In fact, the phonon propagation inside the nanostructures is rather complicated and usually unpredictable. One should always bear this in mind when designing materials. Moreover, at elevated temperatures where diffusion becomes more active, nanoscale precipitates may coarsen and grow, reducing their effectiveness. This challenge needs to be overcome and will require further investigations in this direction, and new insights and techniques already employed in preparing metallurgical nanostructured alloys.158 3.3.3. Mesoscale Structuring. Low- and medium-wavelength phonons carry most of the heat and can be effectively scattered by point defects and second-phase nanostructures, respectively. There are still many long-wavelength phonons remaining, however, and they also contribute to the lattice thermal conductivity. To scatter more of these phonons as well, crystal defects and structure features at the length scale of micrometers or submicrometers is necessary. This is the socalled mesoscale, and it can be engineered into materials through careful powder processing. Typically, the single crystals or cast ingots (Figure 8a) have grains that are rather large (typically > 10 μm in geometric size). A mesoscale sample typically has grain sizes from 0.1 to 3 micrometers, Figure 8b. Rowe et al.159 investigated the thermal conductivities of SiGe alloys with different grain sizes and compared them with that of their single crystals, Figure 8c. They showed that large grain sizes (light blue region where cast ingots locate) are ineffective in decreasing the thermal conductivities. Only grains with relatively small sizes (5 μm and below, light red area where the polycrystalline samples locate) could scatter phonons strongly (∼20% reduction of thermal conductivity with respect to the single crystal). This work highlights the importance of mesoscale grains in impeding the phonon transport of bulk materials and points to further research to better understand mesoscale effects. To fabricate polycrystalline bulk thermoelectric materials, one usually needs to first produce mesoscopic particles and then densify them by either hot pressing or spark plasma sintering. There are two general routes of obtaining the mesoscopic particles, namely, top-down and bottom-up ap-

Figure 7. (a) Schematic ternary diagram, which indicates a strong temperature-dependent solubility of the second-phase B within the matrix A. (b) Typical TEM image that shows B nanostructures embedding in the A matrix, as grown by the nucleation and growth mechanism. (c) Comparison of lattice thermal conductivities between pristine and nanostructured lead chalcogenides.4,27,150 Strong reduction of lattice thermal conductivity is achieved through nanostructuring.

soluble in the parent material when molten (region I) but has very low solubility in the solid state (region II). When the molten mixture (region I) is rapidly cooled to form a solid, the second phase precipitates (region III) since the solid solution limit is exceeded. Moreover, to ensure a fine distribution and good thermal stability, the second phase should have a comparable or higher melting point than the host material. Figure 7b displays typical nanostructuring created by the nucleation and growth mechanism under the observation of transmission electron microscopy (TEM). The sample is PbTe−SrTe composite and features ubiquitous nanocrystals of second-phase SrTe coherently dispersed in the host material of PbTe. This good interface originates from the homogeneous dispersion of SrTe in PbTe enabled by fast cooling as well as their very small lattice mismatch. 12130

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

Figure 8. Schematic diagrams of microstructures of (a) single crystals or cast ingots without grain boundaries or with rather coarse grains (>10 μm in size) and (b) polycrystalline samples with mesoscale grains (hundreds of nanometers to several microns). (c) Experimental study of the thermal conductivity of SiGe alloys (green circles) with different grain sizes as compared to the single crystal.159 Solid line is a guide to the eyes. Light blue and red areas denote where the grain sizes of cast ingot and polycrystalline samples locate, respectively.

Figure 9. Schematic diagram of the integration of multisized defects: (a) atomic-scale point defects, (b) nanoscale precipitates, and (c) mesoscale grains in one single material for all-scale hierarchical architectures of phonon scattering.4 (d) Contributions of phonons with different mean free paths to the cumulative κlat value for PbTe.125 (e) Temperature dependence of lattice thermal conductivity for the all-scale hierarchical architecture SnTe system.88,92 Green color depicts the mesostructuring, red color depicts the approach of solid solution point defects, and blue line depicts the approach of nanostructuring.

proaches.160 In the top-down approach, the mesoscopic particles are acquired by pulverizing the coarse-grained crystals to fine powders with the aid of hand grinding, mechanical alloying, etc. Alternately, as developed by Tang’s group, using one-pot melt spinning or self-propagation high-temperature synthesis (SHS), one could obtain the products with fine grains or small particle size. This has proven to be particularly effective in producing high-performance bulk thermoelectric materials in a very short period.161−164 Inversely, the bottom-up approach means direct control of the mesoscopic particles growth by, such as, solution chemistry methods including hydro- and solvothermal synthesis,165 electrochemical deposition,166 sonochemical synthesis,167 and so on. These solution chemistry methods have been widely used to synthesize a wide variety of thermoelectric semiconductors (e.g., skutterudites,168 metal chalcogenides,167 zinc antimonides169) with controllable size and morphology. 3.3.4. All-Scale Hierarchical Architectures. When point defects, nanostructuring, and mesoscale structuring are all combined into a single thermoelectric material, we call this as an all-scale hierarchical architecture.4 While the thermal conductivity in crystalline materials is limited by the Umklapp scattering process, it can be decreased by defects at different length scales with a broad range of wavelengths and mean free paths (MFP, defined as the average distance traveled by phonons between successive scatterings). Theoretical calculations125 suggest that, in PbTe, over 50% of the κlat value comes from the contribution by phonons with MFP less than 1 nm, which can be ascribed to scattering by atomic-scale point defects, Figure 9c. The remaining contribution to κlat is almost equally divided between phonon modes with MFP of 1−10 nm and phonon modes with MFP 10−1000 nm, Figure 9c), which can be notably impeded by nanoscale precipitates and mesoscale grains, respectively. Similar trends are found in PbTe1−xSex170 and Si.171 All-scale hierarchical architectures achieve the strongest phonon scattering,83,148,155,172 see Figure 9a. All-scale hierarchical architectures were successfully created in lead chalcogenides, have led to significantly enhanced thermoelectric performance, and now have been applied to other lead-free

materials, including SnTe, 43,92 GeTe, 38 AgSbSe 2 , 173 AgCrSe2,174 CoSb3,155,164 Bi2Te3,175 and In2O3.176 Figure 9e shows a comparison of the lattice thermal conductivity of SnTe samples with different defects: mesoscale grains by SPS treatment,88,92 atomic-scale point defects of Cd substitution for Sn,92 and nanoscale CdS precipitates.92 It is apparent that mesostructuring lowers the lattice thermal conductivity of SnTe ingot by about 22% at room temperature, from ∼2.7 to ∼2.1 W m−1 K−1. The atomic point defects and nanoscale precipitates contribute additional phonon scattering, leading to further reduction of κlat by ∼25% and ∼20%, respectively. In all, a total reduction of κlat by more than 50% is achieved by all-scale hierarchical architectures of phonon scattering.

4. DECOUPLING OF ELECTRON AND PHONON TRANSPORT We highlighted above that nanostructuring in bulk materials dramatically reduces the thermal conductivity. However, simultaneously it can also increase the charge carrier scattering due to the additional energy barrier caused by crystallographic mismatch and/or electronic bands misalignment at the interfaces.4,43,83 This has detrimental effects on the carrier mobility and power factor. Therefore, to avoid the mobility loss in nanostructured thermoelectric materials, both coherent interfaces (endotaxy) and energy-matched electronic bands are required. To some extent this situation achieves a decoupling of electron and phonon transport. Below we discuss approaches of decoupling electron and phonon transport in nanostructured bulk thermoelectric materials. 4.1. Strained Endotaxial Nanostructuring

In a high-performing nanostructured thermoelectric material, the phonons should be scattered as much as possible by the nanoprecipitates while the electrons should flow across the interfaces without much scattering, Figure 10a. The electron and phonon scattering are strongly affected by the interface (endotaxy), which is usually determined by the degree of lattice misfit (ε) between the nanoprecipitates and the matrix.157 12131

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

demonstrated, including PbTe−GeTe,38 PbTe−PbS,42 PbTe− MTe (M = Mg, Ca, Sr, Ba),3,4,83,148 PbSe−MSe (M = Ca, Sr, Ba),150 and SnTe−MnTe.89,97 Figure 10 c−f show strained endotaxial nanostructuring present in the PbTe−SrTe system. Figure 10c is a typical phase contrast high-resolution TEM image showing several 2−4 nm SrTe precipitates (depicted by the white arrows, marked as i, ii, iii, and iv, respectively) within the PbTe matrix. An enlarged high-resolution TEM image of precipitate (i) shown in Figure 10d clearly suggests the absence of any dislocations at the boundary highlighted by the dotted line. In Figure 10e, the inverse first Fourier transform reconstructed image of the precipitate (i) depicts a perfect interface endotaxy (coherency) between the precipitate and the matrix. The shear strain map profiles based on geometric phase analysis conducted on the four nanoprecipitates shown in Figure 10c are displayed in Figure 10f.3 Elastic strain is pervasive in and around all precipitates, and there is also plastic strain around the dislocation cores in precipitates iii and iv. The coherent but strained interfaces between the matrix and the nanoprecipitates lead to strong phonon scattering but minimal electron scattering. Strained endotaxial nanostructuring is now becoming a popular approach of designing high-performance bulk thermoelectric materials. 4.2. Matrix/Precipitate Band Alignment

Figure 10. Schematic diagram of electron and phonon transport decoupling in nanostructured bulk materials: phonons (depicted by blue arrows) are strongly scattered by nanoprecipitates, while the electrons (depicted by red arrows) travel across the interfaces without restriction. (b) Schematic of phonons scattered by three possible types of interfaces (coherent, semicoherent, and incoherent) between matrix and nanoprecipitates. Reprinted with permission from ref 157. Copyright 2010 American Chemical Society. (c) High-resolution TEM phase contrast image of several SrTe nanoprecipitates (i, ii, iii, and iv) embedded in the PbTe matrix. (d) Enlarged image of nanoprecipitate i showing coherency at the boundary highlighted by the dotted line. (e) Inverse first Fourier transform image of precipitate i, indicating the absence of dislocations at the grain boundary. (f) Shear strain map profiles for the four nanoprecipitates shown in d. There are elastic strains at and around all precipitates and also plastic strain at and around dislocation cores in precipitates iii and iv. (c−f) Reprinted with permission from ref 3. Copyright 2011 Nature Publishing Group.

Though the scattering of electrons can be minimized by improving the interface coherency in nanostructured materials, the additional energy barrier from misalignment of electronic bands can hinder electron flow. Previous studies on GaAs/InAs (001) heterojunctions indicated that charge carriers can transmit freely when the band offsets are minimal.182 Recently, a similar approach has been developed in three-dimensional bulk nanostructured thermoelectric materials by minimizing the energy offsets of electronic bands between the nanoprecipitates and the matrix.3,4,27,28 The concept of electronic band alignment is illustrated in Figure 11a. The dotted lines depict the energies of the conduction band minimum (CBM) and the valence band maximum (VBM) of two different second phases and a host matrix. Phase A has good valence band alignment with the matrix (an energy difference that is comparable to the magnitude of thermal energy at elevated temperatures of kBT is considered band alignment). Therefore, the matrix/A system

There are generally three types of precipitate−matrix interfaces in light of their coherency: coherent (ε < 1%), semicoherent (1% < ε < 25%), and incoherent (ε > 50%).157,177 These three scenarios embody different atomic configurations and associated local interfacial relaxation mechanisms,157 Figure 10b. In the case of semicoherent or incoherent interfaces, the surrounding atomic order is strongly disrupted, leading to significant scattering of both electrons and phonons. The direct consequence is limited ZT improvement. When the nanoprecipitates and the matrix are coherently (endotaxially) connected but have a certain degree of strain around their interfaces, we call it strained endotaxial nanostructuring. This interface is supposed to facilitate charge transport across it but strongly block/scatter phonon propagation. The presence of strained endotaxial interfaces has been observed in a number of nanostructured materials. In pioneering work by Hsu et al.,1 coherent nanostructured bulk thermoelectric materials LAST were synthesized with a measured ZT of ∼1.7 at 700 K. They contained Ag−Sb-rich nanoprecipitates embedded within the PbTe matrix. Subsequently, the nature of the LAST materials was analyzed in detail by several groups.178−181 Later, other examples were

Figure 11. (a) Schematic representation of the alignment of the valence band (VB) and conduction band (CB) energies of A and B second phase in the matrix. ΔEC and ΔEV denote the energy offsets of CB and VB between the second phase and the matrix, respectively. The small band offsets allow seamless carrier transmission. (b) Temperature dependence of carrier mobilities and lattice thermal conductivity for the PbTe−SrTe system doped with 2% Na.3,4 All samples have similar doping level. The endotaxially nanostructured SrTe with small valence band offsets has negligible effect on the hole mobility while blocking heat-carrying phonons at the matrix/ nanostructure interface. 12132

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

is expected to enable seamless hole transmission across the interface. Likewise, phase B has negligible energy difference of CBM with that of the matrix, and it is expected to allow facile electron transmission. p-Type PbTe embedded with well-dispersed SrTe nanocrystals is a well-known example of valence band alignment.3,4 As shown in Figure 11b, the lattice thermal conductivity of PbTe is remarkably decreased upon the addition of SrTe second phase but the hole mobility is largely unaffected. Specifically, under similar hole doping, the room-temperature lattice thermal conductivity of the control sample (containing no SrTe) is ∼1.8 W m−1 K−1, and it is decreased to ∼1.5 and ∼1.3 W m−1 K−1 by introducing 1% and 2% SrTe, respectively. In contrast, the carrier mobility is almost constant in this process, stabilizing between 340 and 350 cm2 V−1 s−1 at room temperature. The reason for this is 2-fold:3 (1) endotaxial alignment of SrTe and PbTe phases with coherent interfaces, as we elaborated in the previous section, and (2) the alignment of the energies of the two valence bands in the two materials where the band offset is very small at T > 300 K. Actually, the generality of the matrix/precipitate band alignment concept is also evident in the PbS system nanostructured with various metal sulfides (CdS, ZnS, CaS, and SrS).27 Though all these nanostructures are equally effective in achieving a very low lattice thermal conductivity (∼0.6 W m−1 K−1 at 923 K), only CdS helps maintain a high hole mobility due to its smallest valence energy offsets with PbS (∼0.13 eV at 0 K). This is the main reason for the superior ZT of the CdS-containing samples (∼1.3 at 923 K in 2.5% Nadoped PbS with 3% CdS), a record value for PbS.27

Figure 12. (a) Schematic representation of how compositionally alloyed nanostructures can achieve aligned electronic bands. (b) Valence band alignment by compositionally alloyed nanostructuring in p-type PbSe system.28 (c) Room-temperature carrier mobility of 2% Na-doped PbSe nanostructured with CdSe, ZnSe, CdS, and ZnS as a function of second-phase mole fraction.28 (d) Temperature-dependent lattice thermal conductivity of 2% Na-doped PbSe nanostructured with 3% CdSe, ZnSe, CdS, and ZnS.28

These happen via an interfacial reaction between the selenide matrix and the sulfide second phase and are confirmed by the chemical composition analysis using high-resolution energydispersive spectroscopy (EDS).28 The compositionally alloyed new phases CdS1−xSex and ZnS1−xSex then have different band energies and can be employed to create minimal energy offsets with VBM of PbSe, Figure 12b. Because of the good band alignment the compositionally alloyed nanostructures have almost no detrimental effect on the carrier mobility, Figure 12c. In contrast, a decreasing trend of mobility with increasing second-phase fraction is observed for the CdSe- or ZnSe-containing PbSe samples due to incorrect band alignment with PbSe, Figure 12b. Of course, the embedded nanostructures continue to scatter phonons to produce large reductions in lattice thermal conductivity, Figure 12d. Therefore, the compositionally alloyed nanostructuring samples feature superior performance, e.g., Pb0.98Na0.02Se−3% CdS has ZT of ∼1.6 at 923 K, the highest value ever reported for PbSe.28

4.3. Compositionally Alloyed Nanostructures

The ideal nanostructuring case is when the electronic bands of matrix and the second phase are perfectly aligned in energy. In practice it is not easy to find such an appropriate combination of second phases. The CBM or VBM of the possible candidates lie either higher or lower in energy with respect to those of the matrix. In this case we can use the concept of “compositionally alloyed nanostructures”28 to achieve electronic band alignment. The concept is schematically presented in Figure 12a. The CBM of the second phases AB and AC is too high or too low with respect to the matrix to achieve effective conduction band alignment on their own. However, if one could mix them in a proper ratio to form a solid solution (AB1−xCx) then this alloy may have minimal energy offsets with the matrix. Similarly, valence band alignment can be achieved by alloying a proper amount of second-phase XZ in XY as nanostructures (XY1−mZm). In this regard, compositionally alloyed nanostructure is a universal concept of matrix/precipitate band alignment. A good example of applying this idea is p-type PbSe compositionally alloyed with CdS or ZnS.28 As shown in Figure 12b, based on the theoretical band calculation results, the VBMs of CdS (∼0.27 eV) and ZnS (∼0.30 eV) are far below that of PbSe, while those of CdSe (∼0.06 eV) and ZnSe (∼0.13 eV) are higher than that of PbSe. Therefore, none of them could be directly used to align the valence bands. However, when CdS or ZnS is added to PbSe, the following alloying reactions are expected to occur PbSe + CdS → PbSe1 − xSx + CdS1 − x Sex

(11)

PbSe + ZnS → PbSe1 − xSx + ZnS1 − x Sex

(12)

5. DISCOVERIES OF NEW THERMOELECTRIC MATERIALS WITH INTRINSICALLY LOW THERMAL CONDUCTIVITY From the above it is apparent that a number of useful concepts have emerged, which can be integrated in a synergistic manner in bulk materials to greatly enhance thermoelectric performance.183 Clearly, we are in a new era for thermoelectric materials research. In a very special group of materials, high thermoelectric performance can be found in pristine compounds with intrinsically low thermal conductivity. The intrinsically low thermal conductivity means that the most complex approaches to reduce thermal conductivity could be avoided. Below we introduce some promising thermoelectric materials with intrinsically low thermal conductivity, which may arise from an anharmonic and anisotropic bonding,40,56 lattice vibra12133

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

tions,184,185 atomic disorders induced by ion migrations,185 copper ion liquid-like behaviors, 36,186 large molecular weight,55,187 complex crystal structure,188 lone-pair electrons,189,190 etc.

moderate ZT value of approximately 0.8 was measured along the a axis, Figure 13b. These high ZT values observed in SnSe single crystals are not observed in polycrystalline samples of SnSe, which exhibit a ZT peak of ∼1.2 at 973 K. The reason is the big differences in carrier mobility and thermal conductivity between SnSe single- and polycrystalline samples where the three different crystallographic directions are averaged but in a difficult to control manner because of the varying degrees of preferential orientation, Figure 14a.

5.1. Layered SnSe

Recently, a record high ZT of 2.6 in a single binary crystal of SnSe was observed.40 SnSe has an unusual layered structure that derives from a three-dimensional distortion of the NaCl structure. SnSe has two-atom-thick SnSe slabs (along the b-c plane) with strong Sn−Se bonding within the plane of the slabs, which are then linked with weaker Sn−Se bonding along the a direction, see Figure 13a. The two-atom-thick SnSe slabs

Figure 13. (a) SnSe crystal structure viewed along the b axis: gray, Sn atoms; red, Se atoms. Blue box represents a unit cell. (b) ZT values as a function of temperature for SnSe single crystals along different axial directions and polycrystalline SnSe pellets along radial and axial directions.40

Figure 14. (a) Carrier mobility and (b) power factor as a function of temperature for undoped SnSe single crystals40 along different axial directions and polycrystalline SnSe pellets119 along radial and axial directions. (c) Thermal conductivities as a function of temperature of undoped SnSe single crystals40 and polycrystalline SnSe prepared in air and a glovebox.192

are corrugated, creating a zigzag accordion-like projection along the b axis. The easy cleavage in this system is along the (l00) planes.40 The special crystal structure results in impressively low thermal conductivities of 0.7 (300 K) and 0.4 W m−1 K−1 (700 K) for SnSe crystals (b axis). To probe the origin of the intrinsically low thermal conductivity of SnSe, the phonon and Grüneisen dispersions, which reflect an anharmonicity of chemical bonds, have been calculated using first-principles density-functional theory (DFT). The results show that the average Grüneisen parameters along the a, b, and c axes are 4.1, 2.1, and 2.3, respectively. Along the a axis, the maximum longitudinal acoustic Grüneisen parameter is extraordinarily high at ∼7.2. In contrast, the Grüneisen parameters for known low thermal conductivity materials is 2.05 for AgSbTe2,191 3.5 for AgSbSe2,189 and 1.45 for PbTe.191 These materials have lattice thermal conductivities at room temperature of 0.68, 0.48, and 2.4 W m−1 K−1, respectively. The anomalously high Grüneisen parameter of SnSe is a reflection of its crystal structure which contains much distorted SnSe7 polyhedra (due to the lone pair of Sn2+) and the zigzag accordion-like geometry of slabs in the b-c plane. This implies a very soft lattice, in which if mechanically stressed along the b and c directions would not directly change the Sn−Se bond lengths, but instead the zigzag and geometry would allow it to be deformed like a retractable spring. In addition, along the a direction the weaker bonding between SnSe slabs makes it a good stress buffer or a “cushion”, thus dissipating phonon transport laterally. The anomalously high Grüneisen parameter is therefore a consequence of “soft” bonding in SnSe, which presumably leads to the very low thermal conductivity and the surprising thermoelectric performance of the single crystals.40 Single crystals of SnSe show a ZT value of approximately 2.6 at 923 K along the b axis and 2.3 along the c axis, whereas a

The carrier mobility along the b axis for SnSe single crystals is 10 times higher than along the a axis. This strong anisotropic character diminishes the carrier mobility in polycrystalline SnSe, as also observed by Chen et al.119 For example, at room temperature, the hole mobility of polycrystalline SnSe is ∼50 cm2 V−1 s−1 with a carrier density of 3 × 1017 cm−3;119 however, the hole mobility in single crystals is about 250 cm2 V−1 s−1 (at similar carrier density of 3 × 1017 cm−3). At 750 K, at the carrier density of 4 × 1018 cm−3, the hole mobility of polycrystalline SnSe sample is about 10 cm2 V−1 s−1, which is only one-fifth of ∼50 cm2 V−1 s−1 for the SnSe single crystals.40 The low carrier mobility of polycrystalline SnSe results in an electrical conductivity that is five times lower than that of SnSe single crystals. This results in a one-fifth of the power factor (Figure 14b) and lower ZT values (Figure 13b). We expect this number to change from report to report and from sample to sample since the degree of polycrystalline orientation may vary across samples and laboratories. For this reason we suggest such reports of transport properties of polycrystalline SnSe to be accompanied with detailed X-ray and microscopic characterization of the degree of preferential orientation. Very recently, Zhang et al.192 reported a similar thermal conductivity in polycrystalline SnSe to that of SnSe single crystals.40 The maximum power factor in polycrystalline SnSe achieved to date is ∼4 μW cm−1 K−2 at 773 K and a ZT of 1.0.192 Correspondingly, the maximum power factor is ∼6 μW cm−1 K−2 at 773 K, resulting in a ZT of 1.5 at 773 K in undoped SnSe single crystals.40 It should be noted that both these reports have samples of similar carrier concentration and Seebeck coefficient. Therefore, the differences in power factor 12134

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

and ZTs come from the carrier mobility and from the strong anisotropic character of the SnSe crystal structure. Apart from the carrier mobility in some reports,90,119,193−196 the ZT differences between SnSe single crystals and polycrystalline SnSe samples also come from thermal conductivity deviations. Indeed, the explanations of strong anharmonic chemical bonding given for the low thermal conductivity of SnSe single crystals also apply to polycrystalline SnSe samples.40,107 However, accurately measuring the thermal conductivity of SnSe even in single-crystal samples seems to be a challenge. It has been found that polycrystalline samples are more difficult to work with because of surface oxidation and preferential orientation. It is possible that the surface oxidation yields SnO2 second phase, which has a thermal conductivity of ∼98 W m−1 K−1, about 140 times higher than SnSe.197 If the polycrystalline samples are prepared or handled in air this could give higher thermal conductivity values. In our experience polycrystalline samples protected from oxidation have low thermal conductivities similar to single-crystal samples.192 Recently, He et al. gave the same explanation for the thermal conductivity issue.198 They found that the thermal conductivity of polycrystalline SnSe will decrease when the SnSe powders are processed in a glovebox and will closely approach the 0.30 W m−1 K−1 value at 773 K, which is comparable to 0.30 (c axis) and 0.35 W m−1 K−1 (b axis) of SnSe single crystals.40 This oxidation could lead to a thin layer of tin oxides on the surface of the micrometer-sized powders. Indeed, a similar low thermal conductivity of 0.30 W m−1 K−1 at 773 K was observed in polycrystalline SnSe strictly prepared in a glovebox,192 as shown in Figure 14c. To investigate this hypothesis, a series of control experiments was carried out. We find that the thermal conductivity of polycrystalline SnSe is very low when the SnSe powders are processed in a glovebox compared to processing in air and closely approaches the 0.36 W m−1 K−1 value at 773 K, which is comparable to 0.30−0.35 W m−1 K−1 of SnSe single crystals.40 Interestingly, when we processed the SnSe powders with a small added amount of SnO2 (e.g., the SnO2 addition ranged from 0.0038 to 0.02 g of SnO2 per 5 g of SnSe, 0.1−0.6% by weight), the thermal conductivities of the polycrystalline SnSe samples showed a significant increase even though the powders were processed in a glovebox. Besides the preparation conditions, stoichiometric composition is also crucial to the thermal conductivity of SnSe. Serrano-Sanchez et al. confirmed the intrinsically low thermal conductivity of SnSe and provided insights into the stoichiometric ratio and the sample processing conditions.199 The microscopic neutron powder diffraction analysis demonstrated a nearly ideal stoichiometry and a high amount of anharmonicity of the chemical bonds of SnSe, and the lowest thermal conductivity value is close to 0.1 W m−1 K−1 at room temperature for the polycrystalline SnSe.199 As an analogue of SnSe, the thermoelectric properties of SnS are also receiving attention. Experimentally, Tan et al. reported that the low thermal conductivity falls below 0.5 W m−1 K−1 at 873 K and leads to a high ZT of 0.6 in Ag-doped polycrystalline SnS, pointing out that the low-cost SnS is a promising candidate for thermoelectric investigations.200,201 The experimental results are well supported by first-principles theoretical calculations, which indicate that SnS could possess a high Seebeck coefficient and a very low thermal conductivity.202 Bera et al. suggested that SnS is potentially a good thermoelectric material if it can be suitably doped.203 Considering the poor electrical transport properties of SnS, a

high thermoelectric performance could be expected in perfect SnS single crystals with heavy doping. Apart from the intrinsically low thermal conductivity of SnSe, the electronic band structure is also very impressive. It is well known that the Seebeck coefficients could be enhanced through manipulating multiple valence or conduction bands, which are highly sought features in thermoelectric materials. Figure 15a

Figure 15. (a) Schematic of the electronic band structure of SnSe indicating nonparabolic, complex multiband valence states. The Fermi level (green dash line) is moving into the multiband structure after heavily doping. Electrical transport properties as a function of temperature for undoped SnSe and hole-doped SnSe along different crystallographic directions: (b) electrical conductivity; (c) Seebeck coefficients (inset shows the room-temperature Seebeck coefficients comparisons for the lead and tin chalcogenides with similar carrier concentration of ∼4 × 1019 cm−3); (d) power factors including those of the optimized Bi2−xSbxTe3 and PbTe−4SrTe−2Na plotted for comparison.2,4

shows multiple valence bands for SnSe, and it is rare to find a system that exhibits this enhancement while also exhibiting extremely low thermal conductivity and a high electrical conductivity.204,205 Hole doping increases the electrical conductivity from ∼12 to ∼1500 S cm−1 as the carrier concentration increases from ∼1017 to ∼1019 cm−3 at 300 K, as shown in Figure 15b. Meanwhile, the Seebeck coefficient decreases from ∼+500 μV K−1 for undoped SnSe to ∼+160 μV K−1 at 300 K in doped while still keeping a large value of ∼+300 μVK−1 at 773 K. The inset of Figure 15c indicates that the Seebeck coefficient for hole-doped SnSe is clearly much higher than ∼+70 μV K−1 for PbTe,206 ∼+60 μV K−1 for PbSe,207 ∼+50 μV K−1 for PbS,27 and ∼+25 μV K−1 for SnTe92 at similar doping levels. The observed Seebeck coefficient enhancements of hole-doped SnSe can be attributed to the multiband character of the electronic structure.204,205 The combination of vastly increased electrical conductivity and high Seebeck coefficient results in a large power factor of ∼40 μW cm−1 K−2 for hole-doped SnSe (b axis) at 300 K, Figure 15d. The power factors obtained in hole-doped SnSe are as high as those of optimized Bi2−xSbxTe3 materials2 near room temperature and are higher than those of the high-performance hierarchically architectured p-type PbTe−SrTe system in the range of 300−500 K.4 The temperature dependence of total thermal conductivities (κtot) for undoped and hole-doped SnSe are shown in Figure 16a. At room temperature, the values of κtot for undoped SnSe are ∼0.46, 0.70, and 0.68 W m−1 K−1 along the a, b, and c axis 12135

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

Figure 17. BiCuSeO crystal structure and ZT values. (a) Crystal structure along the a axis: yellow, Cu atoms; red, Se atoms; purple, Bi atoms; blue, O atoms. Blue box represents a unit cell. (b) Dimensionless figure-of-merit ZT vs temperature of selected typical current thermoelectric materials: CsBi 4 Te 6 , 29 Bi−Sb−Te, 2 AgPbmSbTem+2 (LAST),1 Mg2SiSn,24 PbTe−SrTe−Na,4 skutterudite,33 BiCuSeO,116,213 half-Heusler,214 SiGe,6 and Zintl phase.55,187 Green lines are p-type, blue lines are n-type materials, and thick red line is BiCuSeO.

Figure 16. (a) Total thermal conductivities, (b) lattice thermal conductivities and (c) ZT values as a function of temperature for undoped SnSe and hole-doped SnSe along different crystallographic directions; (d) ZT values comparisons of hole doped SnSe (b axis) and the current state-of-the-art p-type thermoelectrics, BiSbTe, 2 MgAgSb,34 NaPbmSbTem+2 (SALT),208 PbTe-4SrTe-2Na4 and SnSe (b axis).40

BiCuSeO exceeds 1.0 above 650 K and monotonously increases up to 1.4 at 923 K.116,212,213 The average ZT, quantitatively represented by the area under the ZT vs temperature curve, is the property used to estimate the performance in real application conditions. It can be seen that the average ZT of BiCuSeO outperforms several conventional thermoelectric materials. Moreover, compared to conventional heavy-metalbased thermoelectric materials,5,50,56,96 BiCuSeO oxyselenides are earth abundant and can thus be regarded as promising candidates for vehicle waste heat recovery and solar thermal engine replacement. The promising ZT of the oxyseleneide compound BiCuSeO is mainly due to its intrinsically low thermal conductivity.215−218 Ding et al. calculated the thermal conductivity of BiCuSeO using the lattice thermal conductivity obtained by the Grüneisen parameter and the electronic thermal conductivity using the Wiedemann−Franz relation.185 As shown in Figure 18a, at low temperature the experimental data are far below the calculated results, even after the addition of scattering contributions from the defects. This indicates an extra mechanism may be responsible for the low thermal conductivity of BiCuSeO. Figure 18b shows the heat capacity of BiCuSeO as a function of temperature; the heat capacity at room temperature shows good agreement with the Dulong− Petit law. Saha184 proposed an atomic displacement pattern for the respective lowest frequency optical mode of BiCuSeO, as shown in Figure 18c. In this pattern the arrows are proportional to the amplitude of the atomic motions, since the Bi atoms exhibit a significantly large displacement, indicating higher anharmonic effects.184 Assuming the heat is conducted only by acoustic and quasiacoustic phonons via umklapp scattering processes, the calculations suggest an average Grüneisen parameter of 2.9. Additionally, the calculated thermal conductivity in the out-of-plane direction of BiCuSeO is about two times smaller than that in-plane one, indicating that BiCuSeO has an large anisotropy of heat flow between the inplane and the out-of-plane directions.184 Ding et al. also calculated the phonon spectrum and mode Grüneisen parameters of BiCuSeO.185 First-principles calculations provide direct evidence of the structural in-layer and interlayer off-phase vibration modes, which contribute to the anharmonic vibrations and structural scattering of phonons and

directions, respectively. These already low values continue to decrease with rising temperature, and all fall in the range 0.25− 0.28 W m−1 K−1 at 773 K. After hole doping, the total thermal conductivities are higher due to a large contribution from electronic thermal conductivities. To derive the lattice thermal conductivity, the Lorenz number of 1.5 × 10−8 V2 K−2 was used for the undoped SnSe since it is a nondegenerate semiconductor. To more accurately obtain an estimate of the lattice thermal conductivity of hole-doped SnSe, the Lorenz number L has to be calculated based on a multiband model. Using more appropriate Lorenz numbers, one can see that the lattice thermal conductivity of hole-doped SnSe is comparable to or even lower than undoped SnSe, Figure 16 b. After hole doping, one can find a vast increase in ZT from 0.1-undoped to 0.7hole-doped along the b axis at 300 K while obtaining the ZTmax of 2.0 at 773 K, Figure 16 c. The hole-doped SnSe (b axis) outperforms most of current state-of-the-art p-type materials at 300−773 K, such as Bi2−xSbxTe3,2 MgAgSb,34 NaPbmSbTem+2 (SALT),208 PbTe−4SrTe−2Na4 and SnSe (b axis),40 as shown in Figure 16 d. 5.2. BiCuSeO Oxyselenides

Since BiCuSeO oxyselenides have been reported as potential thermoelectric materials in 2010 by Zhao et al.,209 they have attracted everincreasing attention and have been extensively studied. BiCuSeO crystallizes in a layered ZrCuSiAs structure type, with the tetragonal unit cell a = b = 3.9273 Å, c = 8.9293 Å, Z = 2, and the space group P4/nmm (No.129, PDF 450296).56 Figure 17a shows that BiCuSeO exhibits a twodimensional structure containing fluorite-like Bi2O2 layers alternatively stacked with Cu2Se2 layers along the c axis. The coordination of Bi can also be considered as a distorted square antiprism with four O atoms in one base and four Se atoms in the other. The Bi−Se bond length ∼3.2 Å is longer than the Bi−O bond length ∼2.33 Å, which implies its layered feature. The Cu2Se2 layers consist of slightly distorted CuSe4 tetrahedra with shared Se−Se edges. The Cu2Se2 layer could be considered as a reversed version of Bi2O2 layer, because Cu atoms occupy the Bi sites while Se atoms occupy the O sites.210,211 As shown in Figure 17b, the maximum ZT of 12136

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

Figure 19. (a) α-MgAgSb crystal structure: light blue, Ag atoms; red, Sb atoms; light purple, Mg atoms. (b) ZT values as a function of temperature for doped samples of MgAgSb.34,223

axis, with one-half of the Mg−Sb pseudocubes filled with Ag. The unit cell of the α-MgAgSb phase is doubled (distorted Mg−Sb rocksalt lattice) along both axes to accommodate a different ordering of the Ag atoms. In the pattern of filling of the Mg−Sb cubes (distorted Mg−Sb rocksalt lattice), the Agfilled sites form one-dimensional chains, running in all three primary directions. The presence of complicated phase transitions makes it difficult to obtain pure α-MgAgSb phase when the samples are cooled from 1250 K to room temperature and the impurities degrade the thermoelectric performance. For example, a αMgAgSb sample with high content of secondary phase of Ag3Sb (∼19%) shows a maximum ZT of only ∼0.5 at 425 K.221 To eliminate the impurity phase, Ying et al. carried out a 2-week annealing treatment at 553 K,222 which led to a maximum ZT of ∼0.8 at 525 K. After optimizing the carrier density using In doping on the Sb site, a maximum ZT value of ∼1.1 was achieved at 525 K at an optimized carrier density of 8−9 × 1019 cm−3, Figure 19b. Compared with conventional melting techniques, mechanical alloying has been widely used in the thermoelectric field to nanostructuring grain sizes and purifying material phases.2 Zhao et al. recently used a two-step mechanical alloying method and successfully obtained α-MgAgSb single phase, which shows a ZT of ∼0.9 at 450 K.34 Figure 19b shows that the thermoelectric performance of α-MgAgSb can be significantly enhanced through doping, namely, the maximum ZT values can be achieved ∼1.2, ∼1.4, ∼1.3, and ∼1.3 for MgAg0.97Sb0.99,34 MgAg 0.965 Ni 0.005 Sb 0.99 , 34 MgAg 0.963 Cu 0.007 Sb 0.99 , 223 and Mg0.9925Na0.0075Ag0.97Sb0.99,223 respectively. Experimentally, a high thermoelectric conversion efficiency of 8.5% was reported with a single leg based on α-MgAgSb-based compound operating between 300 and 520 K. Within the temperature range of 300−450 K, the efficiency of a MgAgSb-based device can match that of a Bi2Te3-based device.224 Considering the abundantly available constituent elements, α-MgAgSb is promising for power generation at 300−550 K. The promising thermoelectric properties of MgAgSb are ascribed to the low thermal conductivity. As shown in Figure 20a, the thermal conductivities show a decreasing trend up to 400 K and then increase with rising temperature because of bipolar diffusion. The minimum thermal conductivities at ∼373 K are ∼0.98 and ∼0.95 W m−1 K−1 for MgAgSb and MgAg0.97Sb0.99In0.01, respectively, and can be further reduced by point defects scattering through doping,222 namely, the minimum thermal conductivities at ∼373 K are ∼0.74, ∼0.66, ∼0.71, and ∼0.77 W m −1 K −1 for MgAg 0.97 Sb 0.99 , 34 MgAg 0.965 Ni 0.005 Sb 0.99 , 34 MgAg 0.963 Cu 0.007 Sb 0.99 , 223 and Mg0.9925Na0.0075Ag0.97Sb0.99,223 respectively.

Figure 18. (a) Thermal conductivity as a function of temperature for pristine BiCuSeO, the blue line is the calculated total thermal conductivity, and the black line the experimental data.185,215 (b) Heat capacity as a function of temperature for BiCuSeO; the red line is the data based on the Dulong−Petit law.215,216 (c) Atomic displacement patterns for BiCuSeO; arrows are proportional to the amplitude of the atomic motions.184 (d) Mode Grüneisen parameters of BiCuSeO. Red squares correspond to the mode Grüneisen parameters of acoustic branches, black circles correspond to the frequencies between 2.5 and 6 THz, and the blue diamond represents the phonon modes of O above 6 THz.185

result in a low thermal conductivity for BiCuSeO. As shown in Figure 18d, the mode Grüneisen parameters vary between 0.08 and 6.74. The average Grüneisen parameter is approximately 2.5, which is considerably larger than the experimental value of 1.5 estimated from the average sound velocity in polycrystalline BiCuSeO.219 It is believed that the presence of the heavy element Bi contributes to the large Grüneisen parameter due to its large valence shell formed by lone-pair electrons.56,219 As can be seen from Figure 18 d, the Grüneisen parameters of the lowfrequency modes are very large, which indicates the strong interaction between acoustic phonons and optical phonons.185 These phonon−phonon Umklapp processes can significantly scatter acoustic phonons carrying heat and reduce the lattice thermal conductivity. The anomalously high Grüneisen parameters of BiCuSeO also reflect its crystal structure, which contains distorted edge-sharing CuSe4 tetrahedra and a layered structure with zigzag geometry. This implies a soft lattice, which can scatter phonons easily. 5.3. Half-Heusler MgAgSb

Materials with the half-Heusler structure can possess excellent electrical transport properties that make them of strong interest for thermoelectric applications.220 MgAgSb is compositionally and structurally related to many other half-Heusler materials, whose thermoelectric properties along with ab initio determination of crystal structures were extensively investigated.221 MgAgSb exists in three different crystal structures at 300−700 K, taking the half-Heusler structure at 630−700 K (γ-phase with space group F4̅3m), a Cu2Sb-related structure at 560−630 K (β-phase with space group P4/nmm), and a tetragonal structure at 300−560 K (α-phase with space group I4̅c2). The low Seebeck coefficient of ∼70 μV/K and high electrical conductivity of ∼2.5× 104 S/cm indicate that β-MgAgSb shows a metallic transport behavior; therefore, most reports focus on the thermoelectric properties at 300−560 K of α-MgAgSb. As shown in Figure 19a, the α-MgAgSb phase consists of a distorted Mg−Sb rocksalt lattice, rotated by 45° about the c 12137

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

disordered throughout the structure, and the disordered Cu ions are superionic mobile.231 Figure 21b shows the temperature dependence of ZT values for copper chalcogenides. Yu et al. described Cu2Se as a thermoelectric material with an ordered selenium layer and a disordered copper layer, with a maximum ZT value of ∼1.6 at 973 K.229 Liu et al. described Cu2Se a copper ion liquid-like thermoelectric material with a maximum ZT value of ∼1.5 at 1000 K.36 Following these reports, Zhong et al.232 supposed that Cu1.94Al0.02Se would show anisotropic transport properties and created a highly aligned large lamellae in bulk Cu1.94Al0.02Se using the direct-current hot-pressing process. A high ZT of 2.62 at ∼1030 K was reported in Cu1.94Al0.02Se with lamellar structure, supposing that the superionic properties were concentrated on one direction, and a ZT of 1.4 at ∼1030 K along the direction perpendicular to lamellar structure was reported. Meanwhile, Zhong et al. mentioned that a big spread of ZT values among 50 samples raised concerns about repeatability.232 Surprisingly, Cu1.97S also exhibits an impressive high ZT of 1.7 at 1000 K,41 which was further enhanced to ZT of 2.1 at 1000 K in Cu2S0.52Te0.48.233 The high ZT obtained in Cu2S0.52Te0.48 is attributed to a special microstructure involving mosaic crystals which exhibit the good electrical properties of single crystals and the low thermal conductivity of polycrystalline or nanocrystalline samples. Unlike the random orientation in ordinary polycrystalline materials, blocks in a so-called mosaic crystal exhibit a nearly identical orientation; thereby the bulk appears like a single crystal from a macroscale point of view but contains a number of small angle boundaries. The apparent high performance in copper chalcogenides mainly comes from the intriguing thermal transport properties, especially in the thermal conductivity. As shown in Figure 22a,

Figure 20. (a) Thermal conductivity as a function of temperature for samples of MgAgSb.34,223 (b) Ag+ migration model from Ag site to Mg vacancy and (c) Mg2+ migration model from Mg site to Ag vacancy in α-MgAgSb.185 (b and c) Reprinted with permission from Li et al. Adv. Funct. Mater. 2015, 25, 6478−6488. Copyright 2015 WILEY Publishing.

Li et al. found that unique Ag and Mg ion migrations in αMgAg0.97Sb0.99 cause significant local atomic disorder which is believed to be the main origin of the low thermal conductivity.185 DFT calculations indicate that the activation energy for Ag+ migration from the Ag2 site to the Mg vacancy is ∼0.202 eV, as shown in Figure 20b. The respective activation for Mg2+ migration from the Mg site to the Ag2 vacancy is ∼0.485 eV, as shown in Figure 20c. This reveals a unique type of ion migration between the Ag site and the Mg site in αMgAg0.97Sb0.99, which weakens the chemical bonds, consistent with the softening transverse phonon modes, the low Debye temperature (160 K), and the low mean phonon velocity of 1276 ms−1. 5.4. Copper Chalcogenides

Copper sulfides and selenides have been recognized as interesting thermoelectric materials. The history of copper chalcogenides goes back 1827, when Becquerel observed that burning some sulfur powder atop the copper wire significantly enhanced the generation of electricity.225,226 They were also studied in the 1960s;227 then the general interest waned only to re-emerge in recent years. Ge et al. reported a maximum ZT value 0.5 at 673 K for polycrystalline Cu1.8S fabricated using a combined process of mechanical alloying and spark plasma sintering,228 and then Liu et al.36 and Yu et al.229 independently reported high ZT values of 1.5−1.6 at ∼1000 K for Cu2Se. Interested readers are encouraged to refer to a relevant review by Dennler et al.,186 summarizing the long history of copper chalcogenides as thermoelectric materials. Cu2Se undergoes a phase transition at ∼400 K, where during the cooling process β-Cu2Se turns into α-Cu2Se and this transition is reversible.230 Figure 21a shows the β-Cu2Se crystal structure, where the Se atoms are in a simple face-centeredcubic packing in the space group Fm3̅m, but the Cu ions are

Figure 22. (a) High-temperature thermal conductivity of Cu2Se as a function of temperature.36,41,229,232,233 (b) Specific heat capacity as a function of temperature for copper chalcogenides.36 Dashed blue line shows the expected value of the specific heat at constant pressure Cp in a solid crystal Cu2Se without liquid-like properties, which is usually beyond the Dulong−Petit limit at high temperatures owing to the extra contributions by carriers (Ce) and lattice thermal expansion. B is the bulk modulus, V is the volume per atom, and α is the thermal expansion coefficient.

the thermal conductivities are ultralow at high temperatures, namely, the minimum thermal conductivities at ∼1000 K can reach as low as ∼0.71, ∼0.60 , ∼0.48, and ∼0.39 W m−1 K−1 for Cu2Se,36,229 Cu1.94Al0.02Se,232 Cu1.97S,41 and Cu2S0.52Te0.48,233 respectively. Figure 22b shows that these materials exhibit an unusual heat capacity.36 The theoretical value (Dulong−Petit law) for the high-temperature specific heat at constant volume Cv is 3NkB in a solid crystal, where N is the number of particles and kB is the Boltzmann constant. The lowest Cv theoretical value in a liquid is 2NkB. Cu2Se shows a reduced Cv, approaching 2NkB at high temperatures, which is a typical ion liquid-like behavior,36 as in an ionic battery.234 However, the

Figure 21. (a) β-Cu2Se crystal structure: yellow, Cu atoms; red, Se atoms. (b) ZT values as a function of temperature for copper chalcogenides.36,41,229,232,233 12138

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

physical explanation of this ion liquid-like behavior phenomenon is not clear yet. In an ion liquid-like semiconductor, the lattice thermal conductivity can be reduced below that of a glass by reducing not only the mean free path of phonons but also eliminating lattice vibrations. The idea of using the liquid-like behavior of superionic conductors may be considered an extension of the phonon-glass electron-crystal concept, and such materials could be considered phonon-liquid electron-crystal (PLEC) thermoelectrics. The very property that makes copper chalcogenides exhibit good thermoelectric performance may also be its Achilles heel when it comes to moving these materials to device fabrication and ultimate technological applications, namely, the very low thermal conductivity achievable in these systems is due to the very high mobility of Cu ions in the lattice causing great disorder and strong phonon scattering. Because thermoelectric devices are low-voltage high-current devices operating in a dc model, massive and rapid migration of Cu ions from one end of the sample to the other and severe Se/S vaporization will cause decomposition and massive polarization issues rendering the materials unstable for application. The realization is likely the reason that these materials were not pursued further in the first phase of their investigation in the 1960s. Incidentally, in the photovoltaic field the same issue of strong Cu ion migration and rapid instability during device operation caused the promising and good photovoltaic properties of Cu2S to be abandoned. These issues have been reviewed and comprehensively discussed in recent publications.186,235

centers form Bi−Bi bonds of 3.238 Å. The compound has a lamellar structure with slabs of [Bi4Te6]− alternating with layers of Cs+ ions, Figure 23b. Complex bismuth chalcogenides show very low thermal conductivities. Compared to the state-of-the-art thermoelectric Bi−Sb−Te241 and Bi−Sb,242 it can be readily seen that the thermal conductivity of β-K2Bi8Se13 is low, ranging from 1.0 to 1.5 W m−1 K−1, at temperatures below room temperature. The thermal conductivity of CsBi4Te6 decreases with rising temperature and reaches ∼1.5 W m−1 K−1 at room temperature, Figure 24a. No reports are available that deal rigorously

Figure 24. (a) Thermal conductivity comparisons of complex bismuth chalcogenides, and current state-of-the-art low-temperature thermoelectrics.29,236,240−242 (b) ZT values comparisons of complex bismuth chalcogenides, and current state-of-the-art low-temperature thermoelectrics.29,236,240−242

with origins of the intrinsically low thermal conductivity of these complex bismuth chalcogenides. Clarke reported materials selection guidelines for low thermal conductivity thermal barrier coatings.243 According to those guidelines, a material will have a low thermal conductivity if it satisfies these principal conditions: A large molecular weight, a complex crystal structure, nondirectional bonding, and a large number of different atoms per molecule. Complex bismuth chalcogenides fulfill these conditions. Electrical transport properties and thermal conductivity measurements on β-K2Bi8Se13 single crystals showed a relatively high power factor (∼10 μW cm−1 K−2) at room temperature, giving ZT of 0.22 at room temperature,240 Figure 24b. Extensive doping experiments on CsBi4Te6 give high values of power factor (>30 μW cm−1 K−2) at 100−220 K. The highest power factors of 40−60 μW cm−1 K−2 at 150−180 K were obtained through doping with Sb, Bi, SbI3, and BiI3. A ZT value of 0.8 at 225 K was obtained for 0.06% SbI3-doped CsBi4Te6, which is record high at a temperature gap of 150−275 K.29

5.5. Complex Bismuth Chalcogenides

Since the solid solutions of Bi−Sb−Te−Se were established as the leading thermoelectric materials near room temperature,66,79 there have been continuing efforts to find better thermoelectric materials with these elements, and ternary alkali metal−bismuth−chalcogenides are typical ones that have been found.29,236−240 Compared to Bi−Sb−Te−Se solid solutions, ternary alkali metal−bismuth−chalcogenides possess more complex crystal structures such as β-K2Bi8Se13 and CsBi4Te6. As shown in Figure 23a, β-K2Bi8Se13 has a low-symmetry

5.6. Chalcoantimonates with Lone-Pair Electrons

Just as the so-called lone pair of electrons in Pb2+, Sn2+, Ge2+, Bi3+ and Sb3+ plays a role in defining the thermoelectric properties of the various materials described above (and particularly the thermal conductivity), so do they define the properties of the compounds in the Cu/Sb/Se system. In this class we have three members which allow one to more clearly appreciate the role of the lone pair of electrons on the thermal conductivity, Cu3SbSe4, Cu3SnSe3, and CuSbSe2. The first member is a Sb5+ species and has no lone pair of electrons, while the latter two contain Sb3+ with an active lone pairs of electrons. As shown in Figure 25a, Cu3SbSe4 shows a tetragonal structure with the space group I42m. It derives from the diamond structure in which all atoms have tetrahedral coordination geometry.244 CuSbSe2 has an orthorhombic structure with the space group Pnma190,245 and Cu3SbSe3 a

Figure 23. (a) β-K2Bi8Se13 crystal structure: blue, Bi atoms; pink, Se atoms; gray, Ag atoms. Red box represents a unit cell. (b) Typical CsBi4Te6 crystal structure: blue, Bi atoms; pink, Cs atoms; dark reddish, Te atoms. Red box represents a unit cell.

monoclinic structure that consists of two Bi/Se building blocks (NaCl type and Bi2Se3 type) connected at the K/Bi mixed site (blue); the K atoms are in the channels. These two Bi/Se blocks are infinitely extended along the crystallographic b axis and connected to each other at special mixed-occupancy K/Bi sites. It has a highly anisotropic structure, and it grows along the b axis with a needle or columnar structure.240 CsBi4Te6 features a strong anisotropic structure that contains both formally Bi3+ and Bi2+centers.29 The unusual reduced Bi2+ 12139

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

Figure 25. Typical crystal structures for (a) Cu3SbSe4, (b) CuSbSe2, and (c) Cu3SbSe3: yellow, Cu atoms; red, Se atoms; gray, Sb atoms. Schematic representation of the local atomic environment of Sb in (d) Cu3SbSe4, (e) CuSbSe2, and (f) Cu3SbSe3. Shaded lines represent Sb− Se bonds; dashed lines illustrate the approximate morphology of the Sb lone-pair 5s electron orbital.

Figure 26. (a) Temperature dependence of the thermal conductivity of Cu3SbSe4, CuSbSe2, and Cu3SbSe3.190 (b) Temperature dependence of the thermal conductivity of two I−V−VI2 cubic semiconductors AgSbTe2 and AgBiSe2, of single-crystal PbTe, and of the chalcopyrite compound AgInTe2.191 Electronic contributions to the thermal conductivity are negligible in all samples; thus, the data shown represent the lattice thermal conductivity. Solid line represents the calculated minimum thermal conductivity for AgSbTe2. (c) Temperature dependence of ZT values for Cu3SbSe4 and Cu3SbSe3.249−251 (d) Temperature dependence of ZT values for I−V−VI2 cubic semiconductors.233,253−257

different orthorhombic structure with the space group Pnma,244 Figure 25b and 25c. In Cu3SbSe4, the Sb5+ centers are coordinated by 4 Se atoms with tetrahedral Se−Sb−Se bond angles of 109.5°, Figure 25d, suggesting sp3 hybridization of the Sb valence electron orbitals. In CuSbSe2, however, Sb is in a 3+ formal oxidation state and coordinated by 3 Se atoms in a trigonal pyramidal geometry and an average Se−Sb−Se bond angle of 95.24°, as shown in Figure 25e. In this arrangement only the Sb 5p electrons form bonds with Se, leaving the Sb 5s2 electrons “free” to orient along the missing vertex of the tetrahedron. The configuration is similar for Cu3SbSe3 except the average Se−Sb−Se bond angle is 99.42°. Once again the Sb 5s2 lone pair helps the Sb3+ to form an imperfect tetrahedron if we consider that is occupies the fourth coordination site, Figure 25f. From a structural point of view, one may expect these three compounds to have similar thermal conductivity since their average atomic masses are nearly the same and none of them possess an overly complex crystal structure; this however is not the case. As shown in Figure 26a, Cu3SbSe4, CuSbSe2, and Cu3SbSe3 possess vastly different thermal conductivity.190 In Cu3SbSe4, the thermal conductivity decreases rapidly with rising temperature due to strong phonon−phonon interactions (typical for a crystalline material). In CuSbSe2, thermal conductivity shows the same decreasing trend with temperature, but it is only one-half that of Cu3SbSe4. Remarkably, Cu3SbSe3 has a very low and temperature-independent thermal conductivity even at cryogenic temperatures. The reason for the large thermal conductivity differences is the different intrinsic phonon scattering mechanisms acting in Cu3SbSe4, CuSbSe2, and Cu3SbSe3. These mechanisms are related to the role of the lone pair of s electrons in the Sb3+ centers in CuSbSe2 to Cu3SbSe3. The nature of the lone pair in heavy main group atoms and its chemical and physical properties have been the subject of much study and debate over many decades and remains topical today.246−248 We can say that understanding the behavior of these lone pairs and how they influence the electronic properties is still incomplete. The relationship between lonepair electrons and low thermal conductivity involves the enhancement of anharmonicity in the lattice by the overlapping

wave functions of the lone-pair electrons and nearby valence electrons which induces nonlinear repulsive electrostatic forces. This is well supported by the density functional theory calculations by Zhang et al.,245 showing that Cu3SbSe3 has larger Grüneisen parameters, smaller Debye temperatures, and lower phonon velocities than Cu3SbSe4. The highest degree of anharmonicity should thus be achieved when the lone-pair electron is far removed from the Sb nucleus and not participating in bonding. Therefore, it is reasonable that the thermal conductivity of CuSbSe2 is higher than that of Cu3SbSe3, since the lone electron pair is farther removed from the Sb nucleus in Cu3SbSe3. Another interesting comparison is the case of AgInTe2 versus AgSbTe2,191 where the thermal conductivity of AgInTe2 (lacking lone pair of electrons) is “normal” while that of AgSbTe2 (Sb3+ 5s2 pair) is glasslike, approaching the calculated minimum thermal conductivity value, Figure 26b. Wei et al.249,250 extensively investigated the electrical and thermal transport properties of Cu3SbSe4 between 300 and 673 K. They found that the maximum ZT at 673 K increases from 0.4 for undoped Cu3SbSe4 to 0.7 when Sn doping is used on the Sb site, Figure 26c. For Cu3SbSe3 an extremely low thermal conductivity (0.25−0.30 W m−1 K−1) was observed, a maximum ZT of 0.25 was obtained at 650 K.251 AgSbTe 2 is well known for its very low thermal conductivity.190,191,252 Wang et al.253 investigated the polycrystalline AgSbTe2, which was fabricated using a combined process of mechanical alloying and spark plasma sintering. It was found that stoichiometric AgSbTe2 is a promising thermoelectric material, whose ZT reaches 1.6 at 673 K, benefiting from its extremely low thermal conductivity of 0.30 W m−1 K−1, Figure 26 d. Nielsen et al.189 reported that lonepair electrons also could minimize the lattice thermal conductivity of AgSbSe2, as in AgSbTe2, since the Grüneisen parameter (3.5) of AgSbSe2 is much larger than that (2.05) of 12140

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

AgSbTe2; however, it was found that the electrical transport properties of AgSbSe2 are inferior to AgSbTe2. By improving the electrical transport properties, the ZT values have been enhanced from ∼0.4 for undoped AgSbSe2 to ∼1.0 through introducing Sb vacancies,254 Na doping,233 and Pb doping.255 As an analogue of AgSbSe2, AgBiSe2 is a n-type semiconductor and exhibits very low thermal conductivity; ZT values could be achieved in 0.9−1.0 through Nb doping on the Ag site256 and I doping on the Se site.257

Figure 28. (a) Thermal conductivity as a function of temperature for samples of tetrahedrite Cu12Sb4S13. (b) Sb[CuS3]Sb trigonal bipyramid. (c) Valence electron density map projected onto a (101) plane showing the coordination of Cu12e, Sb, S2a, and S24g atoms. Space group of Cu12Sb4S13 is I-43m and features two Cu sites (12d and 12e), two S sites (2a and 24g), and one Sb site (8c), blue for Cu12e, green for S2a, yellow for S24g, and brown for Sb.251 (b and c) Reprinted with permission from ref 251. Copyright 2015 WILEY Publishing.

5.7. Tetrahedrites

Recently, some natural mineral tetrahedrites were reported as a possible direct source of thermoelectric materials, led by the independent groups of Morelli and Suekuni.258−260 The tetrahedrite phase with general formula Cu12Sb4S13 is one of the most widespread sulfosalt minerals on Earth. It can accommodate a variety of substitutions onto Cu and Sb sites; thus, natural tetrahedrite is commonly written as (Cu, Ag)10(Cu, Zn, Fe, Co, Ni, Cd, Hg)2(Sb, Te, Bi, As)4(S, Se)13. Interest in studying the thermoelectric properties of tetrahedrites has surged because they consist of environmentally benign elements copper and sulfur.261 Figure 27a shows the

drites may result from the synergistic effects of local bonding asymmetry and anharmonic rattling modes.251 As shown in Figure 28 b, a strong local bonding asymmetry is identified inside a Sb[CuS3]Sb atomic cage, in which Cu (blue) forms strong covalent bonding with S (green and yellow) and weak covalent bonding with one of the Sb atoms (brown). The bonding asymmetry of Cu12e causes the low thermal conductivity of tetrahedrites, as indicated in Figure 28c. In addition, the localized anharmonic rattling modes enabled by the lone-pair electrons of Sb also can contribute to resonant scattering of phonons, which is the as same to that of pnicogengroup chalcogenides with lone-pair electrons. Skoug and Morelli identified a correlation between minimal thermal conductivity and the existence of an Sb lone s2 pair in Sbcontaining ternary semiconductors.190 Lone-pair electrons of Sb induce strong lattice anharmonicity;191 the “crescent-moon” shape around Sb is a clear indicator of the existence of electron lone pairs, as shown in Figure 28c.

Figure 27. (a) Crystal structure of tetrahedrite Cu12Sb4S13; (b) ZT values as a function of temperature for doped samples of tetrahedrite Cu12Sb4S13.258,262−265

crystal structure of tetrahedrite Cu12Sb4S13. It possesses a cubic sphalerite-like structure of I-43m symmetry with 6 of the 12 Cu atoms occupying trigonal planar 12e sites and the remaining Cu atoms distributed on tetrahedral 12d sites. The Sb atoms also occupy a tetrahedral site but are bonded to only three S atoms, leading to a void in the structure presumably occupied by the lone pair of electrons. In terms of a simple crystal chemical formula, 4 of the 6 tetrahedral sites are thought to be occupied by monovalent Cu while the other 2 are occupied by Cu2+ ions; the trigonal planar sites are occupied solely by monovalent Cu. Lu et al. found that pure synthetic Cu12Sb4S13 exhibits a ZT value of 0.56 at 673 K,258 as shown in Figure 27b. After optimizing and substituting on Cu or Sb sites, the ZT at 723 K could be further increased to 0.76 for Cu10.5Ni1.5Sb4S13, 0.90 for Cu12Sb3Te4S13, 0.92 for Cu11.5Zn0.5Sb4S13, and 1.1 for Cu10.5Ni1.0Zn0.5Sb4S13.262−265 The promising thermoelectric performance of Cu12Sb4S13 is ascribed to the intrinsically low thermal conductivity. As shown in Figure 28a, the thermal conductivity of Cu12Sb4S13 shows an increasing trend up to 723 K. After doping and alloying, the thermal conductivities are reduced in the entire temperature range; some samples even show the glass-like and temperatureindependent thermal transport behaviors. It seems that the phonon mean free path approaches one interatomic spacing, and the thermal transport reaches the “minimal” thermal conductivity, namely, the thermal conductivity ranges from 0.40 W m−1 K−1 at 300 K to 0.45 W m−1 K−1 at 723 K for Cu10.5Ni1.5Sb4S13.265 The low thermal conductivity of tetrahe-

6. SUMMARY AND OUTLOOK The past two decades have witnessed the rapid growth and several exciting conceptual and performance breakthroughs in thermoelectric research. The development of better and more accessible physical property and thermoelectric measurement techniques and commercially available robust instrumentation especially for reliable high-temperature measurements have played an important role in facilitating in the growth of the thermoelectric community. Furthermore, the power of synthetic solid state chemistry in its ability to deliver new materials and new processing techniques has been critical in this regard. New concepts make it possible to invent new strategies on how to decouple the strong interrelation of thermoelectric transport parameters, raising the ZT values steadily from the historical threshold of around unity to well above 2. As there are no theoretical or thermodynamic limits to the maximum value of ZT, much higher performance can be anticipated for the next generation of thermoelectric devices. Given the many materials yet to be investigated, there is certainly much work ahead and promise for developing higher efficiency thermoelectric materials and devices. As we discussed above, there are generally two main directions to achieving high-performance thermoelectrics: (1) optimizing the known materials by band structure and microstructure manipulation and (2) exploring new compounds with peculiar physical properties that may lead to 12141

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

unexpectedly high ZTs. For the former it is highly desirable to integrate all advanced performance-enhancing concepts in one single material. The latter is also appealing given the power of intuitive thinking63 and the elements of surprising discoveries. In principle, the materials selection process can be accelerated using high-throughput computational techniques to perform initial screenings of the search space and identify the best candidates to be investigated experimentally, namely, the socalled materials genome approach. The materials with intrinsically very low thermal conductivity are also of strong interest, especially when they simultaneously feature unique band structures (for example, the conducting bands are composed of multicarrier pockets) that allow high power factors. Of course, considering the complexity of thermoelectric materials, their future development would benefit greatly from close collaborations between chemists, physicists, and materials scientists. Given the need for alternative energy technologies and materials to ultimately replace the shrinking supply of fossil fuels or extend its duration, energy-related research and the study of clean and reliable thermoelectric energy conversion technology will play an important role.

postdoctoral research fellow at the University of Michigan and Northwestern University. He holds a Charles E. and Emma H. Morrison Professor Chair at Northwestern University.

ACKNOWLEDGMENTS This work was supported in part by a Department of Energy grant to the Revolutionary Materials for Solid State Energy Conversion, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, and Office of Basic Energy Sciences under Award Number DE-SC0001054. Currently, work at Northwestern is supported by the Department of Energy, Office of Science Basic Energy Sciences grant DE-SC0014520. We also acknowledge the use of the transmission electron microscope from the NUANCE Center at Northwestern University. This work was also supported by the “Zhuoyue” program of Beihang University, the Recruitment Program for Young Professionals, and NSFC under Grant No. 51571007 (L.D.Z). We thank especially Professors C. Uher, C. Wolverton, and V. P. Dravid, whose collaboration has been critical to our understanding of many aspects of our novel materials in our laboratory. We also thank Professors D. N. Seidman, T. P. Hogan, J. P. Heremans, E. D. Case, X. F. Tang, H. B. Xu, S. K. Gong, J. Q. He, and Y. L. Pei for plentiful discussions and fruitful collaborations. Of course, most of all, we are grateful to the numerous dedicated graduate students and postdoctoral fellows who have contributed to our thermoelectric research efforts. Their names appear in the various publications cited in this article.

AUTHOR INFORMATION Corresponding Authors

*E-mail: [email protected]. *E-mail: [email protected]. Notes

The authors declare no competing financial interest.

REFERENCES

Biographies

(1) Hsu, K. F.; Loo, S.; Guo, F.; Chen, W.; Dyck, J. S.; Uher, C.; Hogan, T.; Polychroniadis, E.; Kanatzidis, M. G. Cubic AgPbmSbTe2+m: Bulk Thermoelectric Materials with High Figure of Merit. Science 2004, 303, 818−821. (2) Poudel, B.; Hao, Q.; Ma, Y.; Lan, Y.; Minnich, A.; Yu, B.; Yan, X.; Wang, D.; Muto, A.; Vashaee, D.; et al. High-Thermoelectric Performance of Nanostructured Bismuth Antimony Telluride Bulk Alloys. Science 2008, 320, 634−638. (3) 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. (4) 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. (5) Vineis, C. J.; Shakouri, A.; Majumdar, A.; Kanatzidis, M. G. Nanostructured Thermoelectrics: Big Efficiency Gains from Small Features. Adv. Mater. 2010, 22, 3970−3980. (6) Joshi, G.; Lee, H.; Lan, Y.; Wang, X.; Zhu, G.; Wang, D.; Gould, R. W.; Cuff, D. C.; Tang, M. Y.; Dresselhaus, M. S.; et al. Enhanced Thermoelectric Figure-of-Merit in Nanostructured P-Type Silicon Germanium Bulk Alloys. Nano Lett. 2008, 8, 4670−4674. (7) Bux, S. K.; Blair, R. G.; Gogna, P. K.; Lee, H.; Chen, G.; Dresselhaus, M. S.; Kaner, R. B.; Fleurial, J. P. Nanostructured Bulk Silicon as an Effective Thermoelectric Material. Adv. Funct. Mater. 2009, 19, 2445−2452. (8) Boydston, R. W. Thermo-Electric Effect in Single-Crystal Bismuth. Phys. Rev. 1927, 30, 911−921. (9) Henkels, H. W. Thermoelectric Power and Mobility of Carriers in Selenium. Phys. Rev. 1950, 77, 734−736. (10) Soroos, A. Thermoelectric Power of Single Crystal Bismuth near the Melting Point. Phys. Rev. 1932, 41, 516−522. (11) Bidwell, C. C. Electrical Resistance and Thermo-Electric Power of the Alkali Metals. Phys. Rev. 1924, 23, 357−376. (12) Greenwood, N.; Anderson, J. Conductivity and Thermoelectric Effect in Cuprous Oxide. Nature 1949, 164, 346−347.

Gangjian Tan received in 2013 his Ph.D. degree in Materials Science from Wuhan University of Technology, China, under the supervision of Professor Xinfeng Tang, with a focus on nanostructured skutterudite compounds as advanced thermoelectric materials by melt spinning. In June of 2013, he joined Professor Mercouri Kanatzidis’s group at Northwestern University as a postdoctoral research fellow, where his research extends to the development of new thermoelectric chalcogenides. He is interested in the understanding of charge and phonon transport behavior in solids associated with thermoelectric effects. Li-Dong Zhao is a Professor of Materials Science and Engineering at Beihang University, China. He received his B.E. and M.E. degrees in Materials Science from the Liaoning Technical University and his Ph.D. degree in Materials Science from the University of Science and Technology Beijing, China, in 2009. He was a postdoctoral research fellow in the LEMHE-ICMMO (CNRS-UMR 8182) at the University of Paris-Sud from 2009 to 2011 and continued as a postdoctoral research fellow in the Department of Chemistry at Northwestern University from 2011 to 2014. He holds the fellowship of “Recruitment Program for Young Professionals”. His research interests include thermoelectric materials, superconductors, and thermal barrier coatings. Mercouri G. Kanatzidis is a Professor of Chemistry and of Materials Science and Engineering at Northwestern University in Evanston, IL. He also has a senior scientist appointment at Argonne National Laboratory. His interests include the design and synthesis of new materials with emphasis on systems with highly unusual structural/ physical characteristics or those capable of energy conversion, energy detection, environmental remediation, and catalysis. After obtaining his B.Sc. degree from Aristotle University in Greece, he received his Ph.D. degree in Chemistry from the University of Iowa and was a 12142

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

(13) Andrews, J. Thermoelectric Power of Cadmium Oxide. Proc. Phys. Soc. 1947, 59, 990−998. (14) Bidwell, C. C. Electrical and Thermal Properties of Iron Oxide. Phys. Rev. 1917, 10, 756−766. (15) Goldsmid, H.; Douglas, R. The Use of Semiconductors in Thermoelectric Refrigeration. Br. J. Appl. Phys. 1954, 5, 386−390. (16) Herring, C. Theory of the Thermoelectric Power of Semiconductors. Phys. Rev. 1954, 96, 1163−1187. (17) Price, P. Theory of Transport Effects in Semiconductors: Thermoelectricity. Phys. Rev. 1956, 104, 1223−1239. (18) Goldsmid, H. Thermoelectric Applications of Semiconductors. J. Electron. Control 1955, 1, 218−222. (19) Hicks, L.; Dresselhaus, M. Effect of Quantum-Well Structures on the Thermoelectric Figure of Merit. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 47, 12727−12731. (20) Dresselhaus, M. S.; Chen, G.; Tang, M. Y.; Yang, R.; Lee, H.; Wang, D.; Ren, Z.; Fleurial, J. P.; Gogna, P. New Directions for Low Dimensional Thermoelectric Materials. Adv. Mater. 2007, 19, 1043− 1053. (21) Hicks, L.; Dresselhaus, M. Thermoelectric Figure of Merit of a One-Dimensional Conductor. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 47, 16631−16634. (22) Kim, S. I.; Lee, K. H.; Mun, H. A.; Kim, H. S.; Hwang, S. W.; Roh, J. W.; Yang, D. J.; Shin, W. H.; Li, X. S.; Lee, Y. H.; et al. Dense Dislocation Arrays Embedded in Grain Boundaries for HighPerformance Bulk Thermoelectrics. Science 2015, 348, 109−114. (23) 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. (24) Liu, W.; Tan, X.; Yin, K.; Liu, H.; Tang, X.; Shi, J.; Zhang, Q.; Uher, C. Convergence of Conduction Bands as a Means of Enhancing Thermoelectric Performance of N-Type Mg2Si1−xSnx Solid Solutions. Phys. Rev. Lett. 2012, 108, 166601. (25) 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. (26) Pei, Y.; Wang, H.; Snyder, G. Band Engineering of Thermoelectric Materials. Adv. Mater. 2012, 24, 6125−6135. (27) Zhao, L.-D.; He, J.; Hao, S.; Wu, C.-I.; Hogan, T. P.; Wolverton, C.; Dravid, V. P.; Kanatzidis, M. G. Raising the Thermoelectric Performance of P-Type PbS with Endotaxial Nanostructuring and Valence-Band Offset Engineering Using CdS and ZnS. J. Am. Chem. Soc. 2012, 134, 16327−16336. (28) Zhao, L.-D.; Hao, S.; Lo, S.-H.; Wu, C.-I.; Zhou, X.; Lee, Y.; Li, H.; Biswas, K.; Hogan, T. P.; Uher, C.; et al. High Thermoelectric Performance Via Hierarchical Compositionally Alloyed Nanostructures. J. Am. Chem. Soc. 2013, 135, 7364−7370. (29) Chung, D.-Y.; Hogan, T.; Brazis, P.; Rocci-Lane, M.; Kannewurf, C.; Bastea, M.; Uher, C.; Kanatzidis, M. G. CsBi4Te6: A HighPerformance Thermoelectric Material for Low-Temperature Applications. Science 2000, 287, 1024−1027. (30) Saramat, A.; Svensson, G.; Palmqvist, A.; Stiewe, C.; Mueller, E.; Platzek, D.; Williams, S.; Rowe, D.; Bryan, J.; Stucky, G. Large Thermoelectric Figure of Merit at High Temperature in CzochralskiGrown Clathrate Ba8Ga16Ge30. J. Appl. Phys. 2006, 99, 023708. (31) Rhyee, J.-S.; Lee, K. H.; Lee, S. M.; Cho, E.; Kim, S. I.; Lee, E.; Kwon, Y. S.; Shim, J. H.; Kotliar, G. Peierls Distortion as a Route to High Thermoelectric Performance in In4Se3‑δ Crystals. Nature 2009, 459, 965−968. (32) Wang, S.; Li, H.; Lu, R.; Zheng, G.; Tang, X. Metal Nanoparticle Decorated N-Type Bi2Te3-Based Materials with Enhanced Thermoelectric Performances. Nanotechnology 2013, 24, 285702. (33) Shi, X.; Yang, J.; Salvador, J. R.; Chi, M.; Cho, J. Y.; Wang, H.; Bai, S.; Yang, J.; Zhang, W.; Chen, L. Multiple-Filled Skutterudites: High Thermoelectric Figure of Merit through Separately Optimizing Electrical and Thermal Transports. J. Am. Chem. Soc. 2011, 133, 7837−7846.

(34) Zhao, H.; Sui, J.; Tang, Z.; Lan, Y.; Jie, Q.; Kraemer, D.; McEnaney, K.; Guloy, A.; Chen, G.; Ren, Z. High Thermoelectric Performance of MgAgSb-Based Materials. Nano Energy 2014, 7, 97− 103. (35) Rogl, G.; Grytsiv, A.; Bauer, E.; Rogl, P.; Zehetbauer, M. Thermoelectric Properties of Novel Skutterudites with Didymium: DDy(Fe1−xCox)4Sb12 and DDy(Fe1−xNix)4Sb12. Intermetallics 2010, 18, 57−64. (36) 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. (37) Yamada, T.; Yamane, H.; Nagai, H. A Thermoelectric Zintl Phase Na2+xGa2+XSn4−x with Disordered Na Atoms in Helical Tunnels. Adv. Mater. 2015, 27, 4708−4713. (38) Wu, D.; Zhao, L.-D.; Hao, S.; Jiang, Q.; Zheng, F.; Doak, J. W.; Wu, H.; Chi, H.; Gelbstein, Y.; Uher, C.; et al. Origin of the High Performance in GeTe-Based Thermoelectric Materials Upon Bi2Te3 Doping. J. Am. Chem. Soc. 2014, 136, 11412−11419. (39) Korkosz, R. J.; Chasapis, T. C.; Lo, S.-h.; Doak, J. W.; Kim, Y. J.; Wu, C.-I.; Hatzikraniotis, E.; Hogan, T. P.; Seidman, D. N.; Wolverton, C.; et al. High ZT in P-Type (PbTe)1−2x(PbSe)x(PbS)x Thermoelectric Materials. J. Am. Chem. Soc. 2014, 136, 3225−3237. (40) 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. (41) He, Y.; Day, T.; Zhang, T.; Liu, H.; Shi, X.; Chen, L.; Snyder, G. J. High Thermoelectric Performance in Non - Toxic Earth - Abundant Copper Sulfide. Adv. Mater. 2014, 26, 3974−3978. (42) 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 Pseudo-Binary: Extremely Low Thermal Conductivity and Modulated Carrier Concentration. Energy Environ. Sci. 2015, 8, 2056−2068. (43) Tan, G.; Shi, F.; Hao, S.; Chi, H.; Zhao, L.-D.; Uher, C.; Wolverton, C.; Dravid, V. P.; Kanatzidis, M. G. Codoping in SnTe: Enhancement of Thermoelectric Performance through Synergy of Resonance Levels and Band Convergence. J. Am. Chem. Soc. 2015, 137, 5100−5112. (44) Lee, Y.; Lo, S.-H.; Chen, C.; Sun, H.; Chung, D.-Y.; Chasapis, T. C.; Uher, C.; Dravid, V. P.; Kanatzidis, M. G. Contrasting Role of Antimony and Bismuth Dopants on the Thermoelectric Performance of Lead Selenide. Nat. Commun. 2014, 5, 3640. (45) Wood, C. Materials for Thermoelectric Energy Conversion. Rep. Prog. Phys. 1988, 51, 459−539. (46) DiSalvo, F. J. Thermoelectric Cooling and Power Generation. Science 1999, 285, 703−706. (47) Snyder, G. J.; Toberer, E. S. Complex Thermoelectric Materials. Nat. Mater. 2008, 7, 105−114. (48) Sootsman, J. R.; Chung, D. Y.; Kanatzidis, M. G. New and Old Concepts in Thermoelectric Materials. Angew. Chem., Int. Ed. 2009, 48, 8616−8639. (49) Kanatzidis, M. G. Nanostructured Thermoelectrics: The New Paradigm? Chem. Mater. 2010, 22, 648−659. (50) Li, J.-F.; Liu, W.-S.; Zhao, L.-D.; Zhou, M. High-Performance Nanostructured Thermoelectric Materials. NPG Asia Mater. 2010, 2, 152−158. (51) Heremans, J. P.; Wiendlocha, B.; Chamoire, A. M. Resonant Levels in Bulk Thermoelectric Semiconductors. Energy Environ. Sci. 2012, 5, 5510−5530. (52) Liu, W.; Yan, X.; Chen, G.; Ren, Z. Recent Advances in Thermoelectric Nanocomposites. Nano Energy 2012, 1, 42−56. (53) Minnich, A.; Dresselhaus, M.; Ren, Z.; Chen, G. Bulk Nanostructured Thermoelectric Materials: Current Research and Future Prospects. Energy Environ. Sci. 2009, 2, 466−479. (54) Kleinke, H. New Bulk Materials for Thermoelectric Power Generation: Clathrates and Complex Antimonides. Chem. Mater. 2010, 22, 604−611. 12143

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

(55) Toberer, E. S.; May, A. F.; Snyder, G. J. Zintl Chemistry for Designing High Efficiency Thermoelectric Materials. Chem. Mater. 2010, 22, 624−634. (56) Zhao, L.-D.; He, J.; Berardan, D.; Lin, Y.; Li, J.-F.; Nan, C.-W.; Dragoe, N. Bicuseo Oxyselenides: New Promising Thermoelectric Materials. Energy Environ. Sci. 2014, 7, 2900−2924. (57) Koumoto, K.; Wang, Y.; Zhang, R.; Kosuga, A.; Funahashi, R. Oxide Thermoelectric Materials: A Nanostructuring Approach. Annu. Rev. Mater. Res. 2010, 40, 363−394. (58) Zhang, Q.; Sun, Y.; Xu, W.; Zhu, D. Organic Thermoelectric Materials: Emerging Green Energy Materials Converting Heat to Electricity Directly and Efficiently. Adv. Mater. 2014, 26, 6829−6851. (59) Bubnova, O.; Crispin, X. Towards Polymer-Based Organic Thermoelectric Generators. Energy Environ. Sci. 2012, 5, 9345−9362. (60) Medlin, D.; Snyder, G. Interfaces in Bulk Thermoelectric Materials: A Review for Current Opinion in Colloid and Interface Science. Curr. Opin. Colloid Interface Sci. 2009, 14, 226−235. (61) Gaultois, M. W.; Sparks, T. D.; Borg, C. K.; Seshadri, R.; Bonificio, W. D.; Clarke, D. R. Data-Driven Review of Thermoelectric Materials: Performance and Resource Considerations. Chem. Mater. 2013, 25, 2911−2920. (62) Dehkordi, A. M.; Zebarjadi, M.; He, J.; Tritt, T. M. Thermoelectric Power Factor: Enhancement Mechanisms and Strategies for Higher Performance Thermoelectric Materials. Mater. Sci. Eng., R 2015, 97, 1−22. (63) Zeier, W. G.; Zevalkink, A.; Gibbs, Z. M.; Hautier, G.; Kanatzidis, M. G.; Snyder, G. J. Thinking Like a Chemist: Intuition in Thermoelectric Materials. Angew. Chem., Int. Ed. 2016, 55, 6826− 6841. (64) Cutler, M.; Leavy, J.; Fitzpatrick, R. Electronic Transport in Semimetallic Cerium Sulfide. Phys. Rev. 1964, 133, A1143−A1152. (65) Franz, R.; Wiedemann, G. Ueber Die Wärme - Leitungsfähigkeit Der Metalle. Ann. Phys. 1853, 165, 497−531. (66) Ioffe, A. F. Semiconductor Thermoelements and Thermoelectric Cooling; Infoserch Limited: London, 1957. (67) Dashevsky, Z.; Shusterman, S.; Dariel, M.; Drabkin, I. Thermoelectric Efficiency in Graded Indium-Doped PbTe Crystals. J. Appl. Phys. 2002, 92, 1425−1430. (68) Kuznetsov, V.; Kuznetsova, L.; Kaliazin, A.; Rowe, D. High Performance Functionally Graded and Segmented Bi2Te3-Based Materials for Thermoelectric Power Generation. J. Mater. Sci. 2002, 37, 2893−2897. (69) Dashevsky, Z.; Gelbstein, Y.; Edry, I.; Drabkin, I.; Dariel, M. Twenty-Second International Conference on Thermoelectrics, 2003; pp 421−424. (70) 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. (71) 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. (72) Yamini, S. A.; Ikeda, T.; Lalonde, A.; Pei, Y.; Dou, S. X.; Snyder, G. J. Rational Design of P-Type Thermoelectric PbTe: Temperature Dependent Sodium Solubility. J. Mater. Chem. A 2013, 1, 8725−8730. (73) Bergum, K.; Ikeda, T.; Snyder, G. J. Solubility and Microstructure in the Pseudo-Binary PbTe−Ag2Te System. J. Solid State Chem. 2011, 184, 2543−2552. (74) Brebrick, R. F.; Gubner, E. Composition Stability Limits of PbTe. J. Chem. Phys. 1962, 36, 1283−1289. (75) 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. (76) Yamini, S. A.; Mitchell, D. R. G.; Gibbs, Z. M.; Santos, R.; Patterson, V.; Li, S.; Pei, Y. Z.; Dou, S. X.; Jeffrey Snyder, G. Heterogeneous Distribution of Sodium for High Thermoelectric Performance of P-Type Multiphase Lead-Chalcogenides. Adv. Energy Mater. 2015, 5, 1501047.

(77) Chasmar, R.; Stratton, R. The Thermoelectric Figure of Merit and Its Relation to Thermoelectric Generators. J. Electron. Control 1959, 7, 52−72. (78) Tritt, T.; Kanatzidis, M.; Lyon, H.; Mahan, G. Thermoelectric MaterialsNew Directions and Approaches. Mater. Res. Soc. Proc. 1997, 478, 73−84. (79) Rowe, D. M. CRC Handbook of Thermoelectrics; CRC Press, 1995. (80) Goldsmid, H. Thermoelectric Refrigeration; Plenum Press: New York, 1964. (81) Nimtz, G.; Schlicht, B. Narrow-Gap Semiconductors; Springer, 1983. (82) Zhang, J.; Liu, R.; Cheng, N.; Zhang, Y.; Yang, J.; Uher, C.; Shi, X.; Chen, L.; Zhang, W. High-Performance Pseudocubic Thermoelectric Materials from Non - Cubic Chalcopyrite Compounds. Adv. Mater. 2014, 26, 3848−3853. (83) Zhao, L.; Wu, H.; Hao, S.; Wu, C.-I.; Zhou, X.; Biswas, K.; He, J.; Hogan, T. P.; Uher, C.; Wolverton, C.; et al. All-Scale Hierarchical Thermoelectrics: MgTe in PbTe Facilitates Valence Band Convergence and Suppresses Bipolar Thermal Transport for High Performance. Energy Environ. Sci. 2013, 6, 3346−3355. (84) Pei, Y.; LaLonde, A. D.; Heinz, N. A.; Snyder, G. J. High Thermoelectric Figure of Merit in PbTe Alloys Demonstrated in PbTe-CdTe. Adv. Energy Mater. 2012, 2, 670−675. (85) 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. (86) Pei, Y.; Wang, H.; Gibbs, Z. M.; LaLonde, A. D.; Snyder, G. J. Thermopower Enhancement in Pb1‑xMnxTe Alloys and Its Effect on Thermoelectric Efficiency. NPG Asia Mater. 2012, 4, e28. (87) Wang, H.; Gibbs, Z. M.; Takagiwa, Y.; Snyder, G. J. Tuning Bands of PbSe for Better Thermoelectric Efficiency. Energy Environ. Sci. 2014, 7, 804−811. (88) Banik, A.; Shenoy, U. S.; Anand, S.; Waghmare, U. V.; Biswas, K. Mg Alloying in Snte Facilitates Valence Band Convergence and Optimizes Thermoelectric Properties. Chem. Mater. 2015, 27, 581− 587. (89) Tan, G.; Shi, F.; Hao, S.; Chi, H.; Bailey, T. P.; Zhao, L.-D.; Uher, C.; Wolverton, C.; Dravid, V. P.; Kanatzidis, M. G. Valence Band Modification and High Thermoelectric Performance in SnTe Heavily Alloyed with MnTe. J. Am. Chem. Soc. 2015, 137, 11507− 11516. (90) Tan, G.; Shi, F.; Doak, J. W.; Sun, H.; Zhao, L.-D.; Wang, P.; Uher, C.; Wolverton, C.; Dravid, V. P.; Kanatzidis, M. G. Extraordinary Role of Hg in Enhancing the Thermoelectric Performance of P-Type SnTe. Energy Environ. Sci. 2015, 8, 267−277. (91) Zhang, Q.; Liao, B.; Lan, Y.; Lukas, K.; Liu, W.; Esfarjani, K.; Opeil, C.; Broido, D.; Chen, G.; Ren, Z. High Thermoelectric Performance by Resonant Dopant Indium in Nanostructured SnTe. Proc. Natl. Acad. Sci. U. S. A. 2013, 110, 13261−13266. (92) Tan, G.; Zhao, L.-D.; Shi, F.; Doak, J. W.; Lo, S.-H.; Sun, H.; Wolverton, C.; Dravid, V. P.; Uher, C.; Kanatzidis, M. G. High Thermoelectric Performance of P-Type SnTe Via a Synergistic Band Engineering and Nanostructuring Approach. J. Am. Chem. Soc. 2014, 136, 7006−7017. (93) Rogers, L. Valence Band Structure of SnTe. J. Phys. D: Appl. Phys. 1968, 1, 845−852. (94) Sitter, H.; Lischka, K.; Heinrich, H. Structure of the Second Valence Band in PbTe. Phys. Rev. B 1977, 16, 680−687. (95) Veis, A.; Kaidanov, V.; Kuteinikov, R.; Nemov, S.; Rudenko, S.; Ukhanov, Y. I. Investigation of Structure of Valence Band of Lead Selenide. Sov. Phys. Semicond. USSR 1978, 12, 161−163. (96) Zhao, L.-D.; Dravid, V. P.; Kanatzidis, M. G. The Panoscopic Approach to High Performance Thermoelectrics. Energy Environ. Sci. 2014, 7, 251−268. (97) Wu, H.; Chang, C.; Feng, D.; Xiao, Y.; Zhang, X.; Pei, Y.; Zheng, L.; Wu, D.; Gong, S.; Chen, Y.; He, J.; Kanatzidis, M. G.; Zhao, L.-D. Synergistically Optimized Electrical and Thermal Transport 12144

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

Properties of SnTe Via Alloying High-Solubility MnTe. Energy Environ. Sci. 2015, 8, 3298−3312. (98) He, J.; Tan, X.; Xu, J.; Liu, G.-Q.; Shao, H.; Fu, Y.; Wang, X.; Liu, Z.; Xu, J.; Jiang, H.; et al. Valence Band Engineering and Thermoelectric Performance Optimization in SnTe by Mn-Alloying Via a Zone-Melting Method. J. Mater. Chem. A 2015, 3, 19974−19979. (99) 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. (100) Takagiwa, Y.; Pei, Y.; Pomrehn, G.; Snyder, G. J. Dopants Effect on the Band Structure of PbTe Thermoelectric Material. Appl. Phys. Lett. 2012, 101, 092102. (101) Jaworski, C. M.; Kulbachinskii, V.; Heremans, J. P. Resonant Level Formed by Tin in Bi2Te3 and the Enhancement of RoomTemperature Thermoelectric Power. Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 80, 233201. (102) Zhang, Q.; Wang, H.; Liu, W.; Wang, H.; Yu, B.; Zhang, Q.; Tian, Z.; Ni, G.; Lee, S.; Esfarjani, K.; et al. Enhancement of Thermoelectric Figure-of-Merit by Resonant States of Aluminium Doping in Lead Selenide. Energy Environ. Sci. 2012, 5, 5246−5251. (103) Ohta, S.; Nomura, T.; Ohta, H.; Koumoto, K. HighTemperature Carrier Transport and Thermoelectric Properties of Heavily La-or Nb-Doped SrTiO3 Single Crystals. J. Appl. Phys. 2005, 97, 034106. (104) Okuda, T.; Nakanishi, K.; Miyasaka, S.; Tokura, Y. Large Thermoelectric Response of Metallic Perovskites: Sr1−xLaxTiO3 (0< x< 0.1). Phys. Rev. B: Condens. Matter Mater. Phys. 2001, 63, 113104. (105) Jalan, B.; Stemmer, S. Large Seebeck Coefficients and Thermoelectric Power Factor of La-Doped SrTiO3 Thin Films. Appl. Phys. Lett. 2010, 97, 042106. (106) Wunderlich, W.; Ohta, H.; Koumoto, K. Enhanced Effective Mass in Doped SrTiO3 and Related Perovskites. Phys. B 2009, 404, 2202−2212. (107) 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. (108) Kane, E. O. Band Structure of Indium Antimonide. J. Phys. Chem. Solids 1957, 1, 249−261. (109) Sofo, J.; Mahan, G. Electronic Structure of CoSb3: A NarrowBand-Gap Semiconductor. Phys. Rev. B: Condens. Matter Mater. Phys. 1998, 58, 15620−15623. (110) Korringa, J.; Gerritsen, A. The Cooperative Electron Phenomenon in Dilute Alloys. Physica 1953, 19, 457−507. (111) Friedel, J. On Some Electrical and Magnetic Properties of Metallic Solid Solutions. Can. J. Phys. 1956, 34, 1190−1211. (112) Hilsum, C. Simple Empirical Relationship between Mobility and Carrier Concentration. Electron. Lett. 1974, 10, 259−260. (113) Yu, B.; Zebarjadi, M.; Wang, H.; Lukas, K.; Wang, H.; Wang, D.; Opeil, C.; Dresselhaus, M.; Chen, G.; Ren, Z. Enhancement of Thermoelectric Properties by Modulation-Doping in Silicon Germanium Alloy Nanocomposites. Nano Lett. 2012, 12, 2077−2082. (114) Zebarjadi, M.; Joshi, G.; Zhu, G.; Yu, B.; Minnich, A.; Lan, Y.; Wang, X.; Dresselhaus, M.; Ren, Z.; Chen, G. Power Factor Enhancement by Modulation Doping in Bulk Nanocomposites. Nano Lett. 2011, 11, 2225−2230. (115) Zebarjadi, M.; Liao, B.; Esfarjani, K.; Dresselhaus, M.; Chen, G. Enhancing the Thermoelectric Power Factor by Using Invisible Dopants. Adv. Mater. 2013, 25, 1577−1582. (116) Pei, Y.-L.; Wu, H.; Wu, D.; Zheng, F.; He, J. High Thermoelectric Performance Realized in a Bicuseo System by Improving Carrier Mobility through 3D Modulation Doping. J. Am. Chem. Soc. 2014, 136, 13902−13908. (117) Valset, K.; Böttger, P.; Taftø, J.; Finstad, T. Thermoelectric Properties of Cu Doped ZnSb Containing Zn3P2 Particles. J. Appl. Phys. 2012, 111, 023703. (118) Koirala, M.; Zhao, H.; Pokharel, M.; Chen, S.; Dahal, T.; Opeil, C.; Chen, G.; Ren, Z. Thermoelectric Property Enhancement by Cu Nanoparticles in Nanostructured FeSb2. Appl. Phys. Lett. 2013, 102, 213111.

(119) Chen, C.-L.; Wang, H.; Chen, Y.-Y.; Day, T.; Snyder, G. J. Thermoelectric Properties of P-Type Polycrystalline SnSe Doped with Ag. J. Mater. Chem. A 2014, 2, 11171−11176. (120) Zhu, T.; Xu, Z.; He, J.; Shen, J.; Zhu, S.; Hu, L.; Tritt, T. M.; Zhao, X. Hot Deformation Induced Bulk Nanostructuring of Unidirectionally Grown P-Type (Bi,Sb)2Te3 Thermoelectric Materials. J. Mater. Chem. A 2013, 1, 11589−11594. (121) Hu, L.; Wu, H.; Zhu, T.; Fu, C.; He, J.; Ying, P.; Zhao, X. Tuning Multiscale Microstructures to Enhance Thermoelectric Performance of N-Type Bismuth-Telluride-Based Solid Solutions. Adv. Energy Mater. 2015, 5, 1500411. (122) Shen, J.; Hu, L.; Zhu, T.; Zhao, X. The Texture Related Anisotropy of Thermoelectric Properties in Bismuth Telluride Based Polycrystalline Alloys. Appl. Phys. Lett. 2011, 99, 124102. (123) Wang, H.; Sun, X.; Yan, X.; Huo, D.; Li, X.; Li, J.-G.; Ding, X. Fabrication and Thermoelectric Properties of Highly Textured Ca9Co12O28 Ceramic. J. Alloys Compd. 2014, 582, 294−298. (124) Kenfaui, D.; Chateigner, D.; Gomina, M.; Noudem, J. G. Texture, Mechanical and Thermoelectric Properties of Ca3Co4O9 Ceramics. J. Alloys Compd. 2010, 490, 472−479. (125) Qiu, B.; Bao, H.; Zhang, G.; Wu, Y.; Ruan, X. Molecular Dynamics Simulations of Lattice Thermal Conductivity and Spectral Phonon Mean Free Path of Pbte: Bulk and Nanostructures. Comput. Mater. Sci. 2012, 53, 278−285. (126) Klemens, P. Thermal Resistance Due to Point Defects at High Temperatures. Phys. Rev. 1960, 119, 507−509. (127) Callaway, J. Model for Lattice Thermal Conductivity at Low Temperatures. Phys. Rev. 1959, 113, 1046−1051. (128) Callaway, J.; von Baeyer, H. C. Effect of Point Imperfections on Lattice Thermal Conductivity. Phys. Rev. 1960, 120, 1149−1154. (129) Abeles, B. Lattice Thermal Conductivity of Disordered Semiconductor Alloys at High Temperatures. Phys. Rev. 1963, 131, 1906−1911. (130) Tan, G.; Chi, H.; Liu, W.; Zheng, Y.; Tang, X.; He, J.; Uher, C. Toward High Thermoelectric Performance P-Type FeSb2.2Te0.8 Via in Situ Formation of Insb Nanoinclusions. J. Mater. Chem. C 2015, 3, 8372−8380. (131) Tan, G.; Liu, W.; Chi, H.; Su, X.; Wang, S.; Yan, Y.; Tang, X.; Wong-Ng, W.; Uher, C. Realization of High Thermoelectric Performance in P-Type Unfilled Ternary Skutterudites FeSb2+xTe1−x via Band Structure Modification and Significant Point Defect Scattering. Acta Mater. 2013, 61, 7693−7704. (132) Tan, G.; Wang, S.; Tang, X. High Thermoelectric Figure of Merit of P-Type Ternary Unfilled Skutterudite FeSb2Te Via Ge Doping. Sci. Adv. Mater. 2013, 5, 1974−1982. (133) Shi, X.; Yang, J.; Bai, S.; Yang, J.; Wang, H.; Chi, M.; Salvador, J. R.; Zhang, W.; Chen, L. Wong - Ng, W. On the Design of High Efficiency Thermoelectric Clathrates through a Systematic Cross Substitution of Framework Elements. Adv. Funct. Mater. 2010, 20, 755−763. (134) Meisner, G.; Morelli, D.; Hu, S.; Yang, J.; Uher, C. Structure and Lattice Thermal Conductivity of Fractionally Filled Skutterudites: Solid Solutions of Fully Filled and Unfilled End Members. Phys. Rev. Lett. 1998, 80, 3551−3554. (135) Nolas, G.; Cohn, J.; Slack, G.; Schujman, S. Semiconducting Ge Clathrates: Promising Candidates for Thermoelectric Applications. Appl. Phys. Lett. 1998, 73, 178−180. (136) May, A. F.; Singh, D. J.; Snyder, G. J. Influence of Band Structure on the Large Thermoelectric Performance of Lanthanum Telluride. Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 79, 153101. (137) Delaire, O.; May, A. F.; McGuire, M. A.; Porter, W. D.; Lucas, M. S.; Stone, M. B.; Abernathy, D. L.; Ravi, V.; Firdosy, S.; Snyder, G. Phonon Density of States and Heat Capacity of La3‑xTe4. Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 80, 184302. (138) Pei, Y.; Morelli, D. T. Vacancy Phonon Scattering in Thermoelectric In2Te3−InSb Solid Solutions. Appl. Phys. Lett. 2009, 94, 122112. 12145

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

Skutterudites CeFe4−xZnxSb12. J. Solid State Chem. 2012, 187, 316− 322. (156) Tan, G.; Wang, S.; Tang, X.; Li, H.; Uher, C. Preparation and Thermoelectric Properties of Ga-Substituted P-Type Fully Filled Skutterudites CeFe4‑xGaxSb12. J. Solid State Chem. 2012, 196, 203−208. (157) He, J.; Sootsman, J. R.; Girard, S. N.; Zheng, J.-C.; Wen, J.; Zhu, Y.; Kanatzidis, M. G.; Dravid, V. P. On the Origin of Increased Phonon Scattering in Nanostructured PbTe Based Thermoelectric Materials. J. Am. Chem. Soc. 2010, 132, 8669−8675. (158) Fauchais, P.; Montavon, G.; Lima, R.; Marple, B. Engineering a New Class of Thermal Spray Nano-Based Microstructures from Agglomerated Nanostructured Particles, Suspensions and Solutions: An Invited Review. J. Phys. D: Appl. Phys. 2011, 44, 093001. (159) Rowe, D.; Shukla, V.; Savvides, N. Phonon Scattering at Grain Boundaries in Heavily Doped Fine-Grained Silicon-Germanium Alloys. Nature 1981, 290, 765−766. (160) Zhao, Y.; Dyck, J. S.; Burda, C. Toward High-Performance Nanostructured Thermoelectric Materials: The Progress of Bottom-up Solution Chemistry Approaches. J. Mater. Chem. 2011, 21, 17049− 17058. (161) Tan, G.; Zheng, Y.; Tang, X. High Thermoelectric Performance of Nonequilibrium Synthesized CeFe4Sb12 Composite with Multi-Scaled Nanostructures. Appl. Phys. Lett. 2013, 103, 183904. (162) Su, X.; Fu, F.; Yan, Y.; Zheng, G.; Liang, T.; Zhang, Q.; Cheng, X.; Yang, D.; Chi, H.; Tang, X.; et al. Self-Propagating HighTemperature Synthesis for Compound Thermoelectrics and New Criterion for Combustion Processing. Nat. Commun. 2014, 5, 4908. (163) Liang, T.; Su, X.; Yan, Y.; Zheng, G.; Zhang, Q.; Chi, H.; Tang, X.; Uher, C. Ultra-Fast Synthesis and Thermoelectric Properties of Te Doped Skutterudites. J. Mater. Chem. A 2014, 2, 17914−17918. (164) 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. (165) Wang, W.-Z.; Zeng, B.-Q.; Yang, J.; Poudel, B.; Huang, J.; Naughton, M. J.; Ren, Z. Aligned Ultralong Zno Nanobelts and Their Enhanced Field Emission. Adv. Mater. 2006, 18, 3275−3278. (166) Martín-González, M.; Snyder, G. J.; Prieto, A. L.; Gronsky, R.; Sands, T.; Stacy, A. M. Direct Electrodeposition of Highly Dense 50 nm Bi2Te3‑ySey Nanowire Arrays. Nano Lett. 2003, 3, 973−977. (167) Ge, J.-P.; Li, Y.-D. Ultrasonic Synthesis of Nanocrystals of Metal Selenides and Tellurides. J. Mater. Chem. 2003, 13, 911−915. (168) Mi, J.; Zhao, X.; Zhu, T.; Tu, J.; Cao, G. Solvothermal Synthesis and Electrical Transport Properties of Skutterudite CoSb3. J. Alloys Compd. 2006, 417, 269−272. (169) Schlecht, S.; Erk, C.; Yosef, M. Nanoscale Zinc Antimonides: Synthesis and Phase Stability. Inorg. Chem. 2006, 45, 1693−1697. (170) Tian, Z.; Garg, J.; Esfarjani, K.; Shiga, T.; Shiomi, J.; Chen, G. Phonon Conduction in PbSe, PbTe, and PbTe1‑xSex from FirstPrinciples Calculations. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85, 184303. (171) Esfarjani, K.; Chen, G.; Stokes, H. T. Heat Transport in Silicon from First-Principles Calculations. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 84, 085204. (172) Zheng, Y.; Zhang, Q.; Su, X.; Xie, H.; Shu, S.; Chen, T.; Tan, G.; Yan, Y.; Tang, X.; Uher, C.; Snyder, G. J. Mechanically Robust Bisbte Alloys with Superior Thermoelectric Performance: A Case Study of Stable Hierarchical Nanostructured Thermoelectric Materials. Adv. Energy Mater. 2015, 5, 1401391. (173) Guin, S. N.; Negi, D. S.; Datta, R.; Biswas, K. Nanostructuring, Carrier Engineering and Bond Anharmonicity Synergistically Boost the Thermoelectric Performance of P-Type AgSbSe2−ZnSe. J. Mater. Chem. A 2014, 2, 4324−4331. (174) Bhattacharya, S.; Bohra, A.; Basu, R.; Bhatt, R.; Ahmad, S.; Meshram, K.; Debnath, A.; Singh, A.; Sarkar, S. K.; Navneethan, M. High Thermoelectric Performance of (AgCrSe2)0.5(CuCrSe2)0.5 NanoComposites Having All-Scale Natural Hierarchical Architectures. J. Mater. Chem. A 2014, 2, 17122−17129.

(139) Sun, H.; Lu, X.; Chi, H.; Morelli, D. T.; Uher, C. Highly Efficient (In2Te3)x(GeTe)3−3x Thermoelectric Materials: A Substitute for Tags. Phys. Chem. Chem. Phys. 2014, 16, 15570−15575. (140) Tan, G.; Zeier, W. G.; Shi, f.; Wang, P.; Snyder, G. J.; Dravid, V. P.; Kanatzidis, M. G. High Thermoelectric Performance SnTe− In2Te3 Solid Solutions Enabled by Resonant Levels and Strong Vacancy Phonon Scattering. Chem. Mater. 2015, 27, 7801−7811. (141) Tan, G.; Shi, F.; Sun, H.; Zhao, L.-D.; Uher, C.; Dravid, V. P.; Kanatzidis, M. G. SnTe−AgBiTe2 as an Efficient Thermoelectric Material with Low Thermal Conductivity. J. Mater. Chem. A 2014, 2, 20849−20854. (142) Yang, F.; Dames, C. Mean Free Path Spectra as a Tool to Understand Thermal Conductivity in Bulk and Nanostructures. Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, 035437. (143) Li, J.; Tan, Q.; Li, J. F.; Liu, D. W.; Li, F.; Li, Z. Y.; Zou, M.; Wang, K. BiSbTe - Based Nanocomposites with High ZT: The Effect of Sic Nanodispersion on Thermoelectric Properties. Adv. Funct. Mater. 2013, 23, 4317−4323. (144) Min, Y.; Roh, J. W.; Yang, H.; Park, M.; Kim, S. I.; Hwang, S.; Lee, S. M.; Lee, K. H.; Jeong, U. Surfactant - Free Scalable Synthesis of Bi2Te3 and Bi2Se3 Nanoflakes and Enhanced Thermoelectric Properties of Their Nanocomposites. Adv. Mater. 2013, 25, 1425−1429. (145) Shen, J.-J.; Zhu, T.-J.; Zhao, X.-B.; Zhang, S.-N.; Yang, S.-H.; Yin, Z.-Z. Recrystallization Induced in Situ Nanostructures in Bulk Bismuth Antimony Tellurides: A Simple Top Down Route and Improved Thermoelectric Properties. Energy Environ. Sci. 2010, 3, 1519−1523. (146) Hu, L.; Gao, H.; Liu, X.; Xie, H.; Shen, J.; Zhu, T.; Zhao, X. Enhancement in Thermoelectric Performance of Bismuth Telluride Based Alloys by Multi-Scale Microstructural Effects. J. Mater. Chem. 2012, 22, 16484−16490. (147) Androulakis, J.; Lin, C.-H.; Kong, H.-J.; Uher, C.; Wu, C.-I.; Hogan, T.; Cook, B. A.; Caillat, T.; Paraskevopoulos, K. M.; Kanatzidis, M. G. Spinodal Decomposition and Nucleation and Growth as a Means to Bulk Nanostructured Thermoelectrics: Enhanced Performance in Pb1‑xSnxTe-PbS. J. Am. Chem. Soc. 2007, 129, 9780−9788. (148) Biswas, K.; He, J.; Wang, G.; Lo, S.-H.; Uher, C.; Dravid, V. P.; Kanatzidis, M. G. High Thermoelectric Figure of Merit in Nanostructured P-Type PbTe−MTe (M= Ca, Ba). Energy Environ. Sci. 2011, 4, 4675−4684. (149) Ahn, K.; Biswas, K.; He, J.; Dravid, V.; Kanatzidis, M. G. Enhanced Thermoelectric Properties of P-Type Nanostructured PbTe−MTe (M= Cd, Hg) Materials. Energy Environ. Sci. 2013, 6, 1529−1537. (150) Lee, Y.; Lo, S.-H.; Androulakis, J.; Wu, C.-I.; Zhao, L.-D.; Chung, D.-Y.; Hogan, T. P.; Dravid, V. P.; Kanatzidis, M. G. HighPerformance Tellurium-Free Thermoelectrics: All-Scale Hierarchical Structuring of P-Type PbSe−MSe Systems (M= Ca, Sr, Ba). J. Am. Chem. Soc. 2013, 135, 5152−5160. (151) 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. (152) Zhao, L.-D.; Lo, S.-H.; He, J.; Li, H.; Biswas, K.; Androulakis, J.; Wu, C.-I.; Hogan, T. P.; Chung, D.-Y.; Dravid, V. P.; et al. High Performance Thermoelectrics from Earth-Abundant Materials: Enhanced Figure of Merit in PbS by Second Phase Nanostructures. J. Am. Chem. Soc. 2011, 133, 20476−20487. (153) Li, H.; Tang, X.; Zhang, Q.; Uher, C. High Performance InxCeyCo4Sb12 Thermoelectric Materials with in Situ Forming Nanostructured Insb Phase. Appl. Phys. Lett. 2009, 94, 102114. (154) Xiong, Z.; Chen, X.; Huang, X.; Bai, S.; Chen, L. High Thermoelectric Performance of Yb0.26Co4Sb12/yGaSb Nanocomposites Originating from Scattering Electrons of Low Energy. Acta Mater. 2010, 58, 3995−4002. (155) Tan, G.; Wang, S.; Li, H.; Yan, Y.; Tang, X. Enhanced Thermoelectric Performance in Zinc Substituted P-Type Filled 12146

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

(175) Zhang, Q.; Ai, X.; Wang, L.; Chang, Y.; Luo, W.; Jiang, W.; Chen, L. Improved Thermoelectric Performance of Silver Nanoparticles - Dispersed Bi2Te3 Composites Deriving from Hierarchical Two - Phased Heterostructure. Adv. Funct. Mater. 2015, 25, 966−976. (176) Lan, J.-L.; Liu, Y.; Lin, Y.-H.; Nan, C.-W.; Cai, Q.; Yang, X. Enhanced Thermoelectric Performance of In2O3-Based Ceramics via Nanostructuring and Point Defect Engineering. Sci. Rep. 2015, 5, 7783. (177) Wu, H.; Zheng, F.; Wu, D.; Ge, Z.-H.; Liu, X.; He, J. Advanced Electron Microscopy for Thermoelectric Materials. Nano Energy 2015, 13, 626−650. (178) Zhou, M.; Li, J.-F.; Kita, T. Nanostructured AgPbmSbTem+2 System Bulk Materials with Enhanced Thermoelectric Performance. J. Am. Chem. Soc. 2008, 130, 4527−4532. (179) Quarez, E.; Hsu, K.-F.; Pcionek, R.; Frangis, N.; Polychroniadis, E.; Kanatzidis, M. G. Nanostructuring, Compositional Fluctuations, and Atomic Ordering in the Thermoelectric Materials AgPbmSbTe2+m: The Myth of Solid Solutions. J. Am. Chem. Soc. 2005, 127, 9177−9190. (180) Hogan, T. P.; Downey, A.; Short, J.; D’Angelo, J.; Wu, C.-I.; Quarez, E.; Androulakis, J.; Poudeu, P. F.; Sootsman, J. R.; Chung, D.Y.; et al. Nanostructured Thermoelectric Materials and High-Efficiency Power-Generation Modules. J. Electron. Mater. 2007, 36, 704−710. (181) Cook, B. A.; Kramer, M. J.; Harringa, J. L.; Han, M. K.; Chung, D. Y.; Kanatzidis, M. G. Analysis of Nanostructuring in High Figure of - Merit Ag1−xPbmSbTe2+m Thermoelectric Materials. Adv. Funct. Mater. 2009, 19, 1254−1259. (182) Tit, N.; Peressi, M.; Baroni, S. Ab Initio Calculation of the Band Offset at Strained GaAs/InAs (001) Heterojunctions. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 48, 17607−17610. (183) Zhao, L. D.; Dravid, V. P.; Kanatzidis, M. G. The Panoscopic Approach to High Performance Thermoelectrics. Energy Environ. Sci. 2014, 7, 251−268. (184) Saha, S. K. Exploring the Origin of Ultralow Thermal Conductivity in Layered BiOCuSe. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, 41202. (185) Ding, J. X.; Xu, B.; Lin, Y. H.; Nan, C. W.; Liu, W. Lattice Vibration Modes of the Layered Material BiCuSeO and First Principles Study of Its Thermoelectric Properties. New J. Phys. 2015, 17, 083012. (186) Dennler, G.; Chmielowski, R.; Jacob, S.; Capet, F.; Roussel, P.; Zastrow, S.; Nielsch, K.; Opahle, I.; Madsen, G. K. H. Are Binary Copper Sulfides/Selenides Really New and Promising Thermoelectric Materials? Adv. Energy Mater. 2014, 4, 1301581. (187) Brown, S. R.; Kauzlarich, S. M.; Gascoin, F.; Snyder, G. J. Yb14MnSb11: New High Efficiency Thermoelectric Material for Power Generation. Chem. Mater. 2006, 18, 1873−1877. (188) Kurosaki, K.; Kosuga, A.; Muta, H.; Uno, M.; Yamanaka, S. Ag9TlTe5: A High-Performance Thermoelectric Bulk Material with Extremely Low Thermal Conductivity. Appl. Phys. Lett. 2005, 87, 061919. (189) Nielsen, M. D.; Ozolins, V.; Heremans, J. P. Lone Pair Electrons Minimize Lattice Thermal Conductivity. Energy Environ. Sci. 2013, 6, 570−578. (190) Skoug, E. J.; Morelli, D. T. Role of Lone-Pair Electrons in Producing Minimum Thermal Conductivity in Nitrogen-Group Chalcogenide Compounds. Phys. Rev. Lett. 2011, 107, 235901. (191) Morelli, D. T.; Jovovic, V.; Heremans, J. P. Intrinsically Minimal Thermal Conductivity in Cubic I-V-VI2 Semiconductors. Phys. Rev. Lett. 2008, 101, 035901. (192) Zhang, Q.; Chere, E. K.; Sun, J. Y.; Cao, F.; Dahal, K.; Chen, S.; Chen, G.; Ren, Z. F. Studies on Thermoelectric Properties of NType Polycrystalline SnSe1‑xSx by Iodine Doping. Adv. Energy Mater. 2015, 5, 1500360. (193) Goldsmid, H. Conversion Efficiency and Figure-of-Merit; CRC Press, 1995. (194) Bell, L. E. Cooling, Heating, Generating Power, and Recovering Waste Heat with Thermoelectric Systems. Science 2008, 321, 1457−1461.

(195) Sassi, S.; Candolfi, C.; Vaney, J. B.; Ohorodniichuk, V.; Masschelein, P.; Dauscher, A.; Lenoir, B. Assessment of the Thermoelectric Performance of Polycrystalline P-Type SnSe. Appl. Phys. Lett. 2014, 104, 212105. (196) Sassi, S.; Candolfi, C.; Vaney, J. B.; Ohorodniichuk, V.; Masschelein, P.; Dauscher, A.; Lenoir, B. Transport Properties of Polycrystalline P-Type SnSe. Mater. Today Proc. 2015, 2, 690−698. (197) Turkes, P.; Pluntke, C.; Helbig, R. Thermal Conductivity of SnO2 Single Crystals. J. Phys. C: Solid State Phys. 1980, 13, 4941−4951. (198) He, B.; Wiendlocha, B.; Prakash, A.; Heremans, J. P. 2014 MRS Fall Meeting & Exhibit, Boston, MA, 2014; CC7.01. (199) Serrano-Sanchez, F.; Gharsallah, M.; Nemes, N. M.; Mompean, F. J.; Martinez, J. L.; Alonso, J. A. Record Seebeck Coefficient and Extremely Low Thermal Conductivity in Nanostructured SnSe. Appl. Phys. Lett. 2015, 106, 083902. (200) Tan, Q.; Li, J. F. Thermoelectric Properties of Sn-S Bulk Materials Prepared by Mechanical Alloying and Spark Plasma Sintering. J. Electron. Mater. 2014, 43, 2435−2439. (201) Tan, Q.; Zhao, L. D.; Li, J. F.; Wu, C. F.; Wei, T. R.; Xing, Z. B.; Kanatzidis, M. G. Thermoelectrics with Earth Abundant Elements: Low Thermal Conductivity and High Thermopower in Doped SnS. J. Mater. Chem. A 2014, 2, 17302−17306. (202) Parker, D.; Singh, D. J. First Principles Investigations of the Thermoelectric Behavior of Tin Sulfide. J. Appl. Phys. 2010, 108, 1223−1239. (203) Bera, C.; Jacob, S.; Opahle, I.; Gunda, N. S. H.; Chmielowski, R.; Dennler, G.; Madsen, G. K. H. Integrated Computational Materials Discovery of Silver Doped Tin Sulfide as a Thermoelectric Material. Phys. Chem. Chem. Phys. 2014, 16, 19894−19899. (204) Zhao, L. D.; Tan, G. J.; Hao, S. Q.; He, J. Q.; Pei, Y. L.; Chi, H.; Wang, H.; Gong, S. K.; Xu, H. B.; Dravid, V. P.; et al. Ultrahigh Power Factor and Thermoelectric Performance in Hole-Doped SingleCrystal SnSe. Science 2016, 351, 141−144. (205) Peng, K.; Lu, X.; Zhan, H.; Hui, S.; Tang, X.; Wang, G.; Dai, J.; Uher, C.; Wang, G.; Zhou, X. Broad Temperature Plateau for High ZTs in Heavily Doped P-Type SnSe Single Crystals. Energy Environ. Sci. 2016, 9, 454−460. (206) Pei, Y. Z.; Wang, H.; Snyder, G. J. Band Engineering of Thermoelectric Materials. Adv. Mater. 2012, 24, 6125−6135. (207) Zhao, L. D.; Hao, S. Q.; Lo, S. H.; Wu, C. I.; Zhou, X. Y.; 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. (208) Poudeu, P. F. R.; D’Angelo, J.; Downey, A. D.; Short, J. L.; Hogan, T. P.; Kanatzidis, M. G. High Thermoelectric Figure of Merit and Nanostructuring in Bulk P-Type Na1‑xPbmSbyTem+2. Angew. Chem., Int. Ed. 2006, 45, 3835−3839. (209) Zhao, L. D.; Berardan, D.; Pei, Y. L.; Byl, C.; Pinsard-Gaudart, L.; Dragoe, N. Bi1‑xSrxCuSeO Oxyselenides as Promising Thermoelectric Materials. Appl. Phys. Lett. 2010, 97, 092118. (210) Barreteau, C.; Pan, L.; Amzallag, E.; Zhao, L. D.; Berardan, D.; Dragoe, N. Layered Oxychalcogenide in the Bi-Cu-O-Se System as Good Thermoelectric Materials. Semicond. Sci. Technol. 2014, 29, 64001−64010. (211) Barreteau, C.; Pan, L.; Pei, Y. L.; Zhao, L. D.; Berardan, D.; Dragoe, N. Oxychalcogenides as New Efficient P-Type Thermoelectric Materials. Funct. Mater. Lett. 2013, 6, 1340007. (212) Li, J.; Sui, J. H.; Pei, Y. L.; Barreteau, C.; Berardan, D.; Dragoe, N.; Cai, W.; He, J. Q.; Zhao, L. D. A High Thermoelectric Figure of Merit ZT > 1 in Ba Heavily Doped BiCuSeO Oxyselenides. Energy Environ. Sci. 2012, 5, 8543−8547. (213) Sui, J. H.; Li, J.; He, J. Q.; Pei, Y. L.; Berardan, D.; Wu, H. J.; Dragoe, N.; Cai, W.; Zhao, L. D. Texturation Boosts the Thermoelectric Performance of Bicuseo Oxyselenides. Energy Environ. Sci. 2013, 6, 2916−2920. (214) Fu, C. G.; Zhu, T. J.; Liu, Y. T.; Xie, H. H.; Zhao, X. B. Band Engineering of High Performance P-Type FeNbSb Based Half-Heusler 12147

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

Thermoelectric Materials for Figure of Merit ZT > 1. Energy Environ. Sci. 2015, 8, 216−220. (215) Lan, J. L.; Liu, Y. C.; Zhan, B.; Lin, Y. H.; Zhang, B. P.; Yuan, X.; Zhang, W. Q.; Xu, W.; Nan, C. W. Enhanced Thermoelectric Properties of Pb-Doped BiCuSeO Ceramics. Adv. Mater. 2013, 25, 5086−5090. (216) Li, F.; Li, J. F.; Zhao, L. D.; Xiang, K.; Liu, Y.; Zhang, B. P.; Lin, Y. H.; Nan, C. W.; Zhu, H. M. Polycrystalline BiCuSeO Oxide as a Potential Thermoelectric Material. Energy Environ. Sci. 2012, 5, 7188− 7195. (217) Liu, Y.; Zhao, L. D.; Liu, Y. C.; Lan, J. L.; Xu, W.; Li, F.; Zhang, B. P.; Berardan, D.; Dragoe, N.; Lin, Y. H.; et al. Remarkable Enhancement in Thermoelectric Performance of BiCuSeO by Cu Deficiencies. J. Am. Chem. Soc. 2011, 133, 20112−20115. (218) Li, Z.; Xiao, C.; Fan, S. J.; Deng, Y.; Zhang, W. S.; Ye, B. J.; Xie, Y. Dual Vacancies: An Effective Strategy Realizing Synergistic Optimization of Thermoelectric Property in BiCuSeO. J. Am. Chem. Soc. 2015, 137, 6587−6593. (219) Pei, Y. L.; He, J. Q.; Li, J. F.; Li, F.; Liu, Q. J.; Pan, W.; Barreteau, C.; Berardan, D.; Dragoe, N.; Zhao, L. D. High Thermoelectric Performance of Oxyselenides: Intrinsically Low Thermal Conductivity of Ca-Doped BiCuSeO. NPG Asia Mater. 2013, 5, e47. (220) Nolas, G. S.; Poon, J.; Kanatzidis, M. Recent Developments in Bulk Thermoelectric Materials. MRS Bull. 2006, 31, 199−205. (221) Kirkham, M. J.; dos Santos, A. M.; Rawn, C. J.; Lara-Curzio, E.; Sharp, J. W.; Thompson, A. J. Ab Initio Determination of Crystal Structures of the Thermoelectric Material MgAgSb. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85, 144120. (222) Ying, P.; Liu, X.; Fu, C.; Yue, X.; Xie, H.; Zhao, X.; Zhang, W.; Zhu, T. High Performance Alpha-MgAgSb Thermoelectric Materials for Low Temperature Power Generation. Chem. Mater. 2015, 27, 909−913. (223) Shuai, J.; Kim, H. S.; Lan, Y.; Chen, S.; Liu, Y.; Zhao, H.; Sui, J.; Ren, Z. Study on Thermoelectric Performance by Na Doping in Nanostructured Mg1‑xNaxAg0.97Sb0.99. Nano Energy 2015, 11, 640−646. (224) Kraemer, D.; Sui, J.; McEnaney, K.; Zhao, H.; Jie, Q.; Ren, Z. F.; Chen, G. High Thermoelectric Conversion Efficiency of MgAgSbBased Material with Hot-Pressed Contacts. Energy Environ. Sci. 2015, 8, 1299−1308. (225) Becquerel, A. C. De Quelques Phénomènes Electriques Produits Par La Pression et le Clivage des Métaux. Ann. Chim. Phys. 1827, 36, 265−271. (226) Becquerel, A. C. Memoire Sur L’eIectro-chimie et L’emploi de L’electricite Pour Operer des Combinaisons. Ann. Chim. Phys. 1829, 41, 23−26. (227) Purohit, R.; Sharma, B.; Sreedhar, A. Diffusion of Copper in CdS Crystal from Cu2S Layer. J. Appl. Phys. 1969, 40, 4677−4678. (228) Ge, Z. H.; Zhang, B. P.; Chen, Y. X.; Yu, Z. X.; Liu, Y.; Li, J. F. Synthesis and Transport Property of Cu1.8S as a Promising Thermoelectric Compound. Chem. Commun. 2011, 47, 12697−12699. (229) Yu, B.; Liu, W.; Chen, S.; Wang, H.; Wang, H.; Chen, G.; Ren, Z. Thermoelectric Properties of Copper Selenide with Ordered Selenium Layer and Disordered Copper Layer. Nano Energy 2012, 1, 472−478. (230) Yamamoto, K.; Kashida, S. X-Ray Study of the Average Structures of Cu2Se and Cu1.8S in the Room Temperature and the High Temperature Phases. J. Solid State Chem. 1991, 93, 202−211. (231) Oliveria, M.; McMullan, R.; Wuensch, B. Single Crystal Neutron Diffraction Analysis of the Cation Distribution in the HighTemperature Phases A-Cu2‑xS, α-Cu2‑xSe, and α-Ag2Se. Solid State Ionics 1988, 28-30, 1332−1337. (232) Zhong, B.; Zhang, Y.; Li, W.; Chen, Z.; Cui, J.; Li, W.; Xie, Y.; Hao, Q.; He, Q. High Superionic Conduction Arising from Aligned Large Lamellae and Large Figure of Merit in Bulk Cu1.94Al0.02Se. Appl. Phys. Lett. 2014, 105, 123902. (233) Cai, S.; Liu, Z.; Sun, J.; Li, R.; Fei, W.; Sui, J. Enhancement of Thermoelectric Properties by Na Doping in Te-Free P-Type AgSbSe2. Dalton Trans. 2015, 44, 1046−1051.

(234) Armand, M.; Endres, F.; MacFarlane, D. R.; Ohno, H.; Scrosati, B. Ionic-Liquid Materials for the Electrochemical Challenges of the Future. Nat. Mater. 2009, 8, 621−629. (235) Brown, D. R.; Day, T.; Caillat, T.; Snyder, G. J. Chemical Stability of (Ag,Cu)2Se: A Historical Overview. J. Electron. Mater. 2013, 42, 2014−2019. (236) Chung, D. Y.; Hogan, T. P.; Rocci-Lane, M.; Brazis, P.; Ireland, J. R.; Kannewurf, C. R.; Bastea, M.; Uher, C.; Kanatzidis, M. G. A New Thermoelectric Material: CsBi4Te6. J. Am. Chem. Soc. 2004, 126, 6414−6428. (237) Hsu, K. F.; Chung, D. Y.; Lal, S.; Mrotzek, A.; Kyratsi, T.; Hogan, T.; Kanatzidis, M. G. CsMBi3Te6 and CsM2Bi3Te7 (M = Pb,Sn): New Thermoelectric Compounds with Low-Dimensional Structures. J. Am. Chem. Soc. 2002, 124, 2410−2411. (238) Bilc, D. I.; Mahanti, S. D.; Kyratsi, T.; Chung, D. Y.; Kanatzidis, M. G.; Larson, P. Electronic Structure of K2Bi8Se13. Phys. Rev. B: Condens. Matter Mater. Phys. 2005, 71, 085116. (239) Hoang, K.; Tomic, A.; Mahanti, S. D.; Kyratsi, T.; Chung, D.Y.; Tessmer, S. H.; Kanatzidis, M. G. Role of K/Bi Disorder in the Electronic Structure of Beta-K2Bi8Se13. Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 80, 125112. (240) Chung, D.-Y.; Choi, K.-S.; Iordanidis, L.; Schindler, J. L.; Brazis, P. W.; Kannewurf, C. R.; Chen, B.; Hu, S.; Uher, C.; Kanatzidis, M. G. High Thermopower and Low Thermal Conductivity in Semiconducting Ternary K−Bi−Se Compounds. Synthesis and Properties of β-K2Bi8Se13 and K2.5Bi8.5Se14 and Their Sb Analogues. Chem. Mater. 1997, 9, 3060−3071. (241) Vedernikov, M. V.; Kutasov, V. A.; Luk’yanova, L. N.; Konstantinov, P. P. Sixteenth International Conference on Thermoelectrics, Dresden, 1997; pp 56−62. (242) Martin-Lopez, R.; Dauscher, A.; Devaux, X.; Lenoir, B.; Scherrer, H.; Zandona, M. Sixteenth International Conference on Thermoelectrics, Dresden, 1997; pp 184−187. (243) Clarke, D. R. Materials Selection Guidelines for Low Thermal Conductivity Thermal Barrier Coatings. Surf. Coat. Technol. 2003, 163164, 67−74. (244) Zhang, Y.; Ozolins, V.; Morelli, D.; Wolverton, C. Prediction of New Stable Compounds and Promising Thermoelectrics in the Cu-SbSe System. Chem. Mater. 2014, 26, 3427−3435. (245) Zhang, Y.; Skoug, E.; Cain, J.; Ozolins, V.; Morelli, D.; Wolverton, C. First-Principles Description of Anomalously Low Lattice Thermal Conductivity in Thermoelectric Cu-Sb-Se Ternary Semiconductors. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85, 054306. (246) Waghmare, U. V.; Spaldin, N. A.; Kandpal, H. C.; Seshadri, R. First-Principles Indicators of Metallicity and Cation Off-Centricity in the Iv-Vi Rocksalt Chalcogenides of Divalent Ge, Sn, and Pb. Phys. Rev. B: Condens. Matter Mater. Phys. 2003, 67, 125111. (247) Raulot, J. M.; Baldinozzi, G.; Seshadri, R.; Cortona, P. An AbInitio Study of the Role of Lone Pairs in the Structure and InsulatorMetal Transition in SnO and PbO. Solid State Sci. 2002, 4, 467−474. (248) Stoltzfus, M. W.; Woodward, P. M.; Seshadri, R.; Klepeis, J.-H.; Bursten, B. Structure and Bonding in SnWO4, PbWO4, and BiVO4: Lone Pairs vs Inert Pairs. Inorg. Chem. 2007, 46, 3839−3850. (249) Wei, T.-R.; Li, F.; Li, J.-F. Enhanced Thermoelectric Performance of Nonstoichiometric Compounds Cu3‑xSbSe4 by Cu Deficiencies. J. Electron. Mater. 2014, 43, 2229−2238. (250) Wei, T.-R.; Wang, H.; Gibbs, Z. M.; Wu, C.-F.; Snyder, G. J.; Li, J.-F. Thermoelectric Properties of Sn-Doped P-Type Cu3SbSe4: A Compound with Large Effective Mass and Small Band Gap. J. Mater. Chem. A 2014, 2, 13527−13533. (251) Lai, W.; Wang, Y.; Morelli, D. T.; Lu, X. From Bonding Asymmetry to Anharmonic Rattling in Cu12Sb4S13 tetrahedrites: When Lone-Pair Electrons Are Not So Lonely. Adv. Funct. Mater. 2015, 25, 3648−3657. (252) Noda, Y.; Nishida, I. A.; Kang, Y. S.; Niino, M.; Ieee, I. Seventeenth International Conference on Thermoelectrics, Nagoya, 1998; pp 350−353. 12148

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149

Chemical Reviews

Review

(253) Wang, H.; Li, J.-F.; Zou, M.; Sui, T. Synthesis and Transport Property of AgSbTe2 as a Promising Thermoelectric Compound. Appl. Phys. Lett. 2008, 93, 202106. (254) Guin, S. N.; Biswas, K. Sb Deficiencies Control Hole Transport and Boost the Thermoelectric Performance of P-Type AgSbSe2. J. Mater. Chem. C 2015, 3, 10415−10421. (255) Guin, S. N.; Chatterjee, A.; Negi, D. S.; Datta, R.; Biswas, K. High Thermoelectric Performance in Tellurium Free P-Type AgSbSe2. Energy Environ. Sci. 2013, 6, 2603−2608. (256) Pan, L.; Berardan, D.; Dragoe, N. High Thermoelectric Properties of N-Type AgBiSe2. J. Am. Chem. Soc. 2013, 135, 4914− 4917. (257) Guin, S. N.; Srihari, V.; Biswas, K. Promising Thermoelectric Performance in N-Type AgBiSe2: Effect of Aliovalent Anion Doping. J. Mater. Chem. A 2015, 3, 648−655. (258) Lu, X.; Morelli, D. T.; Xia, Y.; Zhou, F.; Ozolins, V.; Chi, H.; Zhou, X.; Uher, C. High Performance Thermoelectricity in EarthAbundant Compounds Based on Natural Mineral Tetrahedrites. Adv. Energy Mater. 2013, 3, 342−348. (259) Suekuni, K.; Tsuruta, K.; Kunii, M.; Nishiate, H.; Nishibori, E.; Maki, S.; Ohta, M.; Yamamoto, A.; Koyano, M. High-Performance Thermoelectric Mineral Cu12‑xNixSb4S13 Tetrahedrite. J. Appl. Phys. 2013, 113, 043712. (260) Suekuni, K.; Tsuruta, K.; Ariga, T.; Koyano, M. Thermoelectric Properties of Mineral Tetrahedrites Cu10Tr2Sb4S13 with Low Thermal Conductivity. Appl. Phys. Express 2012, 5, 051201. (261) Chetty, R.; Bali, A.; Mallik, R. C. Tetrahedrites as Thermoelectric Materials: An Overview. J. Mater. Chem. C 2015, 3, 12364−12378. (262) Lu, X.; Morelli, D. The Effect of Te Substitution for Sb on Thermoelectric Properties of Tetrahedrite. J. Electron. Mater. 2014, 43, 1983−1987. (263) Lu, X.; Morelli, D. T. Natural Mineral Tetrahedrite as a Direct Source of Thermoelectric Materials. Phys. Chem. Chem. Phys. 2013, 15, 5762−5766. (264) Lu, X.; Morelli, D. T. Rapid Synthesis of High-Performance Thermoelectric Materials Directly from Natural Mineral Tetrahedrite. MRS Commun. 2013, 3, 129−133. (265) Lu, X.; Morelli, D. T.; Xia, Y.; Ozolins, V. Increasing the Thermoelectric Figure of Merit of Tetrahedrites by Co-Doping with Nickel and Zinc. Chem. Mater. 2015, 27, 408−413.

12149

DOI: 10.1021/acs.chemrev.6b00255 Chem. Rev. 2016, 116, 12123−12149