Pingjun
Ying
a,
Lennart
Wilkens
a,
Heiko
Reith
a,
Nicolas Perez
Rodriguez
a,
Xiaochen
Hong
ab,
Qiongqiong
Lu
a,
Christian
Hess
ab,
Kornelius
Nielsch
*acd and
Ran
He
*a
aLeibniz Institute for Solid State and Materials Research, 01069 Dresden, Germany. E-mail: k.nielsch@ifw-dresden.de; r.he@ifw-dresden.de
bSchool of Mathematics and Natural Sciences, University of Wuppertal, 42079, Wuppertal, Germany
cInstitute of Materials Science, Technische Universität Dresden, 01062 Dresden, Germany
dInstitute of Applied Physics, Technische Universität Dresden, 01062 Dresden, Germany
First published on 26th April 2022
The applications of thermoelectric (TE) technology around room temperature are monopolized by bismuth telluride (Bi2Te3). However, due to the toxicity and scarcity of tellurium (Te), it is vital to develop a next-generation technology to mitigate the potential bottleneck in raw material supply for a sustainable future. Hereby, we develop a Te-free n-type compound Mg3Sb0.6Bi1.4 for near-room-temperature applications. A higher sintering temperature of up to 1073 K is found to be beneficial for reducing the electrical resistivity, but only if Mg is heavily overcompensated in the initial stoichiometry. The optimizations of processing and doping yield a high average zT of 1.1 in between 300 K and 573 K. Together with the p-type MgAgSb, we demonstrate module-level conversion efficiencies of 3% and 8.5% under temperature differences of 75 K and 260 K, respectively, and concomitantly a maximum cooling of 72 K when the module is used as a cooler. Besides, the module displays exceptional thermal robustness with a < 10% loss of the output power after thermal cycling for ∼32000 times between 323 K and 500 K. These proof-of-principle demonstrations will pave the way for robust, high-performance, and sustainable solid-state power generation and cooling to substitute highly scarce and toxic Bi2Te3.
Broader contextThermoelectric (TE) technology – whereby heat is converted to electrical power and vice versa – holds great potential for cooling and power generation in many applications because its unique solid-state nature enables TE devices to be free from emissions and maintenance, thus providing extraordinary reliability. Realizing this potential requires developing modules that show high performance around room temperature (200 K to 600 K). Currently, almost all commercial modules are based on bismuth telluride (Bi2Te3) because of its unparalleled performance to date. However, Bi2Te3 cannot meet the rapidly increasing demand for TE technology, because tellurium (Te) is very scarce, with the Earth's crust having a concentration of <0.001 ppm. A sustainable solution that avoids the usage of the toxic and scarce element Te is currently missing for commercial TE modules, thus preventing the realization of the full potential of TE technology. This work tackles this challenge by generating a new paradigm for the realization of Te-free, sustainable, high-performance, and robust TE materials and modules using scalable processing approaches, thereby significantly advancing the sustainability of TE technology for both cooling and power generation applications. |
The performance of TE technology, either for cooling or power generation, is largely determined by the figure-of-merit (zT) of the materials, zT = [(S2/ρ)/(κ)]T, where S, ρ, κ, and T are the Seebeck coefficient, electrical resistivity, thermal conductivity, and absolute temperature, respectively. The current state-of-the-art TE materials have shown average zT (zTave) values in the range of 0.8 to 1.2, which can be found in (1) Bi2Te3 around room temperature (200 K to 600 K); (2) skutterudites, Mg2Si, and IV–VI compounds at intermediate temperatures (600 K to 800 K); (3) SiGe alloy and half-Heusler compounds at higher temperatures (800 K to 1200 K).5 Many good TE materials, such as FeSi2, SiGe, and half-Heusler compounds, have peak zT values of 0.6 to 1.5 between 900 K and 1200 K, but only 0.1 to 0.3 at room temperature.6–9 Despite their exceptional robustness and applicability for industrial waste heat recovery, these materials are only suitable for heat sources with much higher temperatures, e.g. >700 K. The application of TE devices around room temperature is still dominated by Bi2Te3 for more than half a century. However, tellurium (Te) is an extremely scarce element with an abundance of less than 0.001 ppm (<0.0000001%) in the Earth's crust and an annual production of less than 500 metric tons globally.10,11 This greatly challenges the near-future growth of the thermoelectric market, and it is imperative to develop TE modules that are free from the scarce and toxic elements while retaining a high performance at around room temperature.
In recent years, Mg3(Sb,Bi)2 intermetallic (so-called Zintl) compounds have attracted great attention because of their non-toxicity and elemental abundance. These compounds were studied as the p-type material for more than 10 years until Tamaki et al. in 2016 first reported the realization of the n-type material by adding excess Mg.12 Subsequent investigations boomed, and a plethora of tailoring strategies have been proposed to advance their thermoelectric performances such as doping by various elements.13 Besides, alloying Mg3Bi2 with Mg3Sb2 not only introduces point defects that can effectively scatter phonons but also modifies the band structure such as the bandgap and the energy positions of several pockets of the conduction band.14 As a rule of thumb, a higher concentration of Bi shifts the peak zT towards a lower temperature, thus favoring applications at near room temperature. Particularly, one important strategy for improving the TE properties of Mg3(Sb,Bi)2 compounds is to reduce the grain-boundary scattering of electrons, which, in general, has been realized through (1) single crystal growth;15 (2) post-annealing after sintering under a Mg-vapor atmosphere;16 (3) introducing secondary phases (mostly elemental metals) to assist the grain growth.17 These approaches successfully mitigated the grain boundary scattering of electrons to certain degrees, so the zT improved accordingly, especially at lower temperatures. However, the growth of single crystals and Mg-vapor annealing are usually time-consuming and can be complex, thus unsuitable for a scaled preparation of materials to facilitate their applications. Moreover, the inclusion of the secondary phase complicates microchemistry, and the mechanism of how these inclusions assist the grain growth remains unclear.
A much more straightforward and scalable strategy to increase the grain size is to apply a higher temperature during sintering, which has been reported to effectively enhance the charge carrier mobility and zT of Mg3(Sb,Bi)2 compounds. However, this approach has been dominantly applied to Sb-rich compounds, mostly Mg3+δSb1.5Bi0.5.13,18 It remains less studied whether this approach is suitable for Bi-rich compounds. This is not a trivial question since the electrical properties of Mg3.2Bi1.5Sb0.498Te0.002Cu0.01 were completely deteriorated by applying a sintering temperature of 1073 K, as reported by Liu et al.19 However, 1073 K is usually a “safe” temperature for sintering Sb-rich compounds.13,18,20 It is thus speculated that Bi-rich compounds are more prone to Mg evaporation upon heating, as can be rationalized by the fact that the melting point of Mg3Bi2 (∼1094 K) is much lower than that of Mg3Sb2 (∼1500 K), so the Mg atoms are more volatile in Bi-rich compounds than in Sb-rich ones. Therefore, it remains a great challenge to establish processing that can sinter Bi-rich Mg3(Sb,Bi)2 compounds at much higher temperatures while minimizing the Mg losses in the sintered body to enable an advanced TE performance.
Hereby, inspired by the work of Wood et al. where the applied Mg-vapor annealing could significantly improve the electrical transport,16 and by the study of Imasato et al. where an n-to-p Seebeck transition was enabled by Mg-vapor annealing,15 we propose to heavily compensate the Mg losses during sintering by adding 20% extra Mg in the initial stoichiometry, yielding a base compositional formula of Mg3.6Sb0.6Bi1.4. This is different from many reports where the extra Mg did not exceed 10%. The high amount of extra Mg will compensate for its evaporation during sintering, which is similar to Mg-vapor annealing but realized in a single step. We found that this way, a higher sintering temperature of up to 1073 K successfully improves the electrical transport properties. Interestingly, we also found that such a high sintering temperature is unfavorable for a composition where the concentration of Mg is only 3.2 (i.e., Mg3.2Sb0.6Bi1.4), since the electrical performance degrades greatly. Using this synthesis approach, the TE properties are further optimized by changing the doping concentration. The peak zT values are enhanced to 0.74 at 300 K and 1.24 at 453 K for Mg3.6Y0.003Sb0.6Bi1.4, resulting in a high average zT of 1.1 between 300 K and 573 K. The improved zT in this work is comparable or even exceeds the state-of-the-art n-type Bi2Te3 and other Bi-rich Mg3(Sb,Bi)2 compounds.19,21,22 Subsequently, the optimized TE materials are translated into modules together with the p-type MgAgSb compound. The fabricated modules show excellent performances for both power generation and cooling. As shown in Fig. 1, for power generation, the conversion efficiencies reach 3% and 8.5% under temperature differences of 75 K and 260 K, respectively. For cooling, the maximum cooling ΔTs reach 52 K, 63 K, and 72 K with heat-rejection sides at 302 K, 325 K, and 347 K, respectively. These are comparable to or exceed the state-of-the-art modules.19,22–43 Furthermore, we perform thermal cycling for a Mg3(Sb,Bi)2/MgAgSb module by exposing the hot side between 323 K and 500 K continuously while maintaining the cold side at 288 K. We found that the loss in the relative output power is less than 10% after ∼32000 times cycling (∼5000 hours equivalence). The cycling performance outperforms the Bi2Te3-based module44 and unveils the great thermal robustness of our module. Our realizations are the first demonstration of an exceptionally robust thermoelectric module that is completely free from toxic and scarce elements for both high-performance power generation and cooling. The Mg-based technology thus promises the great potential to break the room-temperature monopoly of Bi2Te3 that has lasted for more than half a century.
Fig. 1 Module performance for (a) power generation, (b) cooling, and (c) internal resistance (left) and the relative output power (right) upon thermal cycling based on Te-free compounds in this work including p-type MgAg0.97Sb0.99 and n-type Mg3.6Y0.003Sb0.6Bi1.4. The obtained maximum conversion efficiency (ηmax), the cooling temperature (ΔTmax), and the relative output power are compared to other published results.19,22–44 The abbreviation Sku represents skutterudite compounds. The cycling test of Bi2Te3 is based on a commercial Melcor HT9-3-25 module, where the cold side was fixed at 297 K and the hot side varied between 300 K and 473 K.44 |
Fig. 2b shows that the Seebeck coefficient changes slightly with respect to the sintering temperature. The use of a higher sintering temperature, to 1073 K, is thus demonstrated to be effective for improving the power factor (S2ρ−1) by ∼46% at room temperature when compared to the one sintered at 973 K, as shown in Fig. 2c. Notably, such a high sintering temperature is applicable only if Mg is heavily compensated in the initial composition (Mg3.6Y0.003Sb0.6Bi1.4). This is clarified by sintering another compound at 1073 K with less excess Mg, i.e., Mg3.2Y0.003Sb0.6Bi1.4. We found a drastically higher resistivity together with a greatly reduced Seebeck coefficient if the amount of Mg is only 3.2, as shown in Fig. 2. These results manifest the Mg-poor conditions – a known performance killer for the n-type Mg3(Sb,Bi)2. The beneficial effects of the heavily compensated Mg are thus confirmed in maintaining Mg-rich conditions and in improving the electronic conductivity. The improving mechanism in our work is similar to that of Wood et al., where a higher conductivity was realized through post-annealing in Mg vapor for 65 hours.16 Despite the interesting findings, post-annealing not only adds an extra step to complicate and prolong the processing, but may also have insignificant effects if the sintered wafers are thick or if they are sandwiched between contact layers. In this work, we realize a similarly improved conductivity in just one step that will not encounter the thickness or contact issues. This aspect is important for the development of high-performance n-type Bi-rich Mg3(Sb,Bi)2 compounds on a larger scale with higher reproducibility to facilitate their applications.
The variations in the electronic transport properties are correlated with the microstructures. The X-ray diffraction patterns of the as-prepared specimens in this work are shown in the ESI† (Fig. S1). It is surprising to see that the addition of such a high amount of excessive Mg (∼20%) does not introduce impurity peaks. In fact, the peaks of a commonly observed impurity phase, Bi, do not appear in the XRD patterns in this work. This can be explained by the employed higher sintering temperatures that promote the reaction among the constituent elements, and the excessively added Mg guarantees the elemental supply despite its evaporation. In addition, Fig. 2d, e, and f show the scanning electron microscopy (SEM) images of the freshly broken surfaces of the Mg3.6Y0.003Sb0.6Bi1.4 samples sintered at 973 K, 1023 K, and 1073 K, respectively. The enlargement of the grains is found, which is responsible for the reduced electrical resistivity, especially at near room temperature. The observed grain growth stems from an increased atomic diffusion due to the higher sintering temperature. Despite the grain enlargement, the sample with the largest grains in this work still displays a ρ ∼ T0.5 relationship at near room temperature, which indicate the possibility of a subsequently improving the electrical transport upon further grain growth by, for example, adding more Mg in the initial composition and applying a higher sintering temperature. Indeed, a ρ ∼ T1.5 relationship has been found even at room temperature in a grain-boundary-free sample, i.e., a single crystal.14
Fig. 3 Thermoelectric transport properties of Mg3.6YxSb0.6Bi1.4 with x = 0.001, 0.002, 0.003, 0.004, and 0.005. (a) Temperature-dependent (left) electrical resistivity and (right) Seebeck coefficient. (b) Temperature-dependent (left) carrier concentration and (right) mobility. (c) Pisarenko relationship (i.e., Seebeck coefficient vs. the carrier concentration) at 300 K. (d) Temperature-dependent power factor. (e) Temperature-dependent (left) thermal conductivity and (right) zT. (f) Analysis of the (left) power factor and (right) zT at 300 K using the single parabolic band (SPB) model. The power factor and zT are compared to literature studies.19,21,22 |
The decrease in the resistivity is accompanied by the decrease in the Seebeck coefficient, as shown in the right axis in Fig. 3a. Besides, the temperature where the Seebeck coefficient peak appears displays a decreasing trend with a lower content of Y. This indicates an earlier onset of bipolar excitation with the decreased doping level. To understand the variations in both the resistivity and the Seebeck coefficient, we performed the Hall measurement for the Y-doped compounds from room temperature up to 468 K. The temperature-dependent Hall carrier concentration (nH) is found in Fig. 3b. The enhanced Y concentration yields an increase in nH. Furthermore, as is more clearly shown in the ESI† (Fig. S2), the nH curves display an increasing trend for x = 0.001 and x = 0.002 that is correlated with bipolar excitation. Interestingly, while nH is nearly constant for x = 0.003 and x = 0.004, it displays an anomalous decreasing trend for x = 0.005 upon heating. We note that a similar decrease in nH was also reported for the Sb-rich compound, Mg3.2Sb1.5Bi0.49Te0.01.48 The decreased nH with a higher temperature is not very well understood but might be related to the thermodynamic equilibrium concentration of Mg vacancies that depends on the temperature.48 This also agrees with the prerequisite of Mg-rich growth conditions to enable a higher doping efficiency for n-type Mg3(Sb,Bi)2. Furthermore, based on the measured nH and electrical resistivity (ρ), the temperature-dependent Hall mobility (μH) is evaluated using μH =1/ρnHe (e is the positive value of the electron charge). The μH shows a weak dependence on the dopant concentration, as shown in Fig. 3b. Besides, μH decreases with a temperature from ∼180 cm2 V−1 s−2 at 300 K to ∼120 cm2 V−1 s−2 at 450 K due to the acoustic phonon scattering that intensifies upon heating.
To further understand the electronic transport properties, we evaluate the density-of-state effective masses (mDOS) of the charge carriers by plotting the Seebeck coefficient vs. nH, i.e., through the Pisarenko relationship, , where kB is the Boltzmann constant and h is the Planck constant. Note that although this expression is assumed for carriers under the acoustic phonon scattering, it also applies in the presence of grain boundary scattering since both the Seebeck coefficient and the carrier concentration have weak dependences on the grain sizes, at least for the specimens where the grains are not extremely refined. The insensitivity of the carrier concentration and the Seebeck coefficient to the grain boundary scattering was reported for various thermoelectric materials such as Mg3.2Sb1.5Bi0.5 and half-Heusler NbFeSb.13,49 A gradual increase in mDOS from ∼0.54 m0 to ∼0.83 m0 is found by increasing x from 0.001 to 0.005, as shown in Fig. 3c. Note that mDOS = mb (NV)2/3, where mb is the band effective mass and NV is the number of valley degeneracy in the conduction band. Since the mobility shows an insignificant variation with the doping level (except for x = 0.001 at higher temperatures where the single-band assumption is less valid), mb experiences minor change and can be reasonably assumed as the constant. Then, the increased mDOS can be solely ascribed to an enhanced NV by a factor of ∼2, as shown in the ESI† (Fig. S3), indicating a gradual involvement of higher-energy pockets for electron transport when the Fermi level is moved deeper inside the conduction band. This is consistent with the calculation of Imasato et al. where the conduction band of Mg3(Sb,Bi)2 possesses a multi-valley nature, and the differences in the band-edge energies of these pockets are reduced by mixing Bi and Sb.14
Using the electrical resistivity and the Seebeck coefficient, the temperature-dependent power factors (PF, S2ρ−1) of Mg3.6YxSb0.6Bi1.4 are evaluated and shown in Fig. 3d. An improvement in the PF is realized by optimizing the doping level. For example, the PF is ∼1.6 mW m−1 K−2 with x = 0.001 at room temperature; it increases to 2.6 mW m−1 K−2 for x = 0.004, and further reaches a peak value of ∼2.8 mW m−1 K−2 at 363 K. We also plot the PF from other reports among the Bi-rich n-type Mg3(Sb,Bi)2 compounds that are more suitable for applications at near room temperature. In general, the PF realized in this work is better than others except for the recent work of Liu et al. where a slightly higher PF at lower temperatures was achieved.19 This is ascribed to the addition of Cu that effectively assisted the grain growth during sintering, which should also be complementary to the heavy-compensation strategy in this work and will be studied in the future. On the other hand, the achieved PF in this work is higher than that of Liu et al. at above 393 K, so the average PF is similar. In addition, Fig. 3e shows that the thermal conductivity displays an increasing trend by increasing the concentration of Y, which originates from the electronic thermal conductivity due to the reduced electrical resistivity.
The thermal conductivity values of the materials prepared in this work are small, with the lowest value reaching ∼0.8 W m−1 K−2 for Mg3.6Y0.001Sb0.6Bi1.4. Ultimately, the zT value in this work has been improved with the highest values reaching 0.74 at 300 K and 1.24 at 453 K for Mg3.6Y0.003Sb0.6Bi1.4, yielding a high average zT of 1.1 in between 300 K and 573 K. Fig. 3f shows the nH-dependent PF and zT at 300 K using a single parabolic band (SPB) model.50 Note that the modeling assumes a constant effective mass of 0.7 m0 that is averaged among these compounds, and that the charge carriers are scattered by the acoustic phonons only. We found that the experimental zT value agrees well with the SPB model if the carrier concentration is sufficiently high (nH >2 × 1019 cm−2). For the less doped compounds, the experimental zT value is slightly lower than the SPB-model predictions. These results manifest the impacts of the grain boundary scattering and suggest that further zT improvement is feasible by eliminating the grain boundary scattering in a compound with a slightly lower carrier concentration. In comparison to the n-type Bi2Te3 and Bi-rich n-type Mg3(Sb,Bi)2, our work achieves a comparable zT at room temperature,19,21,22,51 showing its great potential for cooling applications. Besides, at temperatures higher than 393 K, our material displays a better zT than many other reports, suggesting its potential to enable a higher conversion efficiency when integrating the as-prepared materials into power-generation and cooling modules.
Fig. 4a shows the measured current–voltage (I–V) curves for the power generation modules under various temperature differences (ΔT). A good linearity is obtained and based on which the module's internal resistance and the open-circuit voltage are evaluated using the slope and intercept, and the results are given in the ESI† (Fig. S7). P is plotted in the right axis of Fig. 4a, where the peak corresponds to the matched impedance condition where the internal resistance and the load resistance are equal. Moreover, the measured current-dependent heat flow (Qmeas.) at a series of ΔT, as well as the simulated temperature distribution, is plotted in the ESI† (Fig. S8 and S9, respectively). The Qmeas. values are only slightly higher (∼5%) than the simulated heat flow (Qsim.), indicating that our module design successfully minimized the thermal radiation energy losses. Subsequently, the current-dependent heat-to-power conversion efficiency (η) is evaluated under a series of ΔT, as shown in Fig. 4b. The ηmax value increases with the ΔT, and a peak ηmax of 8.5% is achieved under a ΔT of 260 K, as shown in Fig. 1a. The obtained efficiency exceeds all the Bi2Te3-based modules that are openly published to the authors’ knowledge. More interestingly, due to the improved zT at near room temperature, we achieved a high efficiency of 3% at a mild ΔT of 75 K. Hence, our obtained module is ideal for application using low-grade heat sources where the heat flow is small.
Because of the module's excellent thermoelectric performance at near room temperature, we further investigate its cooling performance. The theoretical maximum cooling performance is evaluated based on the finite-element simulation (FEM) using COMSOL. The temperature distribution of the modules under the cooling mode is shown in the ESI† (Fig. S9). In the future, we will use FEM to further optimize the geometrical configurations to advance the module-level performance, as has been similarly realized in designing power-generation modules.52 In addition, Fig. 4c gives the simulated coefficient of performance (COP) with a hot-side temperature of 325 K. The COP represents the amount of heat pumped divided by the amount of the supplied electrical power. Unfortunately, the experimental COP cannot be evaluated in this work due to the lack of heat-flow-measurement capability in our measurement setup. The measurement setup will be upgraded in the future with the installation of heaters and heat-flux sensors. Nevertheless, the cooling ΔT can be simulated and measured, and compared, as shown in Fig. 4d. With the increasing current, the cooling ΔT increases until it reaches a maximum and then decreases again. Such a current-dependent ΔT curve stems from the interplay among several effects including heat conduction, the Peltier cooling, the Joule heating, and the Thomson effect. The theoretical maximum cooling ΔT reaches 58 K, 70 K, and 80 K with heat-rejection sides at 302 K, 325 K, and 347 K, respectively, which are slightly higher than the experimental numbers that are 52 K, 63 K, and 72 K, respectively. Considering the high vacuum (∼10−7 mbar) and the relatively low temperature, heat losses from radiation or convection should be negligible. Therefore, we speculate the experimentally reduced cooling ΔT has two possible origins: (1) electrical losses between the contact/thermoelectric materials and the brazing points for wiring. (2) The inaccurate simulation since the low-temperature transport properties are not directly measured, instead they are extrapolated using the data from room temperature and above. In the future, we will address these issues to further understand and boost the cooling performance by, e.g., measuring and improving the low-temperature properties of materials. Nevertheless, the realized cooling is comparable to commercial Bi2Te3-based modules,22 thus showing great promise in developing the Te-free thermoelectric technology for near-room-temperature cooling applications.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2ee00883a |
This journal is © The Royal Society of Chemistry 2022 |