Non-equilibrium processing leads to record high thermoelectric figure of merit in PbTe–SrTe

The broad-based implementation of thermoelectric materials in converting heat to electricity hinges on the achievement of high conversion efficiency. Here we demonstrate a thermoelectric figure of merit ZT of 2.5 at 923 K by the cumulative integration of several performance-enhancing concepts in a single material system. Using non-equilibrium processing we show that hole-doped samples of PbTe can be heavily alloyed with SrTe well beyond its thermodynamic solubility limit of <1 mol%. The much higher levels of Sr alloyed into the PbTe matrix widen the bandgap and create convergence of the two valence bands of PbTe, greatly boosting the power factors with maximal values over 30 μW cm−1 K−2. Exceeding the 5 mol% solubility limit leads to endotaxial SrTe nanostructures which produce extremely low lattice thermal conductivity of 0.5 W m−1 K−1 but preserve high hole mobilities because of the matrix/precipitate valence band alignment. The best composition is hole-doped PbTe–8%SrTe.

T hermoelectric materials technology can play a role in energy-saving alternatives for sustainable power generation, by providing a direct means for heat to electricity conversion [1][2][3][4][5][6] . The efficiency of a thermoelectric material is defined by the figure of merit ZT ¼ S 2 sT/k ¼ S 2 sT/(k el þ k lat ), where S is the Seebeck coefficient, s is the electrical conductivity, T is the temperature and k is the thermal conductivity that is comprised of an electronic part (k el ) and a lattice part (k lat ). PbTe-based materials are the top-performing thermoelectrics in the temperature range of 500-900 K (refs 7-11), thanks to a unique two-valence-band structure: a primary light hole band at the L point (L band) and a second lower-lying heavy hole band at the S point (S band) with an energy separation (DE L À S ) of 0.15-0.20 eV at 300 K between them 12,13 . Recent studies [14][15][16] show that Mg or Mn substitution in PbTe can decrease DE L À S to within a few k B T (k B is the Boltzmann constant). This decrease is referred to as effective band convergence, which makes the contribution of the S band to the charge transport very significant, especially when temperature is higher than 500 K, and produces much higher Seebeck coefficients 9,17 .
In addition to band structure engineering, nanostructuring is an efficient route to enhancing the ZT of PbTe by strongly reducing the thermal conductivity. For instance, previous studies showed only a small equilibrium solid state solubility (less than1 mol%) of SrTe in bulk PbTe (refs 7,8). Samples prepared under thermodynamic control, using slow cooling procedures of molten mixtures (11 K h À 1 from 1,323 to 873 K (ref. 7), dissolve only very small amounts of SrTe. Any excess SrTe emerges as strained endotaxial nanostructures whose valence band is almost aligned with that of PbTe. This matrix/ precipitate valence band alignment mitigates the scattering of holes, while the added interfaces created by SrTe nanostructuring strongly impede phonon propagation 7 . When coupled with the atomic scale point defects created by Na doping and mesoscale structuring using spark plasma sintering (SPS), these combined effects can lead to a ZT of B2.2 at 913 K (ref. 8).
Although the thermodynamic solubility limit of SrTe in PbTe is low (o1 mol%), early studies of PbTe-SrTe thin films 18,19 grown under kinetic control using molecular beam epitaxy showed that up to 15 mol% of SrTe can be dissolved into PbTe. If this higher SrTe alloying fraction can be realized in bulk PbTe, our calculations predict it would lead to valence band convergence, an effect which is negligible in the low solubility samples 7,8 . It should be noted that in previous studies on PbTe-SrTe thin films 19 , no such a valence band convergence effect was observed because the thin films were n-type doped. The theoretical prediction motivated us to pursue kinetically controlled methods for fabricating PbTe samples heavily alloyed with SrTe. Herein, we show that via a non-equilibrium route that involves rapid ice water quenching (estimated cooling rate higher than 100 K/s, based on the time it took before the melt was completely cooled to room temperature, which enables the fine element distribution) of the 1,423 K melt followed by annealing at 873 K (which serves as kinetic control) we can produce hole-doped samples that achieve a 5 mol% solubility of SrTe in bulk PbTe. Exceeding the 5 mol% limit by using excess SrTe then leads to ubiquitous nanostructures.
We show here that the high degree of SrTe alloying produces a material that simultaneously incorporates four practical performance-enhancing mechanisms: valence band convergence of PbTe, a strengthened point defect phonon scattering when SrTe is alloyed, and extensive nanostructuring with good valence band alignment of precipitate/matrix when excess SrTe is added. The result is a highly optimized material with a high ZT of B2.5 at 923 K coupled with a new milestone of average ZT of 1.67 between 300 and 900 K in Na-doped PbTe-8%SrTe. This new approach is quite different from previous studies on PbTe (for example PbTe 0.7 S 0.3 and (PbTe) 0.65 (PbS) 0.25 (PbSe) 0.1 (refs 20,21) where only a few of the above mechanisms are combined. We hope this approach can be equally applied to other thermoelectric materials for improved performance.

Results
SrTe alloying in PbTe and electronic structure modification. Within the detection limit of laboratory X-ray diffraction, all the Pb 0.98 Na 0.02 Te-x%SrTe samples in this study, prepared under non-equilibrium conditions, are single phase compounds crystallizing in the rock salt structure (Supplementary Fig. 1). The lattice parameters shown in Fig. 1a increase almost linearly with increasing SrTe content up to x ¼ 6 and beyond that remain unchanged. The lattice expansion of PbTe upon SrTe addition is consistent with the larger radius of Sr 2 þ (1.32 Å) in comparison to Pb 2 þ (1.19 Å) and confirms the successful substitution of Sr for Pb in PbTe. We also note that water quenching produces a behaviour of the lattice parameter ( Supplementary Fig. 2 band structure ( Supplementary Fig. 3) similar to that obtained by molecular beam epitaxy 18,19 , which further supports the effectiveness of non-equilibrium synthesis in producing a high degree alloying fraction of SrTe in PbTe. By fitting the lattice parameter as a function of SrTe content (on the basis of Vegard's law) one can conclude that the achieved solubility limit of SrTe in these samples is around 5 mol%. This value is significantly higher than the previously achieved o1 mol% under equilibrium synthesis conditions 7,8 . Infrared absorption spectra for low carrier density PbTe-x%SrTe samples show that the electronic absorption edge shifts towards higher energy with x increasing from 0 to 6 and remains unchanged afterwards (Fig. 1b). Figure 1c shows that the optical bandgap increases from 0.23 eV in pristine PbTe to 0.34 eV for samples with x higher than 5. The larger bandgaps are consistent with successful alloying of the much wider bandgap semiconductor SrTe (E g ¼ 3.5 eV, ref. 22) in PbTe.
The presence of band convergence in the PbTe-SrTe systems was probed with Hall effect measurements (Fig. 1d). Typically, for semiconductors with only a single band contributing to the charge transport, the Hall coefficient (R H ) is nearly temperature independent. The case is very different for p-type lead or tin chalcogenides featuring two or more valence bands 14,17,[23][24][25][26] . In these compounds, as the temperature rises, the L band moves to lower energy while the S band remains almost unchanged, resulting in a smaller DE L À S and carrier redistribution between the two sub-bands 14,17,26 . This change in band structure is reflected in the strong temperature dependence of R H which peaks at a temperature when the contribution from the two bands to hole transport is comparable 23,26,27 . The temperature (T H ) corresponding to the peak of R H is a measure of DE L À S and is greatly affected when the band structure is altered for example by alloying 17,23,26 . For pure p-type PbTe with different hole concentrations, R H always peaks at B425 K (refs 14,17; Fig. 1d). However, when PbTe is alloyed with MgTe or MnTe, DE L À S is diminished which is reflected in the lowered T H (refs [14][15][16]. In our current study, we see a similar decreasing trend of T H as the SrTe content is increased ( Supplementary Fig. 4), indicating that the addition of SrTe decreases DE L À S , which is consistent with our theoretical calculations presented below. On the contrary, in our previous report on PbTe-SrTe using equilibrium synthesis, R H peak temperature shows negligible changes with increasing SrTe fraction because of the very low solubility of SrTe in PbTe (ref. 7). Therefore, the high-temperature Hall data serves as another diagnostic to confirm the high solubility of SrTe with PbTe achieved in this work.
We also explored the band structure modification of PbTe by alloying with SrTe using first-principles density functional theory (DFT) electronic structure calculations. Since Na serves only as a dopant in PbTe adjusting the carrier concentration, as has been demonstrated by numerical experiment studies and theoretical calculations [28][29][30] , for the simplicity of calculation, we showed the band structure calculation results of PbTe-SrTe without Na. The DFT band structures of PbTe-SrTe with and without spin-orbit interaction (SOI) are shown in Fig. 1e and Supplementary Fig. 5, respectively, both of which exhibit qualitatively similar trends for E g and DE L À S ( Supplementary Fig. 6). When the Pb atoms are partially replaced by Sr, the conduction band minimum at the L point of PbTe remains roughly constant in energy while both the L and S valence bands move away from the conduction band, leading to an enlarged E g . As the L band decreases in energy much faster than the S band, DE L À S decreases, Fig. 1f. The larger E g and smaller DE L À S favourably impact the thermoelectric power factors of the PbTe-x%SrTe system, as we present below.
For all samples, the electrical conductivity ( Fig. 2a) decreases while the Seebeck coefficient ( Fig. 2b) increases with increasing temperature, typical of degenerate conduction. Because of the similar hole concentrations (Supplementary Table 1), the high-temperature electrical conductivity for all samples is similar as the SrTe content is increased from 0 to 12%; in contrast, the Seebeck coefficient shows a considerable enhancement with increasing SrTe. Specifically, the Seebeck coefficient is 62 mV K À 1 for x ¼ 0 and goes up to greater than 90 mV K À 1 for x46 at 300 K (Supplementary Table 1). Also, the Seebeck coefficient of the Pb 0.98 Na 0.02 Te control sample shows a clear downturn B800 K attributed to bipolar diffusion 31 . In contrast, the SrTecontaining samples, show continuously increasing Seebeck coefficient up to 923 K suggesting negligible bipolar conduction. This behaviour is consistent with the enlarged bandgap ( Fig. 1c) which inhibits the thermal activation of minority carriers at elevated temperatures.
The room temperature Hall carrier concentrations (N p ) and hole mobilities (m H ) for all samples are given in Supplementary  Table 1. We estimated the carrier effective mass (m*) of each sample from the measured S and N p considering a single parabolic band 32 for simplicity. The m* increases gradually from 0.82 m e to 1.42 m e (m e is the free electron mass) with increasing SrTe content from 0 to 6 mol%. We have proposed above that the alloying of SrTe causes convergence between the two valence bands (L and S) of PbTe. The degeneracy (N v ) of the L band of PbTe is 4 while that of the S band is much larger, N v ¼ 12 (refs 9,17). When the L and S bands are sufficiently close in energy (for example, by SrTe alloying), the effective N v becomes 12-16 rather than 4 for the L band or 12 for the S band 9,14 .  convergence effect (increase of N v ) significantly increases m* by employing more degenerate valleys in the hole transport without seriously deteriorating m H (Supplementary Table 1). The valence band convergence caused by the much higher SrTe alloying achieved with the non-equilibrium processing is reflected in the enhancement of the Seebeck coefficient as shown by the Pisarenko relation between S and N p (Fig. 2c). The solid line is the established theoretical Pisarenko plot 14,34 for PbTe. We observe that the Seebeck coefficients for our heavily Sr-alloyed PbTe samples are much higher than predicted by the solid line for PbTe and behave very similar to Mg-or Mn-doped PbTe (refs 14-16; Supplementary Fig. 7). On the contrary, the Seebeck coefficients of the lightly alloyed PbTe-SrTe samples of Biswas et al. 7,8 prepared using equilibrium synthesis fall exactly on the Pisarenko plot for PbTe.
These results further demonstrate that much higher amounts of SrTe successfully alloyed into PbTe contribute to the marked enhancement of the Seebeck coefficient over a broad temperature range (Fig. 2b). Therefore, the power factors (S 2 s) of the Pb 0.98 Na 0.02 Te-x%SrTe samples in this study are remarkably enhanced, with the maximum values exceeding 30 mW cm À 1 K À 2 around 500 K (Fig. 2d). These values are the highest reported for any p-type PbTe (refs 7,8,14-16; Supplementary Fig. 8), including the equilibrium synthesized Pb 0.98 Na 0.02 Te-x%SrTe by Biswas et al. 8 with similar hole densities (o25 mW cm À 1 K À 2 ). We also compared the temperature-dependent power factors of Na-doped PbTe alloyed with 3 mol% MgTe, MnTe and SrTe in this study (within the solid solution region where alloying-induced valence band convergence is taking effect; Supplementary Fig. 9). Although having a similar N p , the 3 mol% SrTe sample displays significantly higher power factor than the other two (note that their Seebeck coefficients at room temperature are quite similar, specifically, 80 mV K À 1 for 3%SrTe, 88 mV K À 1 for 3% MnTe (ref. 16) and 89 mV K À 1 for 3% MgTe (ref. 15)). The reason could be that the smaller radius contrast between Pb 2 þ (1.19 Å) and Sr 2 þ (1.32 Å) (for comparison, the radii of Mg 2 þ and Mn 2 þ are 0.86 and 0.83 Å, respectively) makes the scattering of carriers from crystal defects less significant. As an evidence of this conjecture, we list the room temperature carrier mobilities for the three samples as follows: 69, 65 and 78 cm 2 V À 1 s À 1 for 3% MgTe (ref. 15), MnTe (ref. 16) and SrTe alloyed PbTe, respectively. This additional advantage from crystal chemistry, along with the valence band convergence in the band structure by SrTe alloying, endows the PbTe-SrTe system with superior electrical properties.
Thermal transport and valence band alignment. Both the total (k, Fig. 3a) and lattice thermal conductivity (k lat , Fig. 3b) of PbTe are considerably suppressed upon the addition of SrTe. The 10% SrTe sample exhibits the lowest room temperature k lat of 1.5 W m À 1 K À 1 , which drops to 0.5 W m À 1 K À 1 at 923 K. However, the SrTe nanostructuring has little impact on the carrier mobilities in Pb 0.98 Na 0.02 Te-x%SrTe (x ¼ 0-12). Figure 3c shows hole mobilities (m H ) as a function of temperature for x ¼ 6, 8 and 10 samples (all are nanostructured) with similar hole concentrations of 1.6-1.7 Â 10 20 cm À 3 . Clearly the changes of m H with increasing SrTe amounts are negligible. This is attributed to the favourable energy alignment of the valence bands of PbTe and SrTe as reported earlier 7 . Moreover, for all samples, the mobilities follow a power-law temperature dependence of T À 1 from 300 to 400 K and T À 4 from 400 to 700 K. The fast degradation of mobilities at elevated temperature in p-type PbTe-based compounds arises from the contribution of the heavy holes in the S valence band.
For semiconductors with Umklapp scattering as the dominant phonon scattering mechanism, k lat should vary as 1/T before any bipolar diffusion becomes significant 25,26,32,35 , which is the case in PbTe (ref. 14). k lat starts to deviate from such a relationship when the contribution from bipolar diffusion (k bi ) is significant (Fig. 3d). Figure 3e shows the ratio of k bi to k lat at 923 K as a function of SrTe content for Pb 0.98 Na 0.02 Te-x%SrTe. The high content of SrTe strongly inhibits the bipolar diffusion because of the enlarged bandgap.
Using the Debye-Callaway model 36 we were able to simulate the lattice thermal conductivities of the PbTe-8%SrTe sample (see Supplementary Methods for simulation details and  Supplementary Table 2 where the related parameters for simulation are listed) and figure out the contributions from different scattering mechanisms to the lattice thermal conductivity by examining the relaxation time versus the normalized frequency plot (Supplementary Fig. 10). We find    that the simulated lattice thermal conductivities agree fairly with the experimental ones (Fig. 3b). In addition, the contribution from strain, induced by point defects, to the lattice thermal conductivity is more significant at lower frequency while nanoprecipitates and point defects are more important phonon scattering sources at medium phonon frequencies. We also modelled the room temperature lattice thermal conductivities of Pb 0.98 Na 0.02 Te-x%SrTe using a modified Klemens model 37,38 . This model takes into account both the mass and strain field contrasts and has been demonstrated extremely useful for disordered alloys 25,26,39,40  Microstructural and compositional information was obtained using scanning electron microscopy (SEM), (scanning) transmission electron microscopy ((S)TEM) equipped with energy dispersive spectroscopy and secondary ion mass spectrometry (SIMS). Figure 4a shows a SEM image on a freshly fractured surface of Pb 0.98 Na 0.02 Te-8%SrTe, with grain sizes ranging from 1 to 3 mm on average, as indicated in Fig. 4c. The low magnification TEM image presents high-density nanoprecipitates with darker contrast along the [110] zone axis, as shown in Fig. 4b. The inset shows a selected area electron diffraction pattern that exhibits single crystal diffraction pattern indicating the endotaxial relationship of the precipitates with the PbTe matrix. Moreover, Na, Sr, Pb and Te are all spatially uniformly distributed under SIMS ( Supplementary Fig. 11). In Fig. 4d, the size distribution of the nanoscale precipitates shows a peak of 7 nm along the [110] zone axis. The HRTEM image in Fig. 4e reveals coherent interfaces between precipitates and matrix, highlighted by dashed yellow rectangular box. These interfaces enable effective phonon scattering and minimize hole scattering, and account for the low thermal conductivity, high mobility and high-power factors in our samples.
Our valence band alignment calculations indicate that the valence band energy offset (DE V ) between PbTe and SrTe is only 0.06 eV with the valence band maximum (VBM) of PbTe lying slightly higher than that of SrTe (Fig. 4f). When Sr substitutes for Pb in PbTe, in addition to the band convergence discussed above, the VBM of PbTe is further lowered because of the formation of solid solution Pb 1 À x Sr x Te (Fig. 1e,f; Supplementary Fig. 5). A previous study on PbTe-MgTe by Pei et al. 15 suggested that a shift of the peak R H to lower temperature by 70 K (which is also the case in our study, Supplementary Fig. 4) represents an energy reduction of 0.03 eV between L and S bands. Applying this argument, we estimate that DE V is only 0.03 eV (nearly aligned) between SrTe precipitates and Pb 1 À x Sr x Te matrix, Fig. 4f. This matrix/precipitate band alignment preserves the high hole mobility in nanostructured systems 7,8,46,47 .
With increasing SrTe content dissolved in PbTe, we have demonstrated the highest power factors in PbTe-based systems while preserving the all-scale architectures that give low thermal conductivities. The combination of charge and thermal transport properties presented above leads to high thermoelectric performance for Pb 0.98 Na 0.02 Te-x%SrTe (see their ZT values in Supplementary Fig. 12), with a record high ZT of 2.5 at 923 K for x ¼ 8. This is nearly 80% improvement over the control sample x ¼ 0 (ZT ¼ 1.4) and 15% higher than the sample by Biswas et al. (Fig. 5).
We also note that the magnitude of ZT enhancement in this study with respect to the previous work is especially noticeable at To800 K while at elevated temperature (800-923 K) this enhancement becomes less significant. In PbTe the change of the energy levels of the L and S valence bands with increasing temperature is an unusual phenomenon which causes band convergence. According to recent studies 48,49 , in PbTe itself, the temperature-driven L-S valence band convergence occurs at 700-780 K. Simply applying this argument and considering that further SrTe alloying decreases the energy separation between the two bands ( Fig. 1e; Supplementary Fig. 4), one would expect a slightly lower than 700-780 K band convergence temperature as schematically shown in Supplementary Fig. 13. Therefore, at very high temperature (T4800 K) the thermoelectric performance of PbTe-SrTe system would not benefit from the temperaturedriven band convergence more than PbTe itself does. However, the enlarged bandgap (E g ) of PbTe by alloying with SrTe could play a significant role in improving the thermoelectric performance at elevated temperature (T4800 K), due to the suppression of bipolar diffusion. For example, the Seebeck coefficient of PbTe-SrTe shows no sign of downturn up to 900 K (Fig. 2b) while its bipolar thermal conductivity is largely decreased (Fig. 3c). When the temperature is not high enough to result in temperature-driven band convergence (for example, To600 K), the SrTe alloying-induced band convergence (decrease of DE L À S ) makes the contribution from the heavy S valence band to the Seebeck coefficient become noteworthy. As a result, one can observe a large enhancement of Seebeck coefficient of PbTe-SrTe as compared to pure PbTe, see Pisarenko relationship between Seebeck coefficients and carrier concentrations (Fig. 2c). Correspondingly, the maximum power factors achieved in our PbTe-SrTe samples exceed 30 mW cm À 1 K À 2 at around 600 K (Fig. 2d). They are the highest reported for p-type PbTe. In addition, the ZTs at low temperature are significantly enhanced (Fig. 5).
The average ZT (ZT ave ), which is important for device applications, in the range of 300-900 K, is 1.67 for the high performing sample in the present study. In comparison, ZT ave is 1.05 for PbTe and 1.37 for the lightly alloyed samples of Biswas et al. 8 The high performing Pb 0.98 Na 0.02 Te-8%SrTe sample showed negligible changes of thermoelectric properties during the multiple heating-cooling measurement cycles and after a 15 days of vacuum annealing at 823 K, suggesting an excellent thermal stability ( Supplementary Figs 14 and 15). More solid evidence of this stability comes from the unchanged lattice parameters and bandgaps before and after the vacuum annealing ( Supplementary Fig. 16). We believe that the non-equilibrium solid is in a meta-stable state, but the stability is not compromised at the highest temperatures that we probed in this study. Figure 6 illustrates the difference between equilibrium and non-equilibrium processing in the PbTe-SrTe system. The lower half of Fig. 6 shows how equilibrium processing gives low SrTe solubility (o1%) and leads to band alignment, nanostructuring but not band convergence. This route gives a ZT of 2.2 at 913 K. The upper part of Fig. 6 shows how non-equilibrium processing can give much higher dissolved fraction of SrTe (5 mol%). In addition to the band alignment and nanostructuring, the higher alloying fraction of SrTe converges the two valence bands (L and S) and enlarges the bandgap, leading to much higher power factors and improved performance. This tight integration of all these properties into a single material leads to record high ZT values for Pb 0.98 Na 0.02 Te-8%SrTe of 2.5 at 923 K. We further demonstrate that a peak ZT of 2.5 can be achieved at a relatively lower temperature around 800 K by optimizing the Na concentration ( Supplementary Fig. 17), making PbTe-SrTe system as a robust candidate for thermoelectric power generation.

Methods
Synthesis. High purity Pb wire (99.99%, American Elements, US), Sr chunk (99.9%, Sigma-Aldrich, USA), Na chunk (99.999%, Sigma-Aldrich, US) and Te shot (99.999%, 5N Plus, Canada) were used as the starting materials to synthesize 15 g of Pb 1 À y Na y Te-x%SrTe (y ¼ 0.02, x ¼ 0, 1.5, 3, 4.5, 6, 8, 10, 12; x ¼ 8, y ¼ 1%, 1.25%, 1.5%, 1.75; in mole ratio). Desired amounts of Pb, Sr, Na and Te were weighed and loaded into 13 mm diameter carbon coating silica tubes under an N 2 -filled glove box. The tubes were then evacuated to a residual pressure of B10 À 4 torr, flame-sealed, slowly heated up to 1,423 K in 11 h, soaked at this temperature for 12 h and subsequently ice water quenched to room temperature. The tubes containing the molten samples were periodically shaken to ensure the homogeneity of the compositions. The quenched ingots were further vacuum annealed at 873 K for 3 days. To probe into the bandgap variation of PbTe as a function of SrTe content using infrared spectroscopy (IR), a series of PbTe-x%SrTe ingots without Na doping (B5 g in mass for each, x ¼ 0, 2, 4, 6 and 8, in mole ratio) were also synthesized using the same route. All samples synthesized in this study are all single phase compounds within the detection limit of powder X-ray detection limit (Supplementary Figs 1, 18 and  19). For a typical experiment the following amounts were used: Pb (8.7440 g, 42.2010 mmol), Na (0.0198 g, 0.8612 mmol), Sr (0.3018 g, 3.4450 mmol) and Te (5.9343 g, 46.5072 mmol) were used to prepare 15 g of Pb 0.98 Na 0.02 Te-8%SrTe.
Densification. The annealed cast ingots of Pb 0.98 Na 0.02 Te-x%SrTe were ground into fine powders using agate mortar under an N 2 -filled glove box, put inside a 12.7 mm diameter graphite die and densified by SPS (SPS-211LX, Fuji Electronic Industrial Co., Ltd.) at 823 K for 5 min under an axial compressive stress of 60 MPa in vacuum. Highly dense disk-shaped pellets with dimensions of 12.7 mm diameter and B10 mm thickness were obtained (Supplementary Table 3).
Electron microscopy and X-ray diffraction. (Scanning) transmission electron microscopy (S/TEM) and STEM energy dispersive spectroscopy experiments investigations were carried out using a JEOL 2100F microscope operated at 200 kV. The thin TEM specimens were prepared by conventional methods, including cutting, grinding, dimpling, tripod, with minimal duration of Ar-ion milling with a liquid N 2 cooling stage. Samples pulverized with an agate mortar were used for powder X-ray diffraction. SIMS was performed with PHI TRIFT III ToF-SIMS. The sample was prepared by mechanical polishing on SiC paper with ethanol and then on cloth with oil-based SiC slurry. A 15 kV Gallium ion beam was used in the secondary ion mass spectrometry. After initial spectrum collection, Ga ion was sputtered on a 100 Â 100 mm surface area of the sample for 60 s to clean out the surface contamination. Spectrum and mapping images were collected using a raster size (size of the sputter area on the sample) of 10 Â 10 mm. The powder diffraction patterns were obtained with Cu K a (l ¼ 1.5418 Å) radiation in a reflection geometry on an Inel diffractometer operating at 40 kV and 20 mA and equipped with a position-sensitive detector. The lattice parameter was obtained using a Rietveld refinement method (Fig. 1a). Electrical properties. The obtained SPS-processed pellets were cut into bars with dimensions 12 Â 3 Â 3 mm 3 for simultaneous measurement of the Seebeck coefficient and electrical conductivity using an Ulvac Riko ZEM-3 instrument under a low-pressure helium atmosphere from room temperature to 923 K. Samples were spray coated with boron nitride spray to minimize outgassing except where needed for electrical contact with the thermocouples, heater and voltage probes. The uncertainty of the Seebeck coefficient and electrical conductivity measurements is 5%, which is widely accepted for Ulvac instruments 50 .
Thermal conductivity. Highly dense SPS-processed pellets were cut and polished into a square shape of 6 Â 6 Â 2 mm 3 for thermal diffusivity measurements. The samples were coated with a thin layer of graphite to minimize errors from the emissivity of the material. The total thermal conductivity was calculated from k tot ¼ D Á C p Á d, where the thermal diffusivity coefficient (D) was measured using the laser flash diffusivity method in a Netzsch LFA457, and the density (d) was determined using the dimensions and mass of the sample. Heat capacity (C p ) is estimated from the relation C p /k B per atom atom ¼ 3.07 þ (4.7 Â 10 À 4 Â (T À 300)) 15,16 , which is obtained by fitting the experimental data reported by Blachnik 51 within an uncertainty of 2% above room temperature. The thermal diffusivity data were analysed using a Cowan model with pulse correction. The measured densities of all the samples range between 7.9 and 8.0 g cm À 3 or are above 97% of the theoretical densities (Supplementary Table 3). The uncertainty of the thermal conductivity is estimated to be within 8% (ref. 50 ), considering all the uncertainties from D, C p and d. The lattice thermal conductivity (k lat ) was calculated by subtracting the electrical thermal conductivity (k el ) from k tot using a Wiedemann-Franz relationship k el ¼ L Á s Á T, where L is Lorenz number which can be obtained by fitting the Seebeck coefficient to the reduced chemical potential 32,52 . The values of C p , D and L values for all the samples in this study can be found in Supplementary Figs 20 and 21. The combined uncertainty for all measurements involved in the calculation of ZT is o15%. Unless otherwise noted, all the thermoelectric properties were measured perpendicular to the sintering pressure direction, although no directional anisotropy effects (o3%) were observed in the charge transport properties.
Hall measurements. The high-temperature Hall measurement was performed on a homemade apparatus (University of Michigan) in an argon atmosphere. The Hall resistance was monitored with a Linear Research AC Resistance Bridge (LR-700), with constant magnetic fields of ± 1 T applied by using an Oxford Superconducting magnet. The effective carrier concentration (N p ) was estimated using the relationship N p ¼ 1/eR H , where e is the elemental charge, and R H is the Hall coefficient. The Hall mobility (m H ) was calculated using the relationship m H ¼ sR H with s being the electrical conductivity obtained from the ZEM-3 instrument. It should be mentioned here that the use of two different instruments to deduce carrier mobility data may cause some uncertainties.
Infrared spectroscopy. Room temperature optical diffuse reflectance measurements were performed on finely ground powders to probe optical energy gap of the series. The spectra were collected in the mid-IR range (6,000-400 cm À 1 ) using a Nicolet 6700 FT-IR spectrometer. The reflectance versus wavelength data generated, were used to estimate the bandgap by converting reflectance to absorption data according to Kubelka-Munk equations: a/S ¼ (1 À R) 2 /(2R), where R is the reflectance, a and S are the absorption and scattering coefficients, respectively.
Band structure calculations. DFT calculations of pristine, stoichiometric PbTe, and Sr-doped PbTe were carried out. The calculations were performed using the generalized gradient approximation with PBE 53 functional for the exchange À correlation functional and projector augmented wave potentials as implemented in Vienna ab initio Simulation Package (VASP) 54 . All the atomic positions are relaxed until the forces exerted on the atoms are less than 0.001 eV A À 1 . The calculated lattice constant of PbTe is 6.43 Å, which is very close to experimental result 6.463 Å. All structures are fully relaxed with respect to cell vectors and cell-internal positions. For the Pb and Sr species, the 5d electrons and 4p 5s states are respectively treated as valence states. The total energies were numerically converged to B3 meV per cation using a basis set energy cutoff of 400 eV and dense k-meshes corresponding to 4,000 k-points per reciprocal atom in the Brillouin zone.
To investigate the movements of the conduction band and valence bands (L and S bands) with added Sr ions, we considered 3 Â 3 Â 3 supercells of NaCl-type Pb 27 Te 27 . For the isovalent doping of Sr, we consider a single Sr impurity (B3.7% and 7.4% additions of Sr, Pb 26 Sr 1 Te 27 and Pb 25 Sr 2 Te 27 ) with Sr substituting for Pb. The symmetry of the original primitive cell was changed by the substitution defects in PbTe, thus, for the purposes of a more direct comparison with PbTe we transformed the eigenstates for defect structures into a so-called effective band structure in the primitive Brillouin zone of the parent compound PbTe using a spectral decomposition method 55,56 . Using this approach, we are able to calculate the energy level of the L point and S line and the corresponding energy differences for the supercells with defects. Note that the band structures are calculated by considering both with SOI (Fig. 1e) and non-SOI ( Supplementary  Fig. 5), and similar trends were observed in both cases.
To assess the relative band alignments of the second phase SrTe and MgTe with host phase PbTe we utilize the findings of Van de Walle and Neugebauer, who demonstrated a universal alignment of the electronic transition level of hydrogen in a wide range of materials including semiconductors, insulators and even aqueous solutions 57 . Hence, to infer the band alignment, we compute the energies of H defects in the rock salt compounds of interest, assume alignment between these H energies, and then extract the band alignment of the compounds. To align the valance band maximum position of each system, we consider the defect formation energies of various charge states of interstitial H q (q ¼ À 1, 0, 1) by placing H in the host material with a 128 atom supercell, calculating the total energy of this structure, and subtracting the energy of the corresponding pure host material, hydrogen chemical potential, and electron chemical potential: where E V and E F are VBM and Fermi level (relative to the VBM). To select the most favourable interstitial H binding sites in host materials, multiple binding configurations are calculated. The electrostatic potential correction term DE is calculated by inspecting the potential in the supercell far from the impurity and aligning it with the electrostatic potential in bulk 58 .
Data availability. The authors declare that the data supporting the findings of this study are available within the article and its Supplementary Information files, or from the corresponding authors upon request.