Yu
Yan
a,
Wen
Zhang
a,
Huijun
Kang
*ab,
Zongning
Chen
ab,
Enyu
Guo
ab,
Rongchun
Chen
*a and
Tongmin
Wang
*ab
aKey Laboratory of Solidification Control and Digital Preparation Technology (Liaoning Province), School of Materials Science and Engineering, Dalian University of Technology, Dalian, 116024, China. E-mail: kanghuijun@dlut.edu.cn; crc618@dlut.edu.cn; tmwang@dlut.edu.cn
bNingbo Institute of Dalian University of Technology, Ningbo, 315000, China
First published on 17th November 2025
Mg3(Sb,Bi)2-based thermoelectric (TE) materials have attracted considerable attention due to their earth-abundant constituent elements, excellent environmental compatibility, and outstanding TE properties in the low- and medium-temperature range. However, optimizing TE properties remains challenging for Mg3(Sb,Bi)2, as does device structure design. Given the rapid development of Mg3(Sb,Bi)2-based materials and their devices in recent years, there is an urgent need for a systematic and critical review of the progress, challenges, and prospects. We systematically summarize and analyze various current strategies for improving TE properties, including component modulation, defect engineering, band engineering, grain boundary engineering, and mixed scattering optimization of ionized impurities and phonons. In addition, the key challenges for practical applications are comprehensively analyzed, particularly focusing on connection technology, barrier layer material screening, and device structure design. Finally, we outline the significant potential and current challenges of Mg3(Sb,Bi)2-based materials in low- and medium-temperature waste heat recovery and solid-state cooling applications, and provide future insights regarding TE properties modulation, preparation scale-up, and novel application prospects.
Mg3(Sb,Bi)2 is a typical Zintl compound that crystallizes in a CaAl2Si2-type structure, as shown in Fig. 1a. Its crystal structure features a three-dimensional chemical bonding network comprising Mg2+ cations in octahedral coordination and (Mg2Sb2)2− anionic layers in tetrahedral coordination, establishing a robust basis for excellent TE properties.17–21 As shown in Fig. 1b, the intrinsic band structure of Mg3(Sb,Bi)2 indicates that Mg3Sb2 is an indirect bandgap (Eg) semiconductor (0.4–0.8 eV) with relatively low carrier concentration (nH), whereas Mg3Bi2 is a type-II topological nodal-line semimetal characterized by the absence of Eg and high nH.22–25Fig. 1c shows the density of states (DOS) for Mg3Sb2 and Mg3Bi2. The valence band maximum (VBM) of Mg3Sb2 and Mg3Bi2 is dominated by the Sb/Bi 5p orbitals, while the conduction band minimum (CBM) is dominated by the Mg 3s orbitals.26,27 This band structure indicates that the electronic states of anions and cations dominate the electrical properties of p-type and n-type conduction samples, respectively. The Eg can be regulated by changing the Bi/Sb ratio through Mg3Sb2 and Mg3Bi2 alloying, enabling optimization of electrical transport properties.28–30 However, the volatility and high vapor pressure of Mg during synthesis result in a large number of Mg vacancies, significantly inhibiting electron transport properties.31,32 The breakthrough came when Tamaki et al.33 first synthesized n-type Mg3.2Sb1.5Bi0.49Te0.01 using excess Mg and Te-doping, achieving a remarkable zT of 1.5 at 723 K. This immediately attracted significant research interest. Subsequent studies demonstrated that interstitial Mn doping effectively suppresses the formation of Mg vacancies. This approach led to the first identification of a positive temperature-dependent Eg in Mg3+δSbxBi2−x, directly confirming the suppression of the bipolar effect at medium temperatures. Ultimately, Mg3+δSb1.5Bi0.49Te0.01Mn0.01 achieved an average zT (zTave) of 1.05 between 323 and 523 K. Notably, its favorable mechanical toughness highlights Mg3+δSbxBi2−x materials as promising candidates for TE applications near-room-temperature.34 Further demonstrating potential, a zT value of 0.9 was achieved at 350 K for the n-type Mg3.2Bi1.998−xSbxTe0.002 reported by Mao et al.,35 demonstrating its significance as a substitute for Bi2Te3 in low-grade waste heat recovery applications. Furthermore, in a cooling device with p-type Bi0.5Sb1.5Te3, it achieved a temperature difference (ΔT) of 91 K at a hot-side temperature (Th) of 350 K, exceeding the performance of commercial counterparts. A room-temperature analysis of the electronic quality factor (BET/κl) confirms that n-type Mg3Sb2-based alloys exhibit highly competitive TE properties compared to commercial Bi2Te3 alloys, solidifying their status as strong contenders for TE cooling and power generation near-room-temperature.36 Given this remarkable progress, a systematic and critical review is essential to summarize the advancements, challenges, and prospects of Mg3(Sb,Bi)2.
Currently, defect engineering,37 additive manufacturing,38 and device applications39,40 cover various aspects of TE materials. To facilitate the further development of Mg3(Sb,Bi)2, we systematically summarize the past advancements in TE performance regulation and device structure design. In this review, we commence by introducing the fundamental properties of Mg3(Sb,Bi)2, including crystal and electronic structures, to elucidate the intrinsic foundations of its transport behavior. Next, we systematically categorize strategies for enhancing properties, including component modulation, band engineering, grain boundary (GB) engineering, composite engineering, and mixed scattering optimization of ionized impurities and phonons. We then discuss connection techniques, the screening of barrier layer materials, and the design of device structures for TE modules of Mg3(Sb,Bi)2, addressing critical steps toward practical applications. Ultimately, this review provides valuable insights into the development of Mg3(Sb,Bi)2-based materials and promotes their application in low-grade heat recovery and solid-state cooling.
High-temperature melting/reaction is a traditional synthesis method valued for its simplicity and suitability for large-scale production.41 For instance, Y-doped ingots achieving grain sizes > 100 µm and room-temperature carrier mobility (µH) > 100 cm2 V−1 s−1 were prepared by high-temperature reactive melting with subsequent annealing.42 However, significant Mg volatilization during melting remains a challenge. This can be mitigated by adding excess Mg during sintering, which effectively compensates for Mg loss in conventional n-type Mg3(Bi, Sb)2 alloys.43 Recently, kilogram-level Mg3(Sb,Bi)2 powder was synthesized in merely 5 minutes via the self-propagating high-temperature synthesis method.44 This technique also achieves an oxygen content below 0.3%, offering a viable pathway for the industrial-scale production of powder.
As shown in Fig. 2a, solid-state synthesis in sealed tubes suppresses Mg volatilization and impurity formation through high-temperature reactions.45–47 Using this method, Ponnambalam et al.45 synthesized single-phase p-type Mg3SbBi, achieving a zT value of 0.4 at 825 K. The Mg3.2Bi1.4975Sb0.5Te0.0025 compound, prepared by sealing in tantalum tubes followed by melt annealing and directional solidification, exhibits record grain sizes larger than 500 µm and room-temperature µH of 220 cm2 V−1 s−1 among Mg3(Sb,Bi)2-based materials.48 Although offering low-cost sintering and minimal sintering equipment requirements, the high-temperature synthesis step necessitates specialized equipment and currently lacks scalability for industrial production.
![]() | ||
| Fig. 2 (a) Mg-flux method using a sealed tantalum tube.47 Copyright 2020, Royal Society of Chemistry. (b) Photograph of as-grown Mg3Bi1.49Sb0.5Te0.01 and Mg3.1Li0.003Bi1.49Sb0.5Te0.01 single crystals.51,52 Copyright 2020, Elsevier. Copyright 2024, American Chemical Society. (c) Synthesis process for mechanical alloying and sintering of Mg3(Sb,Bi)2.66,67 Copyright 2024, American Chemical Society. Copyright 2020, Tsinghua University Press. (d) Schematic diagram of the used 3D printer with a picture of the printed TE legs and their assembly into the device.63 Copyright 2025, The American Association for the Advancement of Science. | ||
Single crystals, with fewer grain boundaries, are expected to exhibit excellent near-room-temperature properties in principle, as GB effects can significantly limit carrier mobility µH and TE properties below 500 K.49,50 As shown in Fig. 2b, the Bridgman method enables the controlled growth of high-quality Mg3(Sb,Bi)2 single crystals with precise stoichiometry and low defect densities.51–53 For example, Mg3.1Li0.003Bi1.49Sb0.5Te0.01 achieved an excellent zT of 1.05 at 300 K by suppressing Mg inclusion during single crystal growth, highlighting their competitiveness for cooling applications.52 However, compared to polycrystalline synthesis, this method requires precise temperature control and slow growth rates, leading to higher energy consumption and limited productivity.
As illustrated in Fig. 2c, the combination of mechanical alloying (MA) with either hot pressing (HP) or discharge plasma sintering (SPS) offers a relatively simple and controllable approach for synthesizing Mg3(Sb,Bi)2.54–57 Single-phase n-type Mg3(Sb,Bi)2 were successfully synthesized using a multi-step ball milling process, achieving a remarkable peak zT value of 1.5 at 716 K.33 Particularly in Bi-rich compositions, the nanostructures induced by MA promote interface potential scattering of phonons, significantly suppressing the bipolar effect at high-temperatures.54 Nevertheless, like other specialized techniques, this method faces challenges for large-scale implementation due to high processing requirements.
Despite the many advantages of these preparation methods, developing novel Mg3(Sb,Bi)2-based materials with superior formability and excellent TE performance still faces many obstacles. Given the excellent plasticity of Mg alloys,58,59 deformation processing (e.g., stretching, rolling) offers a promising strategy to further optimize the microstructure and TE properties. Specifically, rolling refines grain boundaries, reduces porosity, and induces texture. These are critical for improving µH and mechanical strength. Furthermore, the enhanced plasticity and mechanical properties of Mg3Bi2 provide a solid foundation for exploring its potential in flexible TE device applications.60–62 Leveraging these properties enables future integration of Mg3Bi2 into flexible platforms, boosting η and device flexibility.
Beyond established methods, emerging additive manufacturing methods such as 3D printing (Fig. 2d) may offer significant potential for fabricating complex and gradient structures in future research, as evidenced in the Bi2(Sb,Te)3 system.63,64 This approach uniquely enables the precise construction of complex geometries and facilitates direct multi-material integration of TE materials, electrodes, and substrates, thereby advancing compact, fully integrated device development.65
![]() | ||
| Fig. 3 The relationship between the peak zT values of Mg3(Sb,Bi)2-based materials and the year of publication.35,51,52,68,74–98 | ||
In recent years, the rapid development of Mg3(Sb,Bi)2 has made it a promising alternative to Bi2Te3. Compared with Bi2Te3, Mg3(Sb,Bi)2 exhibits a lower cost due to its reduced need for the expensive Te element. Meanwhile, n-type Mg3(Sb,Bi)2 shows a good zT value of 0.9 at 300 K,73 making it more promising for solid-state cooling applications compared to Bi2Te3−xSex. Furthermore, Sb-rich compositions have good potential for power generation at medium temperatures, while Bi-rich compositions hold even more potential for cooling applications. The high cooling performance of Mg3Bi2-based alloys motivates future exploration of TE materials for cooling.
| σ = neµ | (1) |
![]() | ||
| Fig. 4 (a) The nH as a function of nominal composition Mg3+xSb1.5Bi0.5.105 (b) Illustration of the reduction in Eg due to alloying with Mg3Bi2.107 Copyright 2020, Elsevier. (c) Schematic nH dependence of zT for Mg3Sb2−xBix using the SPB model.109 Copyright 2020, The American Association for the Advancement of Science. (d) The µH and nH at 300 K as a function of the electronegativity difference |χQ − χMg| in n-type Mg3Sb1.5Bi0.5.74 Copyright 2018, John Wiley and Sons. (e) Doping atoms versus nH of Mg3(Bi,Sb)2-based compounds at room-temperature. (f) Temperature-dependent σ of Mg3(Sb,Bi)2-based materials.52,80,86,89,114,115 Band structures of (g) Mg24Sb16, (h) Pr1Mg23Sb15Se1, and (i) Nd1Mg23Sb15Se1.86 Copyright 2022, Elsevier. | ||
Chalcogen element doping (e.g., S, Se, Te) represents a common approach to enhance the properties of n-type Mg3(Sb,Bi)2. For example, S-doped Mg3Sb1.5Bi0.5 achieved a peak zT of 1.0 at 725 K, though its properties are inferior to Se- and Te-doped systems due to lower nH and µH, as shown in Fig. 4d.74 The decreasing electronegativity difference (|χQ − χMg|) from S to Te leads to enhanced nH and µH, ultimately improving TE properties. Subsequent studies revealed that Se doping offers environmental benefits, obtaining a peak zT of 1.24 at 498 K in Mg3.2Bi1.4Sb0.59Se0.01 and demonstrating its widespread prospect for low-temperature TE applications.87 Additionally, the combination of Te doping with defect control and GB optimization significantly enhances TE properties, underscoring the pivotal role of chalcogen doping in improving σ.112
In addition to anion doping with chalcogen elements, both rare-earth elements doping and single-crystal design enhance the σ of Mg3(Sb,Bi)2.56,78,81,113Fig. 4e demonstrates the higher nH achievable with rare-earth element doping, underscoring its potential for σ enhancement. Fig. 4f compares the temperature-dependent σ of single crystal and rare-earth elements doping.52,80,86,89,114,115 Elements like La and Ce act as promising dopants, enabling nH to approach the theoretical doping limit.89 In addition, the σ of Li–Te co-doped single crystals increased to 13 × 104 S m−1 from 6.2 × 104 S m−1 for Mg3.1Bi1.49Sb0.5Te0.01 at 300 K.52 To further understand the electron transport mechanism, the band structures of Mg24Sb16, Pr1Mg23Sb15Se1, and Nd1Mg23Sb15Se1 were calculated by first-principles calculation.86Fig. 4g shows that the CBM and VBM of pristine Mg3Sb2 are at different locations, which corresponds to indirect Eg semiconductors. As shown in Fig. 4h and i, Pr–Se and Nd–Se significantly increase the DOS near the EF, narrow the Eg, and enhance carrier transport, thereby improving the σ and PF. Ultimately, Mg3.17Pr0.03Sb1.5Bi0.49Se0.01 and Mg3.17Nd0.03Sb1.5Bi0.49Se0.01 achieved peak zT values of 1.67 and 1.74 at 700 K, which are 38% and 43% higher than that of Se-only doped sample, respectively. In conclusion, the optimized treatment of band engineering, doping, and defect control is crucial for improving nH and electrical transport properties in Mg3(Sb,Bi)2. The trade-off between nH optimization and bipolar effect suppression remains a key challenge, especially in high-temperature environments. Therefore, introducing deep-level dopants and leveraging energy filtering effects offer strategies to mitigate detrimental bipolar effects.
![]() | ||
| Fig. 5 (a) Temperature-dependent µH of Mg3.1A0.1Sb1.5Bi0.49Te0.01 (A = Fe, Co, Hf, and Ta).77 (b) Predicted upper-limit of the zT when grain boundary phase resistance is considered negligible (solid line).49 Copyright 2018, Royal Society of Chemistry. (c) SPS displacement as a function of sintering time. (d) Microstructural evolution of Mg3.2Bi1.5Sb0.498Te0.002Cu0.01.85 Copyright 2022, Springer Nature. (e) The EBSD orientation maps for Mg3.2Sb1.5Bi0.49Te0.01 samples sintered at 873 K and 1123 K.129 Copyright 2018, AIP Publishing. (f) and (g) The APT characterization of Nb0.1Mg3.05Sb1.99Te0.01.131 Copyright 2021, John Wiley and Sons. | ||
GB scattering severely limits room-temperature µH in Mg3(Sb,Bi)2.126,127 GB resistance reduces µH, leading to a decrease in the PF, restricting TE properties. Simulation of GB band offsets indicates that eliminating GB resistance could enhance the room-temperature zT of Mg3.2Sb1.5Bi0.49Te0.01 by over 60% (Fig. 5b), emphasizing the significance of GB engineering in enhancing zT.49 Optimizing sintering temperatures and incorporating a small amount of Cu have proven effective in GB engineering. As shown in Fig. 5c and d, optimal sintering produces dense microstructures with large grains and minimal porosity. Significant increase in grain size at 1073 K compared to 723 K sintering. The room-temperature peak zT of Mg3.2Bi1.5Sb0.48Te0.002Cu0.01 increased significantly to 0.9 due to the complete elimination of thermal GB resistance at low temperatures.85 Two-dimensional nanomaterials such as graphene have shown remarkable potential as GB additives. The introduction of graphene (G) can effectively regulate the space charge region at grain boundaries, reducing carrier transport barriers. For instance, in Mg3.24Sb1.5Bi0.49Te0.01/1.0 vol% G, the Eb decreased from 42 meV to 18 meV, enhancing room-temperature µH to 68 cm2 V−1 s−1.93 Additionally, graphene nanoplatelets (GNPs) have been reported to increase interface thermal resistance, reduce the temperature at grain boundaries, and optimize the interfacial S.128 Chemical doping and control of grain size effectively mitigate GB scattering. For example, indium doping significantly suppressed GB scattering and enhanced room-temperature µH.88 Moreover, the addition of transition metals such as Mo and Nb promotes grain growth via GB segregation, which reduces scattering and increases both µH and σ.68,95 For instance, Mo doping increased the average grain size of Mg3.16Mo0.04Sb1.5Bi0.49Te0.01 from 9.46 µm for Mg3.2Sb1.5Bi0.49Te0.01 to 18.07 µm.95 Similarly, introducing Nb impurities into Mg3Sb1.5Bi0.49Te0.01 improved the room-temperature zT by more than 200%, reaching an unprecedented high zT of 2.04 at 798 K.68 Grain size optimization also effectively suppresses GB scattering. As shown in Fig. 5e, the average grain size of Mg3.2Sb1.5Bi0.49Te0.01 increased from 1.0 µm to 7.8 µm,129 significantly reduced GB scattering, enhanced µH, and maintained stable κ.
The synergistic effects of grain size and texturing improved the σ of the layered Mg3.24Sb1.5Bi0.49Te0.01 structure, achieving a maximum zT of 1.71 and a high zTave of 1.12 perpendicular to the pressing direction (⊥P).94 Increasing the sintering temperature further enlarged grain sizes, reducing GB scattering and increasing the room-temperature σ by over 300%.130 Notably, the enhanced µH can be attributed to the weakening of GB scattering due to the grain growth and coarsening. Fig. 5f presents the scanning electron microscope (SEM) image of an atom probe tomography (APT) sample. Fig. 5g illustrates the three-dimensional distribution of Mg, Sb, Nb, and Te. These results indicate that the Nb fails to incorporate into the Mg3Sb2 matrix after doping but instead segregates along the GB in a metallic state. This segregation leads to the formation of a Nb-rich wetting layer, which promotes grain growth, therefore enhancing σ.131 Notably, the elimination of GB scattering in single crystals leads to superior TE properties.132 Demonstrating this advantage, Imasato et al.133 first synthesized n-type Te-doped Mg3Sb2 single crystals by employing a flux method followed by Mg vapor annealing. The elimination of GB resistance led to a remarkable improvement of over 100% in both µW and zT near-room-temperature.
While the aforementioned strategies have successfully enhanced the µH in Mg3(Sb,Bi)2, recent research has shifted toward gaining deeper, atomic-level insights into how defects and carriers interact locally. Fig. 6a–c show high-angle annular dark field (HAADF) and annular bright field (ABF) scanning transmission electron microscopy (STEM) images, along with the integrated differential phase contrast (iDPC) views of Mg3Sb2 single crystals. The iDPC imaging visualizes both Mg and Sb atoms in Mg3Sb2 single crystals.134 Furthermore, Fig. 6f illustrates the Fourier difference electron density maps (Fo − Fc) difference electron density map, confirming the presence of Mg cation vacancy defects and Frenkel defects.135 Meanwhile, the origin of the excellent TE properties is further explored by angle-resolved photoelectron spectroscopy (ARPES). Fig. 6g shows the constant energy diagrams of the ARPES of Mg3Bi2, Mg3Bi1.25Sb0.75, and Mg3Sb2 single crystals. The electronic states close to EF disappear from Mg3Bi2 to Mg3Sb2, and the valence band tops move significantly to higher binding energies, leading to a significant widening of the Eg.47 To elucidate the effect of doping on the local structure, synchrotron X-ray pair distribution function (X-PDF) was performed on Mg3.2(Bi0.7Sb0.3)1.99Te0.01, Mg3.2(Bi0.7Sb0.3)1.99Te0.01–Cu2Se 1.0%, and Mg3.2(Bi0.7Sb0.3)1.99Te0.01–Ag2Se 1.0%, as shown in Fig. 6f–h.136 In comparison to Mg3.2(Bi0.7Sb0.3)1.99Te0.01, the enhanced short-range ordering in the Cu2Se and Ag2Se composites, as evidenced by the narrowing of the characteristic peaks, provides a superior channel for carrier transport. These studies emphasize the crucial role played by GB engineering, defect control, and single-crystal methods in achieving superior µH. Since the cost-effective and controllable approach to achieve defect-free single crystals is still a major challenge, GB engineering combined with defect modulation strategies offers a preferred path to balance performance and cost.
![]() | ||
| Fig. 6 STEM images of Mg3Sb2: (a) HAADF, (b) ABF, and (c) iDPC images along [001] direction.134 Copyright 2021, Elsevier. (d) The Fo − Fc difference electron density map of n-type Mg3Sb2 before refining the interstitial position.135 Copyright 2023, American Chemical Society. (e) Series constant energy maps of Mg3(Sb,Bi)2.47 Copyright 2020, Royal Society of Chemistry Synchrotron X-PDF data as a function of r for (f) Mg3.2(Bi0.7Sb0.3)1.99Te0.01, (g) Mg3.2(Bi0.7Sb0.3)1.99Te0.01–Cu2Se 1.0 wt%, and (h) Mg3.2(Bi0.7Sb0.3)1.99Te0.01–Ag2Se 1.0 wt%.136 Copyright 2025, Elsevier. | ||
![]() | (2) |
![]() | ||
| Fig. 7 (a) Schematic diagrams of the band structures and DOS for band convergence and band anisotropy.138 Copyright 2021, Springer Nature. Electrical properties of (b) n-type and (c) p-type Mg3Sb2.109,138,140 Copyright 2020, The American Association for the Advancement of Science. Copyright 2021, Springer Nature. (d) Strain-dependent m* at conduction band edge.142 Copyright 2020, The American Association for the Advancement of Science. (e) Defect formation energies of the lanthanide substitution on the Mg1 site.157 Copyright 2020, John Wiley and Sons. (f) Temperature-dependent S of Mg3(Sb,Bi)2-based materials.52,80,86,89,114,115 (g) Schematic diagram of Eb near the GB and metal–semiconductor interface region.68 Copyright 2023, Springer Nature. | ||
10,29 have limited the commercial development of Mg3(Sb,Bi)2 TE devices. Strategies have primarily focused on optimizing the band structure and implementing effective doping to enhance nH.
As illustrated in Fig. 7d, under the influence of strain, the m* of the DOS at the CBM of Mg3Sb2 exhibited a linear decreasing trend.142 Theoretically, a reduction in m* leads to a lower S. However, when acoustic phonon scattering dominates, the maximum PF is proportional to the ratio of Nv/m*. Thus, a smaller m* corresponds to a higher µ and a higher PF.107 Meanwhile, biaxial strain can effectively enhance band degeneracy. For example, p-type Mg3Sb2 achieved a maximum S in the plane from 375 µV K−1 under unstrained to 623 µV K−1 under 2.5% biaxial strain at 700 K.143 Furthermore, in p-type Mg3(Sb,Bi)2, the introduction of elements such as Li,144 Zn,31 Cd,32 Ag,145 and Na146 at the cationic sites to replace Mg can effectively regulate the nH and µH. Meanwhile, multiple co-doping at the Mg sites simultaneously enhances both carrier transport and phonon scattering.32,147 For instance, a peak zT value of 0.87 was achieved in Mg1.95Na0.01Zn1Sb2 at 773 K through Na–Zn co-doping.98 First-principles calculations revealed that Ag–Bi co-doping shifts the EF into the VB and narrows the Eg. This facilitates hole generation, increasing nH from 3.50 × 1017 cm−3 to 7.88 × 1019 cm−3. As a result, Mg2.94Ag0.06Sb1.9Bi0.1 achieved an impressive PF of 778.9 µW m−1 K−2 at 723 K.92 Nanostructured bulk materials fabricated using high-energy ball milling combined with SPS exhibit uniform microstructures that enhance S.148
Band structure modulation via doping effectively enhances the S.154–156 The calculated formation energies for various dopants are summarized in Fig. 7e. La and Ce act as effective n-type dopants, exhibiting high donor transition energy levels significantly above the CBM.157Fig. 7f compares the S for single crystals and rare earth element doping.52,80,86,89,114,115 The S of all rare-earth doping increases with temperature, indicating degenerate semiconductor behavior. Moreover, the S of single-crystal Mg3.1Li0.003Bi1.49Sb0.5Te0.01 decreases with rising temperature owing to the presence of bipolar effect.52 The significant enhancement of S in Mn-doped Mg3Sb1.5Bi0.5 was attributed to the presence of multiple conduction bands with high band degeneracy.75 Furthermore, the addition of transition metals such as Nb/Ta can scatter low-energy carriers more efficiently by changing the interfacial barrier (Eb). As shown in Fig. 7g, the interfacial barrier decreases from Eb1 to Eb2, improving the tendency for more high-energy carriers to pass through, which leads to an increase in the S.68 In summary, achieving an optimal S not only requires precise control of nH but also demands a comprehensive strategy involving band engineering, defect chemistry, and interface optimization. Future research should further explore the interplay between band structure modulation and microstructural design to achieve an ideal balance between the S and σ.
![]() | (3) |
![]() | (4) |
![]() | ||
| Fig. 8 Neutron dynamical structure factor S(|Q|, ω) of (a) Mg3Sb2 and (b) Mg3Bi2.161 Copyright 2020, The American Association for the Advancement of Science. (c) Atom-projected phonon DOS of Mg3Sb2 and Mg3Bi2.160 Copyright 2023, American Chemical Society. (d) The Grüneisen parameter, and (e) scattering rate of anharmonic for 300 K of Mg3Sb2.160 Calculated phonon spectra of (f) Mg24Sb16, (g) YMg23Sb15Cl, (h) YMg23Sb14ClTe, and (i) the sound velocity comparison.162 Copyright 2025, Elsevier. | ||
In Mg3.2Bi2, partial substitution of Bi with Sb effectively reduces κl through alloy scattering. The κl reduced from 9.6 W m−1 K−1 for Mg3.2Bi2 to 3.2 W m−1 K−1 for Mg3.2Bi1.298Sb0.7Te0.002, leading to a high zT value of 0.9 at 350 K.35 Additionally, alloy scattering induced by compositional variations reduced the κl of Mg3.2Bi1.4Sb0.595Se0.01 to 0.6–0.7 W m−1 K−1 in the low-temperature range, demonstrating excellent room-temperature TE properties.87 Furthermore, the Sb–Bi lattice disorder introduced during alloying can further reduce κl to 0.44 W m−1 K−1 at 716 K, approaching the Cahill minimum κl.33 Point defects and lattice defects also significantly contribute to κl reduction by scattering phonons of different wavelengths. For instance, introducing heavy-mass Co enhanced phonon scattering, enabling Mg3.15Co0.05(Sb0.3Bi0.7)1.99Te0.01 to achieved a minimum κl of 0.64 W m−1 K−1 at 450 K and a high peak zT of 1.03 at 525 K.79 Furthermore, Mg3.18Y0.02Sb1.5Bi0.49Se0.01 leveraged point defects and dense grain boundaries to enhance phonon scattering, achieving a κl reduction. This effect yielded a theoretical η of 13.8% and an excellent zT of 1.87 under a 450 K temperature gradient.91 Doping with rare-earth elements such as Gd and Pr can further optimize nH and κ. By synergistic tuning the electronegativity difference, mass fluctuations, and point defects, Gd-doped Mg3.065Sb1.3Bi0.7Gd0.015 achieved a high zT of 1.55 at 700 K.82 Similarly, Pr-doping not only increased nH but also reduced κl to 0.42 W m−1 K−1, leading to a peak zT of 1.70 for Mg3.2Pr0.02Sb1.5Bi0.5 at 725 K.76
Nanostructured materials, through the introduction of interfaces and defects, can effectively scatter a wide range of phonons while maintaining favorable electronic transport properties.148,169,170 For instance, nanostructured Mg3Sb2 and Mg3Sb1.8Bi0.2 exhibited respective peak zT enhancements of 54% and 56% at 773 K, primarily due to preferential phonon scattering at interfaces.81,84Fig. 9a demonstrates a transmission electron microscopy (TEM) image of Mg2.94Ag0.06Sb1.9Bi0.1, indicating the presence of nano-precipitates and lattice distortion. Furthermore, lattice distortion was verified by high-resolution transmission electron microscopy (HRTEM) (Fig. 9b). The κl was effectively decreased to 0.48 W m−1 K−1 by phonon scattering at point defects, distorted lattices, and nanoprecipitates, as illustrated in the schematic diagram Fig. 9c.92 Similarly, Bi nanoprecipitation and lattice distortion significantly enhance phonon scattering. Combined with shorter relaxation times, this ultimately reduced κl for Mg3.18Ce0.02SbBi0.97Te0.03 to 0.51 W m−1 K−1 at 673 K.171 In Mg3.24Sb1.5Bi0.49Te0.01/1 vol% G, an equivalent nanoparticle scattering mechanism resulted in a moderate reduction in κl, achieving a peak zT value of 1.30 at 773 K.93 Furthermore, Mn doping suppresses Mg evaporation during the sintering process, helping to decouple the electrical and thermal transport. As shown in Fig. 9d and e, low-magnification TEM imaging reveals the grains fractured along the crystals. Fig. 9f presents an HRTEM image that allows the observation of nanoprecipitates and nanostructured GB grain boundaries, both contributing to lower κl.34 Dislocation and strain engineering have shown remarkable success in further optimizing phonon transport mechanisms. The aggregation and evolution of vacancies caused by Bi defects into dislocations enhanced phonon scattering, resulting in a peak zT of 1.82 for Mg3.16Y0.01Sb1.5Bi0.48 at 473 K.90 The presence of dislocations and strains in Mg3.2Sb1.5Bi0.49Te0.01–4% Mg and Mg3.2In0.005Sb1.5Bi0.49Te0.01 reduced the phonon relaxation time, leading to a remarkable decrease in κl.67,172 Phonon transport is further constrained by strain-induced disorder and lattice anharmonicity, which reduce phonon relaxation times. Additionally, as shown in Fig. 9g, dislocation strain and defect-induced mass fluctuations further suppress phonon transport. Mg3.22Ho0.03Sb1.5Bi0.5 sintered at 1073 K achieved a low κl of 0.47 W m−1 K−1 at 750 K.130
![]() | ||
Fig. 9 HRTEM images of (a) Mg2.94Ag0.06Sb1.9Bi0.1 and (b) rectangular area indicated in (a). (c) Schematic diagram of point defects strengthening phonon scattering.92 Copyright 2025, Elsevier. (d) and (e) Low-magnification TEM image of the sample with Sb content of Mg3.02Mn0.01Sb1.2Bi0.79Te0.01, with an inset depicting the electron diffraction patterns along the [1 1 ] direction. (f) HAADF image of the coarse-grained boundary.34 Copyright 2018, John Wiley and Sons. (g) Schematic diagram demonstrating the combined effect of strain engineering and scattering mechanisms on the modulation of phonon and charge carrier transport.130 (h) The lph dependence of accumulated κl and schematic diagram of all-scale hierarchical structures. Copyright 2024, Elsevier. | ||
Under the Debye model approximation, the lph can be expressed as eqn (5):
| lph = vsτ(x) | (5) |
Fig. 9h illustrates the relationship between the lph and the accumulated κl. Lattice defects of different scales scatter phonons of various frequencies and lead to extremely low κl. To effectively employ a multi-scale defects, the dimensions of the introduced scattering centers should be smaller than the κl and larger than the carrier mean free path.173 These studies reveal the complex mechanism of phonon transport in Mg3(Sb,Bi)2 materials and offer a theoretical basis for reducing κl by modulating phonon scattering. Future studies should integrate multiscale defect optimization with advanced characterization techniques to simultaneously balance electrical/thermal properties and advance TE performance limits.
![]() | (6) |
The manufacture of TE devices entails more intricate scientific and engineering issues than optimization of material properties. Appropriate interface and structural design strategies, including the connection technology, selection of barrier layers, multi-stage optimization, and size matching between n-type and p-type TE legs, are essential for enhancing device performance.150,175–178Fig. 10a illustrates a single-leg structure consisting of a Cu electrode, hot- and cold-side solders, a metal contact layer, a barrier layer, and a TE leg that can be used to test the performance of a single material. The traditional module (Fig. 10b) is suitable for η in terms of a combination of n- and p-type materials over a wide temperature range. However, optimization of leg size matching and interfacial stability between different materials is required. Though thermal expansion compatibility and interfacial contact resistance must be addressed, the segmented module (Fig. 10c) maximizes the η of the device. Additionally, issues such as thermal stability, mechanical properties, device assembly techniques, and manufacturing costs are key obstacles to address for large-scale applications.175,179–181Fig. 10d summarizes the η of Mg3(Sb,Bi)2-based TE devices at different ΔT. Mg3.2Sb1.5Bi0.49Te0.01 + 0.25 vol% TiO2−n single-leg achieved a η of 15% at ΔT of 470 K, which fully demonstrates its excellent potential for application.11,68,182–186
![]() | ||
| Fig. 10 Schematic diagram of the structure of (a) single-leg, (b) traditional module, and (c) segmented module. (d) The η of Mg3(Sb,Bi)2-based TE devices as a function of ΔT.11,68,182–188 | ||
![]() | ||
| Fig. 11 (a) Schematic of conventional soldering technique.197 Copyright 2023, Oxford University Press. (b) Mg3Sb2-based modules assembled by Ag composite paste under low-temperature and low-pressure conditions.197 (c) Mass change and heat flow of Ag composite paste during thermal cycling.197 (d) Analysis of the resistance changes observed before and after soldering with Ag composite paste.197 (e) Schematic diagram of TE device assembly by Ag nanoparticle sintered bonding.183 Copyright 2024, Springer Nature. (f) and (g) Characterization of the nano-Ag sintered connection layer.183 | ||
Compared with the traditional soldering technique, the novel composite Ag paste soldering technique effectively avoids the problem of device failure due to TE arm damage caused by pressure and temperature limitations. As shown in Fig. 11b, the Ag paste was uniformly applied to the electrode before heating and heated at 523 K in argon for 1.5 h. The connection between the TE legs and Cu electrode was realized after heating and exhibits good thermal stability (Fig. 11c).195,196Fig. 11d compares the resistance change before and after soldering of the new Ag composite paste to avoid impairing the TE material and interface, and the resistivity is consistent with the predicted values.197 This low-temperature, low-pressure connection technique utilizing Ag composite adhesive demonstrates effective performance. In addition, the efficiency of n-type single-leg Mg3Sb2-based devices prepared using a low-temperature nano-Ag sintering process is up to 13.3%.183 Nano-Ag sintering was performed at 573 K and 10 MPa, as shown in Fig. 11e. The Ag nanoparticles exhibited an almost spherical morphology and were encapsulated by an organic layer that inhibits oxidation and self-adhesion (Fig. 11f). The sintered nano-Ag network structure was homogeneously dispersed and well-connected at the interface (Fig. 11g). Also, this process can be combined with phase diagram engineering, which exhibits great potential for interconnect packaging of a wide range of high-performance TE power generation devices.
The Fe/Ni layers have been confirmed as proper barrier layers for Mg3(Sb,Bi)2.188,204,205 Taking the Fe contact layer as an example, the microstructure and compositional line scan of the Mg3.2Sb1.5Bi0.49Te0.01/Fe interface show that the mixed region is primarily enriched with Fe and Mg (Fig. 12a). Due to the minimal solubility of Fe and Mg, the excess Mg in the mixed layer is residual and destabilizing. Aging the Mg3.2Sb1.5Bi0.49Te0.01/Fe interface at 523 K for 400 h further elucidated its thermal chemical stability. As shown in Fig. 12b, the thickness of the Fe–Sb layer and the Sb-rich Mg3Sb2 layer increased to 1 µm and 0.7 µm, respectively. The overall thickness of the reaction layer exhibits a favorable linear relationship with the square root of time (Fig. 12c), indicating that the interface layer growth is governed by diffusion. Subsequent studies found that Mg2Cu was a suitable barrier layer for Mg3(Sb,Bi)2 compared to Ni/Fe electrodes.206 The ternary phase diagram of Cu–Mg–Sb(Bi) indicated that Mg2Cu is thermodynamically stable within the Mg–Mg2Cu–Mg3Bi1.5Sb0.5 phase region, as shown in Fig. 12d. Thus, it effectively prevents the diffusion of Cu into Mg3(Sb,Bi)2. The contact resistivity (ρc) of the Mg3.2Bi1.4975Sb0.5Te0.0025 legs, measured using the four-probe method, was as low as 12 mΩ cm2 (Fig. 12e).207Fig. 12f demonstrates the maximum thermal stress (70 MPa) in the contact layer of Cu/Mg2Cu/Mg3.2Bi1.4975Sb0.5Te0.0025 at a constant ΔT and optimized current, indicating that the material is used as an optimal buffer layer to relieve the thermal stresses.
![]() | ||
| Fig. 12 The composition profile of Mg3.2Sb1.5Bi0.49Te0.01/Fe interface after an aging treatment of (a) 0 h and (b) 400 h. (c) The thickness of the reaction layer as a function of aging time.206 Copyright 2020, Elsevier. (d) Mg–Mg2Cu–Mg3Bi1.5Sb0.5 three-phase region. (e) Contact resistance measurement of the interfacial layers. (f) Simulated thermal stress distribution near the barrier layers of the module during operation.207 Copyright 2022, Elsevier. (g) The energy dispersive spectroscopy mapping and line scanning results showing the elemental distribution of Ti foil sintered with Mg3(Sb,Bi)2. (h) Changes in power generation performance of a Mg3Sb2-based single-leg during thermal cycling testing.208 Copyright 2024, Springer Nature. | ||
Through the dynamical transition from a high-temperature substable phase to a low-temperature stable phase, the limitations of thermodynamic equilibrium can be overcome, and ideal materials for TE barrier layers are identified. By comparing the kinetic properties of the reaction and diffusion behavior of Mg3(Sb,Bi)2 and Ti during sintering and operation, the Ti foil/Mg3(Sb,Bi)2 connectors present excellent sintering activity and inert behavior.208 As shown in Fig. 12g, the interface reaction product between the Ti foil and Mg3(Sb,Bi)2 was a MgTiSb compound with a thickness of 1 µm when held at 1023 K for 15 min. Therefore, it can function effectively as a barrier layer for Mg3Sb2-based materials. As illustrated in Fig. 12h, the power generation performance of a Mg3Sb2-based single-leg TE device remained stable during 200 thermal cycles lasting 330 h, with the heat source temperature cycling between 473 K and 748 K, exhibiting excellent stability and reliability. The interface ρc reduced to below 5 µΩ cm2, allowing the Mg3Sb2-based module to achieve a high η of 11% at a ΔT of 440 K, surpassing most advanced medium-temperature TE modules. Advances in novel connection techniques and the development of stabilizing barrier layers have greatly improved the reliability, efficiency, and performance of Mg3(Sb,Bi)2-based TE devices. These strategies address critical issues related to interfacial loss, contact resistance, and thermal stress.
000 thermal cycles.218 These studies further demonstrate the reliability and high η values of the Te-free modules, achieving competitive performance in the low temperature range.
![]() | ||
| Fig. 13 (a) Single-leg η measurement architecture under high ΔT, with operational setup visualization.11 Copyright 2019, Elsevier. (b) Comparison of the maximum η among the single-leg devices in 0.1Nb/Mg3Sb1.5Bi0.49Te0.01 with other TE modules.68 Copyright 2023, Springer Nature. (c) Module fabrication and cooling measurement setup for n-type Mg3.2Bi1.5Sb0.498Te0.002Cu0.01 and p-type α-Mg0.99Cu0.01Ag0.97Sb0.99.85 Copyright 2022, Springer Nature. (d) Cooling performance of p-type MgAg0.97Sb0.99 and n-type Mg3.6Y0.003Sb0.6Bi1.4.218 Copyright 2022, Royal Society of Chemistry. (e) Design of the setup used to characterize the cooling performance of the Mg3Bi2-based double-stage TEC.9 Copyright 2025, Innovation Press. (f) and (g) Solid-state sensor device and applications of Mg3.19Co0.01Sb0.5Bi1.49Te0.01.228 Copyright 2024, Elsevier. | ||
The CTE mismatch between the two TE legs can lead to stress accumulation during repeated thermal cycles, ultimately degrading performance.219,220 For medium- and high-temperature applications, minimizing the CTE mismatch is an important strategy for improving device performance. For example, p-type MgAgSb and n-type Mg3.49Y0.01Sb1.5Bi0.5 exhibited a CTE mismatch of only 0.95%, demonstrating good compatibility between the two leg materials within TE modules.221 Furthermore, Peltier modules mediated by ultrafine grains and nanoporous structures in α-MgAgSb and Mg3.2(Bi,Sb)2 achieved a maximum ΔT of 52 K and a coefficient of performance (COP) of 8.3 under high-temperature conditions.222 For medium- and high-temperature modules, using n- and p-type legs from the same material system can reduce CTE mismatches and simplify electrode material selection.223 For example, n-type Mg3.2Sb1.5Bi0.49Te0.01Mn0.01 and p-type Na0.005Mg2.995Sb1.5Bi0.5 exhibited comparable thermal expansion properties, achieving a theoretical η of approximately 9.5% under a ΔT of 300 K to 773 K.224
Besides power generation, TE cooling is equally important.225–227 Full-scale TE cooling devices with 7 pairs, 31 pairs, and 71 pairs were fabricated by joining an n-type Mg3.2Bi1.4975Sb0.5Te0.0025 layer with an effective and reliable Mg2Cu diffusion barrier to a p-type (Sb0.75Bi0.25)2(Te0.97Se0.03)3.207 Nevertheless, despite the progress in high-performance single-pair and full-scale TE coolers, the development of TE devices entirely free of (Bi, Sb)2Te3-based TE devices remains limited. Inspired by this, a cooling module was constructed by combining n-type Mg3.2Bi1.5Sb0.5 with p-type α-MgAgSb through microstructure design (Fig. 13c). It achieved an η value of 2.8% (ΔT = 95 K), a maximum cooling ΔT of 56.5 K, and a maximum COP (COPmax) of 2.6 (ΔT = 10 K).85 Furthermore, a TE device composed of high-performance n-type Mg3.6Y0.003Sb0.6Bi1.4 and p-type MgAg0.97Sb0.99 compounds exhibits a maximum ΔT of 72 K with Th of 347 K (Fig. 13d).218 A 7-pair TE cooling module fabricated using p-type Bi0.5Sb1.5Te3 and n-type MgAg0.97Sb0.99 reached a maximum ΔT of 100 K at Th = 400 K, exhibiting COPmax of 8.9, 7.5, and 3.2 at Tc of 0 K, 5 K, and 10 K, respectively.48 The subsequently developed multistage refrigeration device follows the design of heat distribution where the cooling capacity of the next stage is larger than the heat dissipation of the previous stage for realizing the cooling stage by stage.9 The temperature field simulation results show that the seven-stage refrigeration device of Mg3.1Sb0.497Bi1.5Te0.003/(Bi,Sb)2Te3 can realize step-by-step refrigeration by the layer-by-layer stacked “pyramid” configuration. As shown in Fig. 13e, the Mg3.1Sb0.497Bi1.5Te0.003/(Bi,Sb)2Te3 two-stage refrigeration device achieved a maximum cooling capacity of 3.0 W and a maximum ΔT of 103.2 K, which is comparable to the commercial Bi2Te3 alloy-based device.
The potential interdisciplinary applications of Mg3(Sb,Bi)2 have also been explored. Fig. 13f shows an infrared thermogram of a respiratory mask equipped with Mg3.19Co0.01Sb0.5Bi1.49Te0.01-based TE sensor, demonstrating its potential for respiratory monitoring applications. Additionally, Fig. 13g reveals that five common liquids have different contact angles on the polyurethane film, indicating the TE sensor can be used for liquid identification.228 These advancements suggest that Mg3(Sb,Bi)2-based modules have the potential to rival, or even surpass, commercial Bi2Te3 systems in performance, offering an environmentally friendly and efficient alternative that represents a significant step forward in the sustainable development of TE technology.
![]() | ||
| Fig. 14 Challenges and optimization aspects of Mg3(Sb,Bi)2.181,186,215 Copyright 2023, John Wiley and Sons. Copyright 2023, Oxford University Press. Copyright 2024, Springer Nature. Copyright 2023, Elsevier. | ||
(1) In terms of fabrication, methods such as mechanical alloying, discharge plasma sintering, and hot-pressing have been widely employed for producing Mg3(Sb,Bi)2 bulk materials. These techniques not only improve material densification but also effectively enhance lattice uniformity and TE properties. Furthermore, advancements in single-crystal fabrication technologies have drastically reduced electrical resistance caused by grain boundary scattering, further improving carrier mobility (µH) and η. Current fabrication techniques are mostly limited to laboratory-scale production, which cannot meet the demands of large-scale industrial manufacturing. Developing cost-effective and environmentally friendly synthesis methods is essential for the commercialization of these materials.
(2) By adjusting the Bi/Sb ratio, the bandgap of Mg3(Sb,Bi)2 can be flexibly tuned to optimize carrier concentration and µH, thereby enhancing electrical transport properties. Additionally, doping with rare-earth elements has further improved the Seebeck coefficient (S) and power factor, with some studies reporting a peak zT exceeding 2.0 across a wide temperature range from room temperature to medium- and high-temperatures. Achieving the ideal balance between electrical conductivity with S and minimizing undesirable effects such as bipolar effects remains a key challenge in realizing the full potential of practical TE applications.
(3) The introduction of multiscale defects significantly reduces lattice thermal conductivity (κl) by effectively suppressing thermal conduction through phonon scattering, resulting in outstanding TE properties in the low- and medium-temperature range. Critically, while enhanced phonon scattering reduces κl, charge carrier mobility µH must be considered to maintain excellent electrical transport properties. Future research should screen Mg3(Sb,Bi)2-based compositions with intrinsically low κl through combined first-principles calculations and high-throughput methods, while exploring emerging phonon transport mechanisms like liquid-like ionic diffusion and topological phonon states.
(4) Mg3(Sb,Bi)2-based TE devices exhibit significant prospects for both power generation and cooling applications. To narrow the gap between theoretical and actual η, researchers have systematically developed interface barrier layer design and connection process optimization strategies. While Fe/Ni, Mg2Cu, and Ti barrier materials exhibit excellent contact resistance, chemical stability, and thermal expansion coefficient matching, long-term reliability under thermal cycling remains challenging due to interfacial reactions. Novel connection processes like low-temperature Ag sintering and composite Ag paste soldering prevent degradation from temperature-pressure coupling while enabling robust bonding, yet face scalability challenges requiring improved uniformity and long-term thermal cycling reliability.
(5) Enhancing the long-term stability of TE modules and optimizing device structural design are central to the practical utilization of TE devices. By optimizing p-type MgAgSb and n-type Mg3(Sb,Bi)2 with matched thermal expansion coefficients, the accumulation of thermal stresses can be significantly reduced. The combination of geometrical optimization and surface coating technology improves the long-term stability of the modules. In addition, emerging manufacturing technologies such as 3D printing are expected to overcome the bottleneck of mass production. With its excellent performance and scalable applications, Mg3(Sb,Bi)2 is expected to be a core candidate for the next generation of TE materials, which will support the advancement of clean energy technology.
| This journal is © The Royal Society of Chemistry 2026 |