Improved figure of merit (z) at low temperatures for superior thermoelectric cooling in Mg3(Bi,Sb)2

The low-temperature thermoelectric performance of Bi-rich n-type Mg3(Bi,Sb)2 was limited by the electron transport scattering at grain boundaries, while removing grain boundaries and bulk crystal growth of Mg-based Zintl phases are challenging due to the volatilities of elemental reactants and their severe corrosions to crucibles at elevated temperatures. Herein, for the first time, we reported a facile growth of coarse-grained Mg3Bi2-xSbx crystals with an average grain size of ~800 μm, leading to a high carrier mobility of 210 cm2 · V−1 · s−1 and a high z of 2.9 × 10−3 K−1 at 300 K. A \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta$$\end{document}ΔT of 68 K at Th of 300 K, and a power generation efficiency of 5.8% below 450 K have been demonstrated for Mg3Bi1.5Sb0.5- and Mg3Bi1.25Sb0.75-based thermoelectric modules, respectively, which represent the cutting-edge advances in the near-room temperature thermoelectrics. In addition, the developed grain growth approach can be potentially extended to broad Zintl phases and other Mg-based alloys and compounds.

The low-temperature thermoelectric performance of Bi-rich n-type Mg 3 (Bi,Sb) 2 was limited by the electron transport scattering at grain boundaries, while removing grain boundaries and bulk crystal growth of Mg-based Zintl phases are challenging due to the volatilities of elemental reactants and their severe corrosions to crucibles at elevated temperatures.Herein, for the first time, we reported a facile growth of coarse-grained Mg 3 Bi 2-x Sb x crystals with an average grain size of ~800 μm, leading to a high carrier mobility of 210 cm 2 • V −1 • s −1 and a high z of 2.9 × 10 −3 K −1 at 300 K.A ΔT of 68 K at T h of 300 K, and a power generation efficiency of 5.8% below 450 K have been demonstrated for Mg 3 Bi 1.5 Sb 0.5 -and Mg 3 Bi 1.25 Sb 0.75 -based thermoelectric modules, respectively, which represent the cutting-edge advances in the near-room temperature thermoelectrics.In addition, the developed grain growth approach can be potentially extended to broad Zintl phases and other Mgbased alloys and compounds.
Thermoelectrics (TEs) can directly convert heat into electricity, or vice versa, and enable broad applications in solid-state cooling and power generation [1][2][3] .Generally, the performance of thermoelectric device is governed by the materials' figure of merit z, which can be expressed as z = S 2 σ= κ L + κ e À Á , where S, σ, κ L , and κ e are the Seebeck coefficient, electrical conductivity, lattice thermal conductivity, and electronic thermal conductivity, respectively 4 .Here, the Seebeck coefficient (S) and electrical conductivity (σ) are coupled through the carrier concentration (n) and the band effective mass (m * ) 3 .Among the strategies to decouple S and σ for enhanced power factor (PF = S 2 σ), various strategies have been employed, including the resonant level doping 5 and band convergence 6 .However, for material systems with intrinsically low lattice thermal conductivity, such as SnSe(S) 7 , Cu 2 Se(S,Te) 8 , MgAgSb 9 , etc., the grain growth has been proposed as an effective strategy to enhance the carrier mobility and figure of merit z, through reduced grainboundary scattering and having nearly no impacts on Seebeck coefficient.
The high thermoelectric performance of n-type Mg 3 (Bi,Sb) 2 originated from its large degeneracy number of the conduction band up to six 10,11 and the intrinsically low lattice thermal conductivity of 0.497 W • m −1 • K −1 for Mg 3 Sb 2 and 0.349 W • m −1 • K −1 for Mg 3 Bi 2 10 .Among these, Bi-rich alloys possess relatively higher z values around room temperatures (150-550 K), given the reason that the small effective mass of narrow bandgap material favors high carrier mobility at low temperatures 12,13 .It is revealed that the carrier scattering mechanism of n-type Mg 3 (Bi,Sb) 2 was dominated by grain boundaries [14][15][16][17][18] and several methods on single crystal growth have been reported 16,[19][20][21] .However, unfortunately, the above methods involve complicated apparatus 19 and multiple experimental steps 16,21 , and the impurities introduced by high-temperature processes can severely degrade TE properties.Moreover, the irregular shapes and limited single crystal sizes (lamellar shape with a thickness of 0.2-1.2mm) add extra difficulties in thermoelectric device fabrications.Therefore, the facile growth of high-quality Mg 3 (Bi,Sb) 2 crystals remains to be challenging.
Grain growth usually performs at high temperatures, posing difficulties in deliberate control of the content of volatile reactants and meanwhile the raw materials could react with the crucible, leading to deviation from the stoichiometry of the target composition.To solve this problem, we hereby developed a universal approach to synthesize a variety of coarse-grained alloys or compounds containing active elements, e.g.Mg and Sb, represented by the Zintl phase Mg 3 (Bi,Sb) 2 .Following this method, we have prepared bulk polycrystal n-type Bi-rich Mg 3 Bi 2-x Sb x (x = 0.5, 0.75) with grain size up to 1.0 mm.Such grain size encompasses almost twothird of the length of a thermoelectric leg in the devices.The chemical composition, defects and grain size of Mg 3 Bi 2-x Sb x materials were successfully regulated to approach the limit of carrier mobility for a single crystal.Combined with the intrinsically low lattice thermal conductivity 22,23 , the z values of both Mg 3 Bi 1.5 Sb 0.5 and Mg 3 Bi 1.25 Sb 0.75 increased significantly below and at room temperatures, which is beneficial for the cooling purpose.Moreover, the asgrown coarse-grained bulk crystals (with a dimension of +12:7mm × 13mm) also exhibit excellent uniformity and mechanical performance, bringing great advantages for subsequent thermoelectric device fabrication in this work.This technique is also applicable for the crystal growth of MgAgSb 9 , Zintl phase 24 , and other Mg-based alloys 25 , which contain reactive and volatile elements (such as alkaline metals, and rare earth elements).

Growth of coarse-grained bulk crystal and microstructure characterizations
To better control the synthesis of Mg 3 (Bi,Sb) 2 crystals, a ZrO 2 /Ta/ SiO 2 three-layer nested-crucible (Supplementary Fig. 1) was developed to seal and load the highly reactive and corrosive elements, Mg and Sb (details available in Methods).A temperature procedure for the synthesis and growth of bulk Mg 3 Bi 2-x Sb x (x = 0.5, 0.75) was designed based on the pseudo-binary phase diagram of Mg-Bi 0.75 Sb 0.25 26,27 , as shown in Fig. 1a.The crystal growth process involved the following steps: The samples were placed in a rocking furnace and initially heated to 1273 K to melt the precursors.The molten sample was continuously swung for 2 h to ensure a complete reaction and the formation of a uniform liquid phase.Subsequently, the samples were slowly cooled down to 1193 K at a rate of 2 K/h and held at this temperature for 24 h.At 1193 K, the high-temperature Mg 3 (Bi,Sb) 2 phase began to solidify from the melt state.Nucleation is likely initiated at specific sites in the crucible, such as the junction of the sidewall and the bottom of the crucible.The 24-hour holding period facilitated the preservation of high-quality nuclei for subsequent grain growth, while low-quality nuclei gradually shrank and melted away.Next, the crucible was further cooled down to 1073 K at a rate of 2 K/h and annealed for 24 h.At this stage, the amount of the high-temperature phase (cubic phase) gradually increased by consuming the liquid phase during cooling.The small amount of remaining liquid phase between the grains promoted grain growth during the annealing process, as depicted in the inset of Fig. 1a.Following this, the crucible was cooled down to 973 K at a rate of 2 K/h and held at this temperature for 72 h.Throughout this step, the high-temperature phase completely transformed into the uniform hexagonal phase (room temperature phase) of Mg 3 Bi 2-x Sb x (x = 0.5, 0.75).The longtime annealing helped eliminate defects in the material.Finally, the crucible was allowed to naturally cool down to room temperature.The obtained coarse-grained Mg 3.15 Bi 1.4975 Sb 0.5 Te 0.0025 and Mg 3.15 Bi 1.2475 Sb 0.75 Te 0.0025 samples are referred to as CGBi1.5 and CGBi1.25 in this context, respectively.The measured densities of as-grown CGBi1.5 and CGBi1.25 materials are 5.36 g/cm 3 and 5.14 g/cm 3 , respectively, agreeing well with the experimental and theoretical values 17,19 (Supplementary Table 1).
EBSD was employed to characterize the grain morphology of the as-grown CGBi1.5 bulk material.Contrast with the small grain sizes of  19,30,31 .The inset shows the percentage of improvements in z compared with the data reported by Mao et al. 31 d Comparison of the maximum ΔT of real cooling devices, including our CGBi1.5/Bi0.5 Sb 1.5 Te 3 cooling module and new material-based modules 30,[35][36][37] .Inset shows the as-fabricated modules.~20 μm observed in the spark plasma sintered (SPS) sample (Supplementary Fig. 2), the obtained bulk ingot has an average grain size of ~800 μm, and with some regions reaching ~1.0 mm as shown in Fig. 1b.Here the large grain sizes in CGBi1.5 bulk material exhibited a typical feature of the coarse grain 28 .Considering the dimensions (~1.3 × 1.3 × 3 mm 3 ) of thermoelectric legs in our cooling device, ideally the transport of electrons in one single leg can be scattered by only two grain boundaries, which is analogy to the scenarios in a single crystal.Moreover, the as-grown CGBi1.5 bulk materials show a Vickers hardness of 0.54 GPa, close to that of 0.61 GPa for Mg 3.065 Sb 1.3 Bi 0.7 Gd 0.015 as reported in the literature 29 .The compressive strength of the CGBi1.5 sample is greater than 300 MPa, because no cracks occur under the loading stress of 300 MPa, which is consistent with the literature results 29 , as shown in Supplementary Fig. 3.The above results affirm the excellent mechanical properties of our coarsegrained materials, fulfilling the needs for TE device manufacturing.
The obtained bulk material was cut and polished for the measurement of thermoelectric properties.Figure 1c and Supplementary Fig. 4 show that the CGBi1.5 sample exhibited a peak z of 2.92 × 10 −3 K −1 at 290 K, which was the highest among literatures for Mg 3 Bi 1.5 Sb 0.5 materials to the best of our knowledge 19,30,31 .Most of the thermoelectric cooling devices work in the temperature range of 150-300 K, where the coarse-grained CGBi1.5 demonstrated remarkable advantage over other Mg 3 (Bi,Sb) 2 based materials, e.g., the average z is 31.8%higher than the data reported by Mao et al. 31 , indicating a great potential for thermoelectric cooling applications [32][33][34] , as shown in the inset of Fig. 1c.A 7-pair thermoelectric cooling module was fabricated based on above CGBi1.5 and the commercial p-type Bi 0.5 Sb 1.5 Te 3 .A cooling temperature difference (ΔT) of 68 K was obtained at the hotside temperature (T h ) of 300 K, which is the largest among the recently reported full-scale new material-based real cooling devices such as Mg 3 (Bi,Sb) 2 /Bi 2 Te 3 30 , Mg 3 (Bi,Sb) 2 /MgAgSb 35,36 , and SnSe/Bi 2 Te 3 37 (Fig. 1d).It is noted that the temperature difference was underestimated by ~2 K due to the use of alumina ceramic plates.
The Rietveld refinement against experimental powder X-ray diffraction (XRD) data identified the formation of phase-pure hexagonal Mg 3.15 Bi 1.4975 Sb 0.5 Te 0.0025 (space group P 3m1, No. 164, a = 4.6652(1) Å, c = 7.4058(2) Å; Fig. 2a).The Rietveld refinement also suggested a site occupancy factor (SOF) of Mg1 of 92 ± 1%, i.e., corresponding to the stoichiometry of Mg 2.92(1) Bi 1.4975 Sb 0.5 Te 0.0025 , which is in the similar range and slightly higher than that of 89 ± 2% as reported for Mg 3.2 Sb 1.5 Bi 0.49 Te 0.01 52 , revealing the possible enhanced n-type carrier transportation in our sample.The phase purity of the as-grown CGBi1.25 sample was also confirmed by the XRD as shown in Supplementary Fig. 5.To confirm the composition uniformity of the material, EDS analyses were performed for a total of nine regions from a 12.7 diameter CGBi1.5 pellet.Figure 2b shows the probability distribution of the atomic ratios across these regions, which confirms the uniform elemental distributions across the whole area and is consistent with the elemental mapping results as shown in Supplementary Fig. 6.The high-resolution scanning transmission electron microscopy (STEM) image in Fig. 2c showed the regular arrangement of atoms in the

Thermoelectric property characterization
As shown in Fig. 3a, the increasing of electrical resistivity with temperature shows metallic behavior of CGBi1.5 and CGBi1.25 samples with a respective carrier concentration of 2.08 × 10 19 cm −3 and 2.2 × 10 19 cm −3 .Owing to the reduced grain-boundary scattering, the electrical resistivity of CGBi1.5 samples at room temperature was measured to be 14.3 μΩ m, decreased by ~25% compared to that of SPS-processed polycrystalline Mg 3.2 Bi 1.4975 Sb 0.5 Te 0.0025 samples (19.2 μΩ • m) with a similar carrier concentration 30 .The literature data of grain sizes (d) and corresponding carrier mobilities (μ H ) of Bi-rich ntype Mg 3 (Bi,Sb) 2 materials 12,19,36 were collected and compared with the measured data in this work, as shown in Fig. 3b.A grain-boundary dominated transportation model was adopted to unveil the μ H ∼ d relation (details available in Methods).The experiment data in this work was in good agreement with the fitting curve, and the height of the potential barrier of 170 meV was obtained for the grain-boundary scattering.As the grain size increases, the effect of grain-boundary on the mobility decreases rapidly, and the CGBi1.5 coarse grain with an average grain size of ~800 μm exhibits a high carrier mobility (210 cm 2 • V −1 • s −1 ) approaching the limit for a single crystal 19 , indicating the high quality of the crystallization of the coarse-grained bulk materials.
Given the difference in density of states' effective mass (m * d ) with varied Bi/Sb ratio, the weighted mobility (μ w ) can better describe the  16, 19-21, 30, 31, 53-55 .d Temperature-dependent Seebeck coefficient with the calculated bandgap in comparison to the literatures 12,13 shown in the inset.e Temperature-dependent thermal conductivity, lattice and bipolar thermal conductivity.f z values compared with state-of-the-art n-type Bi 2 Te 2.7 Se 0.3 materials 40,41 .inherent carrier transport properties 38 .As shown in Fig. 3c, the weighted mobility of coarse-grained CGBi1.5 crystals (276 cm 2 • V −1 • s −1 ) was the highest among the n-type Mg 3 (Bi,Sb) 2 materials.Despite the similar μ H , the μ w of CGBi1.5 sample was higher than that of the single crystal with the same Bi/Sb ratio.Since the grain size nearly has no impact on the effective mass of conduction band, the Seebeck coefficient of coarse-grained CGBi1.5 crystals was very close to that of polycrystals with the same carrier concentration 30 .As shown in the inset of Fig. 3d, the Goldsmid-Sharp bandgap E g was calculated from the maximum Seebeck coefficient and the temperature at which it occurs 39 .The Seebeck coefficient of CGBi1.5 and CGBi1.25 reached their peak values at 400 K and 525 K, and the bandgaps are 0.19 eV and 0.28 eV, respectively, which was consistent with polycrystalline materials with the same Bi/Sb ratio 12,13 , further confirming the excellent uniformity and high quality of coarse-grained crystals.Owing to a wider bandgap, the Seebeck coefficient of CGBi1.25 was less affected at elevated temperature.The enhanced electrical conductivity together with the high Seebeck coefficient enabled the high power factor of 33 μW • cm −1 • K −2 for CGBi1.5, and 27.5 μW • cm −1 • K −2 for CGBi1.25 at 300 K, respectively, as shown in Supplementary Fig. 7.
Figure 3e shows the combination of lattice and bipolar thermal conductivity κ L + κ b À Á as the function of temperatures, and the value of 0.80 W • m −1 • K −1 for CGBi1.5 at 300 K is comparable to that of reported values, i.e., 0.796 W • m −1 • K −1 for Mg 3 Bi 1.49 Sb 0.5 Te 0.01 single crystal 19 and 0.79 W • m −1 • K −1 for Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 polycrystal 31 .The results indicated that the grain-boundary scattering can pose significant impact on the carrier transportation through the potential barrier located at the grain-boundary, while having less effect on the phonons.This mainly because the phonon is less coupled with the electric field and its mean free path is much smaller than the grain size.
The improved electrical transport and intrinsically low lattice thermal conductivity significantly improved the z to 2.9 × 10 −3 K −1 for CGBi1.5 at 300 K, as can be seen from Fig. 3f, which is comparable to state-of-the-art commercial n-type Bi 2 Te 2.7 Se 0.3 materials (z ~2.94 × 10 −3 K −1 at 300 K) and outperformed most of the values reported in the literatures 40,41 .Thermoelectric properties of commercial Bi 2 Te 3 materials were shown in Supplementary Fig. 8.As for CGBi1.25, the peak z reaches ~2.7 × 10 −3 K −1 at the temperature of 350 K.The electrical transport properties of CGBi1.5 samples in different parts and from different samples have been tested (Supplementary Fig. 9 and 10), attesting to the uniform thermoelectric performance achieved through our method.Supplementary Fig. 11 reveals the good reproducibility of this synthesis method represented by the thermoelectric properties of several batches of CGBi1.5 samples.The z of CGBi1.5 in 150-300 K is comparable to that of commercial n-type Bi 2 Te 2.7 Se 0.3 , indicating its potentially high cooling performance.

Performance of thermoelectric modules
For CGBi1.5/Bi 0.5 Sb 1.5 Te 3 -based module, the measured ΔT max reached 68 K, 74 K, and 82 K with the hot-side temperatures of 300 K, 325 K, and 350 K, respectively, as shown in Fig. 4a and Supplementary Fig. 12, which are the highest among all the cooling devices with measurable cooling power density based on Mg 3 Bi 1.5 Sb 0.5 30,35   , Mg 3 Bi 1.4 Sb 0.6 36 and SnSe 37 materials (Supplementary Fig. 13), and also comparable to the state-of-the-art commercial Bi 2 Te 3 cooling module.The difference in ΔT max curve between the modules reported in this work and commercial Bi 2 Te 3 -based devices can be attributed to the higher nearroom temperature thermoelectric properties of Bi 2 Te 3 and the advanced contact layer fabrication technique in the leading thermoelectric enterprise.It is worth noting that Mao et al. reported a similar ΔT of 69 K at T h = 300 K with a unicouple TE cooling device based on , SnSe 37 and Mg 3 (Bi,Sb) 2 materials 35,36,[45][46][47] .f Aging time dependent maximum cooling temperature differences of CGBi1.5/Bi0.5 Sb 1.5 Te 3 devices and electrical resistivity of CGBi1.5 samples.Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and Bi 0.5 Sb 1.5 Te 3 , while the relevant cooling capacity (Q c ) and coefficient of performance (COP) were hard to obtain because of the large uncertainty in the measurement of small heat flow.As for CGBi1.25/Bi0.5 Sb 1.5 Te 3 module, despite the relatively low ΔT of 62 K at T h = 300 K, the higher z values of CGBi1.25 at 350-400 K lead to larger ΔT at high temperatures, i.e., 85 K at T h of 350 K, and 100 K at T h of 400 K, as shown in Fig. 4a and S12.Therefore, these two kinds of Mg 3 Bi 2-x Sb x /Bi 0.5 Sb 1.5 Te 3 (x = 0.5, 0.75) TE modules could substitute commercial Bi 2 Te 3 modules over a wide temperature range.
As shown in Fig. 4b, the thermoelectric module created a larger q c at a smaller ΔT between the two sides of the module because less cooling heat was dissipated through thermal conduct in the materials.The heat flow increased gradually with the working current and reached its maximum cooling power density (q c,max ).With the increase of temperature difference from 0 K to 10 K, the q c,max decreased from 1.07 W cm −2 to 0.94 W cm −2 at the optimal current of around 2.3 A. In this work, the normalized q c,max of 1.0 W cm −2 at ΔT = 5 K is 25% higher than that of the cooling module made of SPS-processed Mg 3 Bi 1.5 Sb 0.5 and commercial Bi 0.5 Sb 1.5 Te 3 materials, which is originated from the higher power factor of CGBi1.5 samples (Supplementary Fig. 7).COP, another important parameter of cooling module, defined as the ratio of Q c to the input power p, is critical to evaluate the energy conversion efficiency.The measured maximum COP of CGBi1.5/Bi0.5 Sb 1.5 Te 3 module reached 8.9, 7.5, and 3.2, at the temperature difference of 0 K, 5 K, and 10 K, respectively, as shown in Fig. 4c.These values were higher than the Mg 3 (Bi,Sb) 2 /MgAgSb module 35 (COP max = 2.6 at ΔT of 10 K) and comparable to the state-of-the-art commercial Bi 2 Te 3 module (COP max = 8.8 at ΔT of 0 K) (as shown in Supplementary Fig. 14).
It is known that industrial waste heat can be classified into six categories based on the sources, as illustrated in Fig. 4d.Among them, waste heat below 473 K occupies ~31% of total waste heat.As for the power generations of CGBi1.25/Bi0.5 Sb 1.5 Te 3 modules, a conversion efficiency as high as 5.8% was demonstrated within the temperature range of 280-450 K, as shown in Fig. 4e.At the ΔT of 100 K, both Mg 3 Bi 2-x Sb x /Bi 0.5 Sb 1.5 Te 3 (x = 0.5, 0.75) modules exhibit similarly high conversion efficiencies, which is originated from their comparable average z values in the temperature range of 280-390 K, as indicated in Fig. 4e and Supplementary Fig. 4. The efficiencies at a variety of ΔT from 20 to 170 K was the highest among those of reported modules based on Bi 2 Te 3 system [42][43][44] , Mg 3 (Bi,Sb) 2 /MgAgSb 35,36,45 , SnSe/Bi 2 Te 3 37 and so on 46,47 .
Additionally, a ΔT of 25 K in this module could produce an output power of 4 mW with a conversion efficiency >1%, which can use the tiny temperature difference between human skin and the environment to support the operation of wearable electronic devices 48 .With the waste heat source of 450 K, an output power density up to 0.14 W • cm −2 can be produced by this high-efficiency CGBi1.25/Bi0.5 Sb 1.5 Te 3 based module (η~5.8%).The open-circuit voltage (V oc ) and internal resistance (R in ) as a function of T h of Mg 3 Bi 2-x Sb x /Bi 0.5 Sb 1.5 Te 3 (x = 0.5, 0.75) modules are shown in Supplementary Fig. 15 and 16.The measured V oc is consistent well with the prediction from the TE properties.The discrepancy between the experimental and predicted values of R in was ascribed to the contact resistance (13.5 μΩ • cm 2 for Mg 2 Cu/CGBi1.5 interfaces, as shown in Supplementary Fig. 17).Therefore, high thermoelectric performance of the as-grown bulk Bi-rich Mg 3 Bi 2-x Sb x (x = 0.5, 0.75) materials was confirmed, revealing great potentials for both cooling and power generation at near-room temperatures.
As shown in Fig. 4f, the maximum temperature differences of CGBi1.5/Bi0.5 Sb 1.5 Te 3 devices were only reduced by ca. 1 K after being stored in the glove box for nearly 6 months.In addition, we also observed that with reduced Mg1 vacancies, the electrical resistivity of coarse-grained CGBi1.5 bulk materials stored in the glove box remained unchanged over 16 months, which, to the best of our knowledge, has been the best record to date.These reinforce the excellent stabilities of both the thermoelectric materials and electrodes of Mg 3 Bi 2-x Sb x /Bi 0.5 Sb 1.5 Te 3 (x = 0.5, 0.75) modules.The results also provide feasibility to further improve device-level stabilities through subsequent device packaging strategies, underpinning potential realworld applications.
In conclusion, two coarse-grained Mg 3 Bi 2-x Sb x (x = 0.5, 0.75) bulk crystals were synthesized and grown through a high-temperature melting approach.The as-grown CGBi1.5 crystals exhibit a high carrier mobility and low lattice thermal conductivity, leading to a high z of ~2.9 × 10 −3 K −1 at room temperature.Meanwhile, this material also exhibits the highest average z value of ~2.57× 10 −3 K −1 in the temperature range of 150-300 K among n-type Mg 3 (Bi,Sb) 2 materials with the same Bi/Sb ratio.The 7-pair cooling module made of n-type CGBi1.5 and p-type commercial Bi 0.5 Sb 1.5 Te 3 yielded a high ΔT of 68 K with the hot-side temperature of 300 K.Besides that, the CGBi1.25/Bi0.5 Sb 1.5 Te 3 devices displayed a high power-generation efficiency of 5.8% within the temperature range of 280-450 K.This work realized the growth of coarse-grain Mg 3 Bi 2-x Sb x bulk crystals through the nested-crucible, and facilitate its applications in high-performance thermoelectric cooling and power-generation modules.Moreover, the grain growth approach for enhancing thermoelectric properties, as developed in this study, shows promise for broader applications in Zintl phases, along with other Mg-based alloys and compounds.

Materials synthesis
High purity Mg shots, Bi shots, Sb shots, and Te pieces were weighted according to the composition of Mg 3.15 Bi 1.4975 Sb 0.5 Te 0.0025 and Mg 3.15 Bi 1.2475 Sb 0.75 Te 0.0025 .Raw materials with a total mass of 12 g were loaded into ZrO 2 crucibles in a glove box under Ar atmosphere with O 2 and water level below 0.1 ppm.The crucibles were embed into Ta tubes and the Ta tubes were sealed through arc-melting under Ar atmosphere.Then they were put into quartz tubes and sealed in vacuum.ZrO 2 crucible can protect Ta tube from the reaction of melting Sb element, and the sealed Ta tube prevent the reaction of Mg vaper with the quartz.The crystal growth follows the procedures, within a rocking furnace the samples were firstly heated up to 1273 K.The molten sample was swung continuously for 2 h in the rocking furnace.Then the samples were slowly cooled down to 1193 K at a rate of 2 K/h and kept for 24 h.Next, the crucible was cooled down to 1073 K at a rate of 2 K/h and annealed for 24 h.In succession, the crucible was cooled down to 973 K at a rate of 2 K/h and maintained for 72 h.Finally, the crucible was cooled down to room temperature naturally.

Phase and microstructure characterization
Powder XRD was performed using a Bruker D2 Phaser diffractometer with the Cu Kα radiation (λ = 1.5418Å, voltage of 40 kV, and emission current of 40 mA) in Bragg-Brentano geometry.The XRD patterns were collected at room temperature over a 2 theta range of 10-80 o with a step size of 0.006°for 40 min.Rietveld refinement against ascollected XRD data for CGBi1.5 was performed using GSAS via the EXPGUI interface 49 .Considering the very small amount of Te doping, a published Mg 3 Bi 2 structure 27 was modified by accommodating 25% of Sb and 75% of Bi at the 2d Bi site, i.e. with the fixed respective SOF of 0.25 and 0.75; the structure was employed as an initial crystallographic model.A constraint of the same coordinates and thermal displacement factors was then applied to Bi and Sb atoms.During the refinement, isotropic atom thermal displacements were applied to the heavy Sb and Bi atoms, while anisotropic types were set for the relatively light Mg atoms.SOFs of Mg2 and Mg1 were carefully refined, whereas the refinement of the former led to a SOF(Mg2) unchanged from unity, indicating the unfavored conditions for Mg2 vacancies within the anionic layers.However, a SOF(Mg1) lower than 1.0 was consistently suggested, and the stoichiometry of the sample was finally evaluated.
Results of the refinement are available in Supplementary Tables 2 and 3 in the SI file and Fig. 2a in the main manuscript.The microstructures, grain size, and element distribution were conducted by scanning electron microscopy (SEM, S8100, Hitachi) equipped with an energy dispersive spectrometer (EDS) and electron back-scattering diffraction (EBSD) detector.STEM was conducted using JEOL ARM200F.

Measurement of thermoelectric properties
For high-temperature (>300 K) measurement, both electrical and thermal transport properties were measured at the same direction.The obtained ingots were cut into bars with a dimension of 2.5 × 2.5 × 10 mm 3 for the simultaneous measurements of Seebeck coefficient and electrical resistivity (Linseis, LSR-3, Germany).Disks with diameter of 12.7 mm and thickness of 1.5 mm were cut from the ingots for the thermal diffusivity measurement.The thermal conductivity (κ) was calculated by κ = ρC p D, where ρ was density measured by Archimedes drainage method, C p was the specific heat estimated according to the Dulong-Petit law, and D was the thermal diffusivity measured by LFA 1000 (Linseis, Germany).The measurement uncertainties are 3%, 5%, and 6% for electrical resistivity, Seebeck coefficient, and thermal diffusivity, respectively, which yield an error in zT of ~20%.
For low temperature (150-370 K) measurement, all the transport properties were measured at one bar with a dimension of 3 × 1.5 × 7 mm 3 .The thermal conductivity and Seebeck coefficient were measured simultaneously (Supplementary Fig. 18) using TTMS (Thermal Transport Measurement System, MultiFields Technology), and the electrical resistivity was measured by PPMS (Physical Properties Measurement System, Quantum Design).The bars were polished to 3 × 0.5 × 7 mm 3 for the Hall coefficient (R H ) measurement under a reversible magnetic field (±2 T) using PPMS.The Hall carrier concentration ðn H Þ and the Hall carrier mobility (μ H ) were calculated by n H = 1=eR H and μ H = σ= en H À Á , respectively, where e is the electronic charge and σ is the electrical conductivity.

Mechanical property measurement
Vickers hardness was tested at room temperature with an applied load of 10 gf maintained for 10 s on a microhardness tester (MH500, China).The compressive strength was performed on samples with dimensions of 2 × 2 × 4 mm 3 using a universal test machine (Instron 5967, American).

Module preparation and cooling performance measurement
The sintering process of the materials, interfacial layers and contact layers was shown in our previous work 30 .The test method of contact resistivity was also the same as reported before 30 , and the measured contact resistivity in this work was 13.5 μΩ cm 2 (Supplementary Fig. 17).The obtained sandwich pellets were diced into squares with a dimension of 2.0 × 2.0 × 3 mm 3 using a dicing machine (DS830, Heyan technology).These legs were polished to 1.3 × 1.3 × 3 mm 3 to eliminate the influence of surface corrosion layer and then annealed at 613 K for 2 h.The p-type Bi 2 Te 3 materials were cut into squares with an optimized area of 1.25 × 1.25 mm 2 .The fabrication of thermoelectric modules and the cooling performance measurement process can also be found in previous report 30 .The power-generation performance was characterized using Mini-PEM (ADVANCE RIKO, Japan).The hot-side temperature T h of the module was controlled by a heater, and the cold-side temperature T c was controlled by the flowing water.Both T h and T c were measured by thermocouples of the Mini-PEM.The output power (P) and cold-side heat flow (Q c ) were recorded using the Mini-PEM.Therefore, the power-generation efficiency was calculated by η = Calculation of weighted mobility.The single parabolic band (SPB) model is expressed as follows: 50 In the above equations, k B is the Boltzmann constant, e the electron charge, h the Plank constant, r the scattering factor, m * d the density of state effective mass, η the reduced Fermi level, x the reduced carrier energy, and F n η ð Þ the Fermi-Dirac integral.The weighted mobility μ w is calculated by 38 where m e is the electron mass.
The fitting of grain size and carrier mobility.According to Matthiessen's rule, the total mobility μ can be depicted as: given the carriers are scattered by grain boundary, ionized impurity, acoustic phonons, alloy disorder and so on.Considering that grain size only affects grain-boundary scattering, other terms in Eq. 6 except μ À1 GB can be approximated as constants C. Expression of μ GB is given by: 51 where d is the grain size, m * the band mass, and E b the height of the potential barrier.The curve of μ H ∼ d was fitted according to Eqs. 6 and 7.The literature data used for the fitting can be found in Supplementary Table 4.
Calculation of cooling performance of thermoelectric devices.The cooling capacity (Q c ) of a thermoelectric cooling module can be expressed as: where S np , T h , T c , I, R, R c , K are the total Seebeck coefficient of device, the hot-side temperature of device, the cold-side temperature of device, the electric current, the total electrical resistance of thermoelectric materials on the legs, the total contact resistance from the interface layer, and the thermal conductance through the thermoelectric legs, respectively.S np , R, R c and K can be expressed as: Article https://doi.org/10.1038/s41467-023-40648-5 Nature Communications | (2023) 14:4932 where N, S n , S p , h, A n , A p , ρ n , ρ p , ρ nc , ρ pc , κ n , κ p are the number of the thermocouples, the Seebeck coefficient of n-type and p-type materials, the height of the thermoelectric legs, the cross-section area of n-type and p-type legs, the electrical resistivity of n-type and p-type materials, the contact resistivity of n-type and p-type legs, the thermal conductivity of n-type and p-type thermoelectric materials, respectively.The input power (p) on the device can be expressed as: The COP is defined as the ratio of Q c to p: The cooling power density (q c ) is given by: where S ceramic is the area of ceramic plate.

Fig. 1 |
Fig. 1 | The synthesis and growth mechanism of coarse-grained Mg 3 Bi 2-x Sb x (x = 0.5, 0.75) bulk material and characterization of the thermoelectric performances for CGBi1.5 and the relevant TE modules.a Schematic of the pseudobinary Mg-Bi 0.75 Sb 0.25 phase diagram, the synthesis procedures, and the setup of the crucibles.The picture of as-grown coarse-grained Mg 3 Bi 2-x Sb x bulk crystals was shown on the right.b The electron back-scattering diffraction (EBSD) image of as-grown coarse-grained crystals, and the photograph of a single CGBi1.5 leg.c The measured z of as-grown CGBi1.5 sample as a function of temperature, compared to literatures19,30,31 .The inset shows the percentage of improvements in z compared with the data reported by Mao et al.31 d Comparison of the maximum ΔT of real cooling devices, including our CGBi1.5/Bi0.5 Sb 1.5 Te 3 cooling module and new material-based modules30,[35][36][37] .Inset shows the as-fabricated modules.

Fig. 2 |
Fig. 2 | Crystal and microstructure characterizations of CGBi1.5. a Rietveld refinement against the X-ray diffraction patterns at 300 K for CGBi1.5.b The probability distribution of elemental compositions at different regions of a 12.7 mm diameter CGBi1.5 pellet.c The high-resolution scanning transmission electron microscopy (STEM) image of CGBi1.5, with the fast Fourier transform (FFT) diffractogram shown in the inset.d The STEM image of CGBi1.5 and line profile along the dashed blue line marked in the image.

Fig. 4 |
Fig. 4 | Cooling and power-generation performances of Mg 3 Bi 2-x Sb x / Bi 0.5 Sb 1.5 Te 3 (x = 0.5, 0.75) based TE modules.a The comparison of maximum cooling temperature differences between Mg 3 Bi 2-x Sb x /Bi 0.5 Sb 1.5 Te 3 (x = 0.5, 0.75) TE modules and commercial Bi 2 Te 3 -based module with a variety of hot-side temperatures T h .b The cooling power density (q c ) as a function of current for CGBi1.5/Bi 0.5 Sb 1.5 Te 3 based modules, in comparison to the 7-pair SPS Mg 3 Bi 1.5 Sb 0.5 / Bi 0.5 Sb 1.5 Te 3 cooling module with the normalized cross-section.c COP for CGBi1.5/Bi 0.5 Sb 1.5 Te 3 -based modules at temperature differences of 0 K, 5 K, and 10 K, respectively.All the dashed dot lines represent the theoretical prediction and the open circles indicate experimental data in b, c. d The distribution of industrial waste heat, where the waste heat at T < 473 K (mostly cooling medium and waste water and steam) takes up ~31% 56 .e The measured conversion efficiencies as a function of temperature difference for the Mg 3 Bi 2-x Sb x /Bi 0.5 Sb 1.5 Te 3 (x = 0.5, 0.75) TE modules compared to that of commercial Bi 2 Te 3 and results from literatures based on Bi 2 Te 3 42-44