Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
N-TYPE MG3.2BI2-BASED MATERIALS FOR THERMOELECTRIC COOLING APPLICATION
Document Type and Number:
WIPO Patent Application WO/2021/007126
Kind Code:
A1
Abstract:
New thermoelectric materials, such as Mg3Bi2-based Zintl phase compounds are described, where the semi-metallic Mg3.2Bi2 show an unexpectedly large Seebeck coefficient at 350 K and enhanced thermoelectric performances.

Inventors:
REN ZHIFENG (US)
MAO JUN (US)
Application Number:
PCT/US2020/040818
Publication Date:
January 14, 2021
Filing Date:
July 03, 2020
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
UNIV HOUSTON SYSTEM (US)
International Classes:
H01L35/14; H01L35/16; H01L35/18
Domestic Patent References:
WO2017143213A12017-08-24
Foreign References:
US20170250334A12017-08-31
US10121953B22018-11-06
US20100307551A12010-12-09
US20090211619A12009-08-27
Attorney, Agent or Firm:
MAAG, Gregory L. et al. (US)
Download PDF:
Claims:
CLAIMS

What is claimed is:

1. A thermoelectric material, comprising:

a Mg3Bi2-based Zintl phase compound, wherein said Mg3 2Bi2 Zintl phase compound comprises a ZT of about 0.9 at about 350 K.

2. The thermoelectric material of claim 1 , wherein said Mg3Bi2-based Zintl phase compound comprises a Seebeck coefficient of about -100 mV K 1 at about 350 K.

3. The thermoelectric material of claim 1 , wherein said compound comprises a electron-to-hole weighted mobility ratio greater than 8.

4. The thermoelectric material of claim 1 wherein said compound comprises an enhanced thermoelectric performance.

5. The thermoelectric material of claim 1 , wherein said Mg3Bi2-based Zintl phase compound comprises a partial substitution of Bi, wherein said Mg3 2Bi2-based Zintl phase compound further comprises at least one of Sb, or Te.

6. The thermoelectric material of claim 5, wherein said compound is Mg3 2Bii.498Sbo.5Teo.oo2·

7. The thermoelectric material of claim 1 , wherein the Mg3Bi2-based Zintl phase compound is a n-type Mg32Bi2-based Zintl phase compound.

8. The thermoelectric material of claim 6, wherein said Mg3Bi2-based Zintl phase compound comprises a unicouple of n-type Mg3.2Bii.498Sbo.5Teo.oo2 and p-type Bio.5Sb1.5Te3.

9. The thermoelectric material of claim 8, wherein said unicouple comprises a temperature difference of ~91 K and a hot-side temperature of ~350 K.

10. The thermoelectric material of claim 9, wherein said unicouple comprises a thermoelectric legs, and wherein said thermoelectric legs comprise an enhanced cooling performance.

11. The thermoelectric material of claim 1 , further comprising a contact material.

12. The thermoelectric material of claim 11 , wherein said contact material is at least one of Fe, Ni, or Cu.

13. The thermoelectric material of claim 12, wherein said contact material is Fe.

14. The thermoelectric material of claim 10, wherein said thermoelectric legs comprise Fe/Mg3.2Bii.498Sbo.5Teo.oo2/Fe.

15. The thermoelectric material of claim 15, wherein the thermoelectric material is a pure phase material.

16. The thermoelectric material of claim 5, wherein a power factor of the thermoelectric material is between about 20 pW cm 1 K2 and about 25 pW cm 1 K2 over a temperature range of about 200 K to 350 K.

17. The thermoelectric material of claim 5, wherein an electrical resistivity of the thermoelectric material is from more than 0.58pQ*m to less than 9 pQ*m from about 2 K to about 300K.

18. The thermoelectric material of claim 5, wherein a thermal conductivity of the thermoelectric material is about 2.8 W m 1 K 1 at about 300 K.

19. The thermoelectric material of claim 6, wherein said Mg3Bi2-based Zintl phase compound comprises a Seebeck coefficient of about -107 mV K 1 at about 350 K.

20. The thermoelectric material of claim 1 , wherein said compound has a narrow band gap of between about 0.16 and about 0.26 eV at room temperature.

Description:
N-TYPE MG 3.2 BI 2 -BASED MATERIALS FOR THERMOELECTRIC

COOLING APPLICATION

CROSS-REFERENCE TO RELATED APPLICATIONS

[0001]The present application claims benefit of U.S. Provisional Patent Application No. 62/870,884, filed July 5, 2019, and entitled “N-Type MG 3 2 Bi 2 -Based Materials for Thermoelectric Cooling Application,” which is hereby incorporated herein by reference in its entirety for all purposes.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

[0002] N/A

FIELD OF THE INVENTION

[0003] The disclosure herein is related to the field of n-type Mg 3.2 Bi2-based compounds with improved thermoelectric performance, wherein n-type Mg 32 Bi 2 -based materials provide an improved thermoelectric cooling.

BACKGROUND

[0004] The goal of utilizing waste heat for electricity via the Seebeck effect and cooling via the Peltier effect have motivated research on thermoelectrics for decades. When applying a thermoelectric material at a certain temperature, the electron-hole pairs could be thermally activated (i.e., bipolar conduction), greatly deteriorating the thermoelectric performance. Hence, embodiments herein address such needs and demonstrate that the electron-to-hole weighted mobility ratio {A), which can be used to quantify the asymmetry of transport properties between the conduction and valence bands, plays a pivotal role in bipolar conduction. When A has an extreme value, either A » 1 or A « 1 , the effect of bipolar conduction can be substantially suppressed. The disclosure herein addresses such needs, and discloses that it has been found that the n-type semi- metallic Mg 32 Bi 2 (nominal composition) with a large A above eight exhibits an unexpectedly large Seebeck coefficient above -100 mV K 1 at 300 K and very promising thermoelectric performance. The optimized Mg 3.2 Bi 298 Sb 0.7 Te 0.002 (nominal composition) can achieve a peak ZT above 0.9 at 350 K. The n-type Mg 3.2 Bi 2 -based materials show great applicability for thermoelectric cooling. When paired with p-type Bio .5 Sb 1.5 Te 3 , a large temperature difference of ~91 K can be realized at the hot-side temperature of ~350 K, which outperforms the commercial thermoelectric legs consisting of p-type Bio .5 Sb 1.5 Te 3 and n-type Bi 2 Te 2.7 Se 0 . 3.

BRIEF SUMMARY

[0005] The disclosure herein relates in some embodiments to a thermoelectric material, comprising a Mg 3 Bi 2 -based Zintl phase compound, wherein the Mg 3.2 Bi 2 Zintl phase compound comprises a ZT of about 0.9 at about 350 K. In some embodiments, of the thermoelectric material, the Mg 3 Bi 2 -based Zintl phase compound comprises a Seebeck coefficient of about -100 mV K 1 at about 350 K. In other embodiments of the thermoelectric material, the Mg 3 Bi 2 -based Zintl phase compound comprises a electron- to-hole weighted mobility ratio greater than 8. In some further embodiments, the Mg 3 Bi 2 - based Zintl phase compound comprises an enhanced thermoelectric performance.

[0006] In another embodiment of the thermoelectric material disclosed herein, the Mg 3 Bi 2 -based Zintl phase compound comprises a partial substitution of Bi, wherein the Mg 3.2 Bi 2 -based Zintl phase compound further comprises at least one of Sb, or Te. In another embodiment of the thermoelectric material the compound is Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and in a further embodiment of the thermoelectric material the Mg 3 Bi 2 -based Zintl phase compound is a n-type Mg 3.2 Bi 2 -based Zintl phase compound, and in a still further embodiment of the thermoelectric material the Mg 3 Bi 2 -based Zintl phase compound comprises a unicouple of n-type Mg 3.2 Bi 1.498 Sbo .5 Teo .002 and p-type Bio .5 Sb 1.5 Te 3 . In some embodiments of the thermoelectric material, the unicouple comprises a temperature difference of ~91 K and a hot-side temperature of ~350 K.

[0007] In other embodiments, the unicouple comprises a thermoelectric legs, and wherein the thermoelectric legs comprise an enhanced cooling performance, and in some embodiments the thermoelectric material comprising a further contact material.

[0008] In other embodiments of the thermoelectric material, the contact material is at least one of Fe, Ni, or Cu; in another embodiment the contact material is Fe; and in a further embodiment thermoelectric legs comprise Fe/Mg 3.2 Bi-i . Sbo . Teo . oo /Fe. In some embodiments of the thermoelectric material, the thermoelectric material is a pure phase material. In other embodiments, the power factor of the thermoelectric material is between about 20 pW cm 1 K 2 and about 25 pW cm 1 K 2 over a temperature range of about 200 K to 350 K; electrical resistivity of the thermoelectric material is from more than 0.58rW * th to less than 9 rW * th from about 2 K to about 300K; and in some embodiments the thermal conductivity of the thermoelectric material is about 2.8 W m 1 K 1 at about 300 K. In another embodiment of the thermoelectric material, the Mg 3 Bi 2 -based Zintl phase compound comprises a Seebeck coefficient of about -107 mV K 1 at about 350 K; and in a further embodiment of the thermoelectric material, the compound has a narrow band gap of between about 0.16 and about 0.26 eV at room temperature.

BRIEF DESCRIPTION OF THE DRAWINGS

[0009] For a detailed description of the disclosed embodiments, reference will now be made to the accompanying drawings, wherein:

[0010] Figure 1 (a) depicts the electron concentration and electron-to-hole weighted mobility ratio dependent Seebeck coefficient of a n-type semi-metal with a zero bandgap at 300 K;

[0011] Figure 1 (b) depicts the electron concentration and electron-to-hole weighted mobility ratio dependent Seebeck coefficient of a n-type semiconductor with a bandgap of 0.08 eV at 300 K;

[0012] Figure 1 (c) depicts the electron concentration and electron-to-hole weighted mobility ratio dependent Seebeck coefficient of a n-type semiconductor with a bandgap of 0.3 eV at 300 K;

[0013] Figure 2. Depicts (a) the carrier concentration (the upper panel) and mobility (the lower panel) of Mg 3 2 Bi 2 . (b) Electrical resistivity, (c) Seebeck coefficient, and (d) power factor of Mg 3.2 Bi 2. xTe x (x = 0, 0.002, and 0.004). The Seebeck coefficient of polycrystalline Bi is plotted in (c) for comparison;

[0014] Figure 3. Depicts (a) the calculated band structure of Mg 3 Bi 2 with the band overlap shifted to -0.1 eV. (b) Partial Seebeck coefficients of electrons and holes, (c) electron and hole mobilities, and (d) electron-to-hole weighted mobility ratio for Mg 3 2 Bi 2 ; [0015] Figure 4. Depicts the electronic thermoelectric properties of Mg 3.2 Bi 1.998- x Sb x Te 0.002 · (a) Electrical resistivity, (b) Seebeck coefficient, (c) power factor, (d) carrier concentration, and (e) mobility (f) Estimated bandgap of undoped Mg 3.2 Bi 2-x Sb x ;

[0016] Figure 5. Depicts (a) the lattice thermal conductivity and (b) ZT of Mg 3.2 BBii 2 1.998- x Sb x Te 0.002 · (c) Experimental setup for the thermoelectric cooling measurement with a unicouple consisting of p-type Bi 0.5 Sb 1.5 Te 3 and n-type Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 · (d) Electrical current dependent temperature difference between the hot and cold sides at the hot-side temperature of ~350 K. The inset shows the comparison of hot-side temperature dependent maximum temperature difference (DT) between this work and the commercial thermoelectric (TE) legs;

[0017] Figure 6. Depicts the hermoelectric properties of Mg 3+5 Bi 2 ( <5 = 0, 0.05, 0.1 , and 0.2). (a) Electrical resistivity, (b) Seebeck coefficient, (c) power factor, (d) thermal conductivity, (e) lattice thermal conductivity, and (f) ZT wherein the Mg 3 Bi 2 shows strong p-type conduction, while Mg 3 Bi 2 with excess Mg in the nominal composition of Mg 3+5 Bi 2 ( <5 = 0.5, 0.1 , and 0.2) is n-type (Fig. S1b). The Seebeck coefficient of p-type Mg 3 Bi 2 is only ~50 mV K ' 1 at 350 K but it is as high as ~-100 mV K '1 for n- type Mg 3+d Bi 2 . Such a noticeable difference in Seebeck coefficient can be partially ascribed to the large electron-to-hole weighted mobility ratio. As a result, the power factor of p-type Mg 3 Bi 2 is much lower compared to that of the n-type Mg 3+5 Bi 2 . In addition, it is very interesting to note that all of the n-type Mg 3+d Bi 2 samples show high peak thermal conductivity of more than 20 W m ' 1 K ' 1 around 20 K and it is only ~3.8 W m '1 K '1 for p-type Mg 3 Bi 2 . Such a difference in thermal conductivity should be ascribed to the strong phonon scattering by Mg vacancies in p-type Mg 3 Bi 2 ;

[0018] Figure 7. Depicts the thermoelectric properties of polycrystalline bismuth (a) Electrical resistivity, (b) Seebeck coefficient, (c) thermal conductivity, and (d) ZT;

[0019] Figure 8. (a) depicts the carrier concentration and (b) mobility of Mg 3.2 Bi 2-x Te x (x = 0, 0.002, and 0.004);

[0020] Figure 9. (a) depicts the thermal conductivity and (b) ZT of Mg 3.2 Bi 2-x Te x (x = 0, 0.002, and 0.004);

[0021] Figure 10. (a) depicts calculated band structure, (b) the band structure with the band overlap energy shifted to -0.1 eV;

[0022] Figure 11. depicts the calculated band overlap of Mg 3 2 Bi 2 from the two-band modelling;

[0023] Figure 12 depicts the calculated Fermi surfaces of (a) conduction and (b) valence bands of Mg 3 Bi 2 ; [0024] Figure 13 depicts the temperature-dependent electrical resistivity of Mg 3.2 Bi 2-x Sb x (x = 0, 0.1 , 0.2, 0.3, 0.4, 0.5, 0.6, and 0.7);

[0025] Figure 14 depicts the calculated average power factor of Mg3 Bh . 998- x Sb x Te0.002 between 2 K and 350 K;

[0026] Figure 15 depicts the thermal conductivity of Mg3 Bh . 998- x Sb x Te0 . 002;

[0027] Figure 16 depicts the thermoelectric properties of Bio . sSb^Tes. (a) Electrical resistivity,

(b) Seebeck coefficient, (c) thermal conductivity, and (d) figure of merit; and

[0028] Figure 17 depicts the contact resistance measurement of Fe/Mgs . ^h ^SbosTeo^ / Fe.

NOTATION AND NOMENCLATURE

[0029] Certain terms are used throughout the following description and claims to refer to particular system components. As one skilled in the art will appreciate, computer companies may refer to a component by different names. This document does not intend to distinguish between components that differ in name but not function. In the following discussion and in the claims, the terms“including” and“comprising” are used in an open-ended fashion, and thus should be interpreted to mean“including, but not limited to... Also, the term “couple” or“couples” is intended to mean either an indirect, direct, optical or wireless electrical connection. Thus, if a first device couples to a second device, that connection may be through a direct electrical connection, through an indirect electrical connection via other devices and connections, through an optical electrical connection, or through a wireless electrical connection.

DETAILED DESCRIPTION OF THE DISCLOSED EXEMPLARY EMBODIMENTS

[0030] The following discussion is directed to various exemplary embodiments of the invention. Although one or more of these embodiments may be preferred, the embodiments disclosed herein are merely examples and should not be interpreted, or otherwise used, as limiting the scope of the disclosure, including the claims. In addition, one skilled in the art will understand that the following description has broad application, and the discussion of any embodiment is meant only to be exemplary of that embodiment, and not intended to intimate that the scope of the disclosure, including the claims, is limited to that embodiment. In the following discussion and in the claims, the terms“including” and“comprising” are used in an open-ended fashion, and thus should be interpreted to mean“including, but not limited to... As used herein, the term“about,” when used in conjunction with a percentage or other numerical amount, means plus or minus 10% of that percentage or other numerical amount. For example, the term“about 80%,” would encompass 80% plus or minus 8%. References cited herein are incorporated in their entirety by such reference.

[0031] The present disclosure provides for new thermoelectric materials, i.e., Mg 3 Bi 2 - based Zintl phase compounds, wherein it has been discovered that the semi-metallic Mg 3 2 Bi 2 (nominal composition) described herein show an unexpectedly large Seebeck coefficient of -107 mV K '1 at 350 K and enhanced thermoelectric performance.

[0032] Such a high Seebeck coefficient may be attributed to the very large electron-to- hole weighted mobility ratio above 8. By partial substitution of Bi with Sb, the thermoelectric performance of Mg 3 2 Bi 2 -based materials was substantially improved. The optimized Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 shows a peak ZT of ~0.9 at 350 K. The n-type Mg 3 2 Bi 2 -based materials are proposed for thermoelectric cooling applications. A unicouple of n-type Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and p-type Bio . 5Sb1 . 5Te3 was thus prepared and a large temperature difference of ~91 K can be realized with the hot-side temperature of ~350 K. Such a cooling performance is even better than the commercial legs consisting of n-type Bi 2 Te 2. 7Se 0 . 3 and p-type Bio . 5Sb1 . 5Te3 (Mao, J.et al. , Science, Vol. 365, Issue 6452, pp. 495-498, 2019; Zu et al., Sci. Adv. 2019; 5, 5:eaav5813, are incorporated herein in its entirety).

[0033] Again, embodiments of Mg 3.2 Bi 2 -based thermoelectric materials demonstrate a high ZT of ~0.9 at 350 K. Such a high performance is better than the state-of-the-art n- type Bi 2 Te 2.7 Se 0 . 3 , and further minimization of the contact resistance remains to be a challenge for the nanostructured Bi 2 Te 2.7 Se 0 . 3 of the prior art. In comparison, it has been found herein, that Fe is a good contact material for n-type Mg 3.2 Bi 2 -based materials in order to fabricate the Fe/Mg 3.2 Bii .498 Sbo .5 Teo . oo 2 /Fe thermoelectric legs (Many other elements can also be used as contact materials such as Ni, Cu, etc.). Additionally, the Mg 3.2 Bi 2 -based materials are much cheaper than Bi 2 Te 2.7 Se 0 . 3 due to the higher price of Te. Therefore, replacing Bi 2 Te 2.7 Se 0.3 with the Mg 3.2 Bi 2 -based materials can significantly reduce the cost of thermoelectric module. [0034] Due to the excellent cooling performance of the n-type Mg3 .2 Bi 2 -based materials disclosed herein, and their advantages over Bi 2 Te 2.7 Se 0 . 3 , the n-type Mg 3 2 Bi 2 -based materials paired with the existing p-type Bi 2 Te 3 -based materials can potentially change the landscape of cooling applications by using thermoelectric modules.

[0035] The n-type Mg 3 2 Bi 2 -based materials with various compositions were thus prepared herein, and their thermoelectric properties characterized. The unicouple of n- type Mg 3.2 Bi-i . 498Sbo .5 Te 0.002 and p-type Bio .5 Sb 1.5 Te 3 were also prepared. The thermoelectric cooling measurement shows that a large temperature difference of ~91 K can thus be realized when the hot-side temperature is around 350 K.

[0036] In order to elucidate the technologies described herein: solid-state thermoelectric energy conversion holds great promises for power generation as well as thermal management 11 41 . The energy conversion efficiency of a thermoelectric module largely relies on the performance of the materials used, which is governed by the dimensionless figure of merit ( ZT ). ZT = S 2 T/K , where S is the Seebeck coefficient, a the electrical conductivity, k the thermal conductivity, and T the absolute temperature 15 11 ]

[0037] When a thermoelectric material is applied at a certain temperature, electrons in the valence band can be thermally activated into the conduction band, freeing an equal number of holes in the valence band. Since electrons and holes contribute oppositely to the Seebeck coefficient, it will be greatly diminished by such a bipolar conduction, which will also deteriorate the thermoelectric performance. Because the activation of electrons requires sufficient thermal energy (which depends on temperature) to cross the bandgap, the bandgap energy (E g ) plays a decisive role in bipolar conduction 1121 . As a result, most of the good thermoelectric materials reported thus far have been semiconductors (with a finite bandgap).

[0038] The optimum bandgap for thermoelectric materials has been predicted to be in the range between 6k B T and 10/c B T (/. e. , in the range between ~0.16 and ~0.26 eV at room temperature), where k B is the Boltzmann constant 113 151 . It should be noted that the conventional wisdom based on the assumption that the band effective masses ( m ) and carrier mobilities ( m ) are comparable for electrons and holes, /.e., m e * = m b * and m B = [13] . Flowever, in a more complicated scenario where the transport properties are asymmetrical, /.e., , the optimal bandgap will be different.

[0039] Such asymmetrical thermoelectric properties can be quantified by the electron-to- hole weighted mobility ratio where N is the valley

degeneracy. Therefore, the effect of bipolar conduction on thermoelectric performance when the transport properties are asymmetrical is of importance. When bipolar conduction occurs, the contributions from electrons and holes should both be considered and the Seebeck coefficient can be expressed as:

where S e and S h are the partial Seebeck coefficients for electrons and holes, respectively, and s e and sh are the partial electrical conductivities for electrons and holes, respectively. By applying the two-band model (/.e., considering one conduction band and one valence band), Eq.(1 ) can be rewritten as [11]

where F is the Fermi integral, e is the elementary charge, r is the scattering parameter, is the conduction band reduced Fermi energy, and is the reduced bandgap energy

[0040] According to Eq.(2), the Seebeck coefficient depends on the carrier concentration , bandgap electron-to-hole weighted mobility ratio (A), and carrier

scattering mechanism (r). It is well known that increasing the bandgap energy (i.e., increasing and increasing the concentration of majority carriers ( e.g ., increasing for n-type materials) can effectively suppress the bipolar conduction. In contrast, the importance of the electron-to-hole weighted mobility ratio ( A ) and carrier scattering mechanism (r) on bipolar conduction has not yet been fully recognized. In fact, preferential scattering of minority carriers (i.e., tuning the scattering mechanism) has been considered as a possible strategy for suppressing bipolar conduction, embodwe will mainly limit the discussion to the relationship among the electron-to-hole weighted mobility ratio, the Seebeck coefficient, and bipolar conduction.

[0041] Early studies addressing the weighted mobility ratio were mainly focused on its impact on the Goldsmid-Sharp bandgap 116 181 . To clearly demonstrate its effect on the Seebeck coefficient, the electron concentration and the electron-to-hole weighted mobility ratio dependent Seebeck coefficient were calculated herein for a n-type semi metal with a zero bandgap at 300 K.

[0042] As shown in Figure. 1 , the Seebeck coefficient is nearly independent of the electron-to-hole weighted mobility ratio when electron concentration is high (e.g., above 10 20 cm 3 ). In this case, the concentration of majority carriers (i.e., electrons) exponentially outnumbers that of the minority carriers (i.e., holes) and the Seebeck coefficient will be dominated by the contribution of electrons. The electron-to-hole weighted mobility ratio plays a vital role in determining the Seebeck coefficient when the electron concentration is relatively low (e.g., around 5 c 10 18 cm 3 ). At such a low electron concentration, the bipolar conduction is usually very noticeable, especially for a semi-metal. Indeed, the Seebeck coefficient is rather low (with an absolute value of less than 30 mV K 1 ) when the weighted mobility ratio is around unity, which indicates that the transport properties are symmetrical between the conduction and valence bands. This means that the contribution from electrons has been almost offset by that from holes, i.e., a strong bipolar conduction effect. In a sharp contrast, a high Seebeck coefficient (with an absolute value above 100 mV K 1 ) can be achieved when the transport properties are highly asymmetrical, i.e., A » 1 or A « 1. Specifically, a large negative Seebeck coefficient can be obtained for n-type materials when A » 1 and a large positive Seebeck coefficient can be realized by p-type materials when A « 1. In other words, highly asymmetrical transport properties are quite beneficial for suppressing bipolar conduction in the conduction type with good thermoelectric performance. Conversely, noticeable bipolar conduction will occur in the conduction type with inferior thermoelectric performance. Eventually, this will lead to highly asymmetrical thermoelectric performance between the n- and p-type of a given material (US 7,465,871 ; US 8,293, 168; and US 9,01 1 ,763 are herein incorporated in their entirety by reference).

[0043] Modelling was also conducted for semiconductors with different bandgap energies as shown in Figure 1 (b) and Figure 1 (c). Similarly, the weighted mobility ratio plays a pivotal role in the Seebeck coefficient of the semiconductors with a very narrow bandgap (Figure 1 (b)). When the bandgap energy is large enough (e.g., E g > 10k B T), bipolar conduction will be negligible and Seebeck coefficient is independent of the weighted mobility ratio (Figure 1 (b)). It should be reiterated that all of the modelling was performed at 300 K. At elevated temperature (e.g., 700 K or above), noticeable bipolar conduction will occur (for semi-metals and semiconductors with a very narrow bandgap) even when the weighted mobility ratio has an extreme value. Therefore, these models indicate that the semi-metals or semiconductors with a very narrow bandgap (e.g., E g less than ~0.1 eV) could demonstrate exceptional thermoelectric performance up to room temperature when the weighted mobility ratio has an extreme value.

[0044] Recently, n-type Mg 3 Sb 2 -based Zintl phases have emerged as highly promising thermoelectric materials [22 36] . Unlike the semiconducting Mg 3 Sb 2 , the isostructural Mg 3 Bi 2 is a semi-metal [37 39] . Experimentally, Mg 3 Bi 2 prepared with the exact stoichiometry shows p-type conduction (Figure 6 (b), which is in good agreement with the previous reports [37, 38] . Such a persistent p-type conduction should be attributed to the presence of a high concentration of Mg vacancies, which is similar to Mg 3 Sb 2 [22, 31 ] .

[0045] Flerein, n-type Mg 3 Bi 2 samples with excess Mg in nominal composition was successfully synthesized, i.e., Mg 3+d Bi 2 (d = 0.05, 0.1 , and 0.2, throughout the disclosure all the compositions are nominal to compensate the loss of Mg during sample preparation) and it was found that the thermoelectric properties are very similar to one another (Figure 6 (a)-(f), hence in one embodiment the base composition of Mg 3.2 Bi 2 was chosen for study.

[0046] The measured carrier concentration of Mg 3.2 Bi 2 is ~2.1×10 19 cm -3 at 10 K (Figure 2a, the upper panel) and it increases noticeably with temperature when the temperature is above 150 K. In addition, the carrier mobility of Mg 3.2 Bi 2 is more than 200 cm 2 V -1 s -1 when the temperature is below 300 K and it is as high as ~4198 cm 2 V -1 s -1 at 10 K (Figure 2a, the lower panel).

[0047] The high carrier concentration and mobility in conjunction with the low bipolar conduction temperature clearly indicate the semi-metallic characteristic of Mg 3.2 Bi 2 . Benefiting from the high carrier concentration and mobility, the electrical resistivity of Mg 3.2 Bi 2 is only ~9 mW m at 300 K and it is as low as ~0.58 mW m at 2 K (Figure 2 (b)). Surprisingly, the Seebeck coefficient of Mg 3.2 Bi 2 is more than -80 mV K -1 over a broad temperature range between 130 K and 350 K and a high Seebeck coefficient of ~-105 mV K -1 can be achieved at room temperature.

[0048] Elemental bismuth is a classic semi-metal with good thermoelectric performance (Figure 7 (a)-(d)) [40-44] . High-purity single-crystalline bismuth shows n-type conduction with a large room-temperature Seebeck coefficient around -100 mV K -1 along the trigonal axis and around -50 mV K -1 in the direction perpendicular to the trigonal axis [41, 42] .

[0049] Polycrystalline bismuth (light green symbols in Figure. 2(c) exhibits a Seebeck coefficient about -70 µV K -1 at 300 K, which is lower than that of the polycrystalline Mg 3.2 Bi 2 sample prepared in this work. Combined with the low electrical resistivity and the unexpectedly large Seebeck coefficient, Mg 3.2 Bi 2 can achieve a power factor ~10 µW cm -1 K -2 across a broad temperature range (Figure 2(d)). Additionally, minute Te doping can effectively modulate the carrier concentration of Mg 3.2 Bi 2 (Figure 8(a)-(b)), thus reducing the electrical resistivity and slightly enhancing the Seebeck coefficient (Figure 2b and Figure 2c, respectively). As a result, substantially higher power factors on the order of ~20 µW cm -1 K -2 over a broad temperature range can be realized in Te- doped Mg 3.2 Bi 2 (Figure 2d). Finally, a peak ZT around 0.3 at 350 K can be obtained by Mg 3.2 Bi 1.998 Te 0.002 and the ZT is above 0.1 down to 150 K (Figure 9(b)). The average ZT

  (calculated by the integration method) between 100 K and 350 K is ~0.2 for Mg 3.2 Bi 1.998 Te 0.002 and ~0.2 for Mg 3.2 B 1 . 996 Te 0.004 ·

[0050] To understand the electronic properties of Mg 3 Bi 2 , its band structure was calculated as shown in Figure 10 (a). The calculated band overlap energy between the conduction band minimum and valence band maximum is around -0.4 eV. However, the estimated band overlap from the two-band modelling is about -0.1 eV (Figure 1 1 ), therefore the calculated band structure with the band overlap energy shifted to -0.1 eV is shown in Figure 3a. It can be seen that the conduction band minimum is located along the L-M line with the valley degeneracy of 6 (Figure 12 (a)) and the valence band maximum locates around G point with the valley degeneracy of 1 (Figure 12(b)). It should be noted that different band structures were previously reported for Mg 3 Bi 2 [23, 39, 451 and the difference can be attributed to the choice of pseudo-potential and whether the spin-orbital coupling is considered 139, 45] . As mentioned above, the unexpectedly large Seebeck coefficient in the n-type semi-metallic Mg 3.2 Bi 2 is essential for its good thermoelectric performance. To reveal the underlying mechanism of the high Seebeck coefficient, a two-band model (i.e., one conduction band along the L-M line and one valence band around G point) was applied to calculate the partial Seebeck coefficients and carrier mobilities for Mg 3.2 Bi 2 [11] .

[0051] The partial Seebeck coefficients for electrons and holes of Mg 3.2 Bi 2 were calculated as shown in Figure 3b. The partial Seebeck coefficient of electrons is noticeably larger than that of holes. The difference in the partial Seebeck coefficients should be attributed to the disparity in the density-of-states effective mass

[0052] The calculated is ~0.530m 0 for the conduction band and ~0.276m 0 for the

valence band (Table 1 ).

[0054] The density-of-states effective mass depends on the band effective mass (m ) as well as the valley degeneracy ( N ) according to the relation of m = /V 2/ [12] . The different valley degeneracies in the conduction and valence bands (Figure 9(a)-(b)) can account for the disparity in the density-of-states effective mass. The partial carrier mobilities of Mg 3 2 Bi 2 (in the range between 100 K and 350 K) were further calculated as shown in Figure 3 (c). It should be pointed out that the calculated partial mobilities have relatively large uncertainties when temperature is below 100 K due to the small absolute value of Seebeck coefficient, therefore these values are not shown in Figure 3 (c). It can be seen that the electron mobility is noticeably higher than the hole mobility. When acoustic phonon scattering is the predominant scattering process, the carrier mobility is proportional to (m * ) 3/2 (w * ) \ where hz is the inertial effective mass 15, 12, 461 .

[0055] Therefore, low band effective mass and inertial effective mass are very beneficial for achieving a high carrier mobility. In fact, the calculated band effective mass is ~0.161 m 0 for electrons (m 0 is the free electron mass) and ~0.276m 0 for holes and the inertial effective mass (w * ) is ~0.133m 0 for electrons and ~0.259m 0 for holes.

[0056] Therefore, the higher mobility of electrons compared to that of holes is very reasonable. As a result, the electron-to-hole mobility ratio is around 3 and the weighted mobility ratio is above 8 (Figure 3 (d)). Such a large electron-to-hole weighted mobility ratio contributes to the high Seebeck coefficient of n-type Mg 3. 2Bi 2. Meanwhile, it also partially explains why the p-type Mg 3 Bi 2 shows a much lower Seebeck coefficient and inferior thermoelectric performance (Figure 6). In other words, the transport properties of Mg 3 Bi 2 are highly asymmetrical between the conduction and valence bands, thus leading to the noticeable difference in thermoelectric performance between the n- and p-type samples. It should also be noted that the electron-to-hole mobility ratio of single crystalline bismuth is ~9.19 along the trigonal axis [41 , 42, 47] and the density-of-states effective mass is ~0.1 13m 0 for electrons and ~0.093m 0 for holes [42] . As a result, the estimated electron-to-hole weighted mobility ratio is as large as ~1 1 along the trigonal axis for bismuth, which explains its high Seebeck coefficient along this direction. The highly asymmetrical transport properties of these semi-metals disclosed herein render embodiments thereof highly promising thermoelectric materials. The relatively high thermal conductivity (~2.8 W m 1 K 1 at 300 K, Figure 9(a)) of Mg 3 2 Bi 2 is the main limit to its thermoelectric performance. In fact, partial substitution of Bi with Sb in Mg 3 2 Bi 2 should substantially reduce the lattice thermal conductivity and improve the thermoelectric performance. Based upon the composition [g 3.2 Bi 1.998 Te 0.002 , solid solutions with the nominal composition of Mg 3.2 Bi 1.298 Sb 0.7 Te 0.002 (x = 0, 0.1 , 0.3, 0.5, and 0.7) were further prepared.

[0057]The electrical resistivity of Mg 3.2 Bi 1.998-x Sb x Te 0.002 gradually increases with increase Sb concentration (Figure 4a). The electrical resistivity is ~8.2 mW m at 300 K and ~2.2 mW m around 2 K for Mg3 .2 Bi 1.998 Te 0.002 · In comparison, it is ~24.4 mW m at 300 K and ~9.9 mW m around 2 K for Mg3 .2 Bi 1.298 Sb 0.7 Te 0.002 · The increased electrical resistivity can be attributed to the reduced carrier concentration (Figure 4 (d)) and carrier mobility (Figure 4(e)) after Sb alloying.

[0058] Similarly, the Seebeck coefficient is greatly enhanced with Sb alloying, as shown in Figure 4(b). The room-temperature Seebeck coefficient is ~-129 mV K 1 for Mg3 .2 Bii . 99 8 Te 0.002 and it is ~-229 mV K 1 for Mg3 .2 Bi 1.298 Sb 0.7 Te 0.002 , an enhancement of ~78%. The Seebeck coefficient enhancement should be ascribed to the reduced carrier concentration (Figure 4 (d) as well as the modified band structure after Sb alloying. Since Mg 3 Bi 2 is a semi-metal while Mg 3 Sb 2 is a semiconductor, it is very reasonable to expect a band structure transition from semi-metallic to semiconducting in Mg 3 2 Bii . 998- xSb x Te 0.002 solid solutions. To probe the variation in band structure after Sb alloying, the temperature-dependent electrical resistivity of undoped Mg 3.2 Bi 2.x Sb x (x = 0, 0.1 , 0.2, 0.3, 0.4, 0.5, 0.6, and 0.7) was studied. The temperature-dependent electrical resistivity of Mg 3.2 Bi 2 - x Sb x can be described by the expression of where E g is the band overlap energy (for a semi-metal) or the bandgap energy (for a semiconductor). Therefore, E g was determined for Mg 3.2 Bi 2-x Sb x as shown in Figure 4f (and derived by Figure13 (a)-(h)). Clearly, there is a semi-metallic to semiconducting transition with increasing Sb concentration.

[0059] The band overlap energy is about -0.013 eV for Mg 3.2 Bi 2 , and Mg 3.2 Bi 1.9 Sb 0.1 shows a nearly zero-gap (~0.005 eV). All of the samples with Sb concentration of more than 5% (i.e., x is above 0.1) are semiconducting and the bandgap increases with greater Sb concentration. The reported bandgap of Mg 3.2 Bi 0.9 Sb 1.1 is ~0.24 eV [35] and that of Mg 3.2 Bi 1.3 Sb 0.7 is estimated to be ~0.147 eV in this work. Considering the lower Sb concentration in this composition, the estimated bandgap seems to be reasonable. It is also noted that the band overlap of Mg 3.2 Bi 2 that estimated from the temperature- dependent electrical resistivity (around -0.013 eV) is smaller than the estimated value from the two-band modelling (around -0.1 eV, Figure 3 (a) and Figure 11) and the value derived from the calculated band structure (about -0.4 eV, Figure 10 (a)). It suggests that the band overlap energy of Mg 3 Bi 2 should be small. To unambiguously resolve this discrepancy, further characterization of the high quality single-crystalline Mg 3 Bi 2 is necessary.

[0060] Such a semi-metal to semiconductor transition after Sb alloying can also be understood from the variation in the temperature dependence of carrier concentration (Figure 4d). The carrier concentration of Mg 3.2 Bi 1.998 Te 0.002 and Mg 3.2 Bi 1.898 Sb 0.1 Te 0.002 increases noticeably with temperature above 125 K, which is caused by the activation of electron-hole pairs due to the bipolar conduction. Such a phenomenon is greatly suppressed in Mg 3.2 Bi 1.698 Sb 0.3 Te 0.002 and finally minimized in Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and Mg 3.2 Bi 1.298 Sb 0.7 Te 0.002 due to the increased bandgap energy. In addition, the reduced carrier concentration in Mg 3.2 Bi 1.998-x Sb x Te 0.002 with increasing Sb concentration also indicates the downward shifting of the Fermi energies due to the variation in band structure (i.e., upward movement of conduction band edges that opens the bandgap). Despite the semi-metal to semiconductor transition after Sb alloying, the power factors among Mg 3.2 Bi 1.998-x Sb x Te 0.002 samples studied are comparable (Figure 4 (c)). A slight enhancement in power factor can be observed in Mg 3.2 Bi 1.898 Sb 0.1 Te 0.002 (with a bandgap of ~0.005 eV) and Mg 3.2 Bi 1.698 Sb 0.3 Te 0.002 (with a bandgap of ~0.06 eV) compared to that of Mg 3.2 Bi 1.998 Te 0.002 . When the Sb concentration is above 25% (i.e., x is above 0.5), power factors at lower temperature are notably reduced due to the substantially reduced carrier mobility (Figure 4e). The average power factors of Mg 3.2 Bi 1.998-x Sb x Te 0.002 are further calculated by the integration method between 2 K and 350 K as shown in Figure 14. High average power factors (above 16 µW cm -1 K -2 ) can be achieved in the semi-metallic Mg 3.2 Bi 1.998 Te 0.002 and semiconducting Mg 3.2 Bi 1.998- x Sb x Te 0.002 (x < 0.7) with very a narrow bandgap (less than 0.12 eV). As indicated by the above-mentioned modellings, semi-metals and semiconductors with a very narrow bandgap can achieved promising thermoelectric performance when the transport properties are highly asymmetrical (i.e., the large electron-to-hole weighted mobility ratio in Mg 3.2 Bi 2 -based materials).

[0061] A substantial reduction in the thermal conductivity of Mg 3.2 Bi 1.998-x Sb x Te 0.002 was observed with increasing Sb concentration (Figure 15). Such a reduction should be mainly attributed to the significantly reduced lattice thermal conductivity ( ^ lat ) as shown in Figure 5a. The peak lattice thermal conductivity is as high as ~9.6 W m -1 K -1 for Mg 3.2 Bi 1.998 Te 0.002 and it is only ~3.5 W m -1 K -1 for Mg 3.2 Bi 1.298 Sb 0.7 Te 0.002 , a reduction of ~65%. Benefitting from the largely reduced lattice thermal conductivity, substantial enhancements in the thermoelectric figure of merit have been realized (Figure 5b). The ZT at 350 K is around 0.3 for Mg 3.2 Bi 1.998 Te 0.002 and it is more than 0.9 for Mg 3.2 Bi 1.298 Sb 0.7 Te 0.002 . Although the thermoelectric performance of Mg 3.2 Bi 1.998- x Sb x Te 0.002 decreases with reduced temperature, both Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and Mg 3.2 Bi 1.298 Sb 0.7 Te 0.002 maintain a ZT above 0.1 down to 120 K. The average ZT (calculated by the integration method) between 100 K and 350 K is ~0.4 for Mg 3.2 Bi 1.298 Sb 0.7 Te 0.002 , ~0.5 for Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 , and ~0.2 for Mg 3.2 Bi 1.998 Te 0.002 . Therefore, the n-type Mg 3.2 Bi 2 -based materials are very promising for application in low- temperature thermoelectric cooling.

[0062] To validate such a potential, a unicouple consisting of n-type Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and p-type Bi 0.5 Sb 1.5 Te 3 (its low-temperature thermoelectric performance can be found at Figure 16) was prepared and a cooling measurement was conducted. The experimental setup for the cooling measurement is shown in Figure 5 (c), where the thermoelectric legs have been separately soldered onto two copper

  blocks and bridged by a copper plate. The copper blocks were placed on a heater to maintain the hot-side temperature around 350 K during the cooling measurement. A silicon wafer was placed in between the copper blocks and heater for electrical insulation. By applying electrical current to the unicouple, heat is removed from the top copper plate (i.e., the cold-side) and dissipates into the bottom copper blocks (i.e., the hot-side) due to the Peltier effect [2] . As a result, a temperature difference will exist between the cold and hot sides. As shown in Figure 5d, the temperature difference between the hot and cold sides first increases with increasing applied electrical current and then saturates at high electrical current. Finally, a large temperature difference of ~91 K can be obtained with the hot-side temperature maintained around 350 K. Comparison of the hot-side temperature dependent maximum temperature difference (DT) between this unicouple and the commercial thermoelectric legs is further shown in the inset of Figure 5d. Clearly, the cooling performance of the unicouple consisting of n- type Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and p-type Bi 0.5 Sb 1.5 Te 3 is better than the commercial thermoelectric legs consisting of n-type Bi 2 Te 2.7 Se 0.3 and p-type Bi 0.5 Sb 1.5 Te 3 . It should be noted that the dimensions of the n-type Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and p-type Bi 0.5 Sb 1.5 Te 3 legs are optimized according to their thermoelectric properties at 350 K, therefore the cooling performance of the unicouple at the hot-side temperature of 300 and 325 K is not optimized. It is worth pointing out that the contact to nanostructured n- type Bi 2 Te 2.7 Se 0.3 is notoriously difficult [49] , while Fe is a good contact material for n-type Mg 3.2 Bi 2 -based materials (Figure 17) and it is easy to fabricate the thermoelectric leg of Fe/Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 /Fe.

[0063] Additionally, the Mg 3.2 Bi 2 -based materials are much cheaper than Bi 2 Te 2.7 Se 0.3 due to the higher price of Te. Therefore, replacing Bi 2 Te 2.7 Se 0.3 with the Mg 3 Bi 2 -based materials can significantly reduce the cost of thermoelectric module. Due to the excellent cooling performance of the n-type Mg 3.2 Bi 2 -based materials and their advantages over Bi 2 Te 2.7 Se 0.3 , the n-type Mg 3.2 Bi 2 -based materials can potentially change the landscape of cooling applications by using thermoelectric modules. Conclusions: In this work, it has been demonstrated that the weighted mobility ratio, which can be used to quantify the asymmetry of thermoelectric transport properties between the conduction and valence bands, can play a decisive role in bipolar

  conduction. When the weighted mobility ratio has an extreme value (either A >> 1 or A <<1), the effect of bipolar conduction on the Seebeck coefficient can be substantially minimized. Therefore, semi-metals and semiconductors with a very narrow bandgap can possibly achieve a large Seebeck coefficient around ambient temperature when the weighted mobility ratio has an extreme value. Experimentally, it has been demonstrated herein that the semi-metallic Mg 3.2 Bi 2 shows an unexpectedly high Seebeck coefficient of above -100 µV K -1 around 300 K and promising thermoelectric performance. The optimized Mg 3.2 Bi 1.298 Sb 0.7 Te 0.002 can realize a peak ZT above 0.9 at 350 K and ZT above ~0.1 down to the temperature of 120 K.

[0064] These results show that n-type Mg 3.2 Bi 2 -based materials are highly promising for thermoelectric cooling applications. A unicouple consisting of n-type Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and p-type Bi 0.5 Sb 1.5 Te 3 can realize a large temperature difference of ~91 K with the hot-side temperature of ~350 K, which is better than the commercial thermoelectric legs. Importantly, it has been demonstrated herein that semi- metals and semiconductors (with a very narrow bandgap) with extreme weighted mobility ratios could exhibit outstanding thermoelectric performance. This work greatly expands the scope for the future discovery of promising thermoelectric materials.

Experimental Section

[0065] Synthesis: High-purity magnesium turnings (Mg, 99.98%; Alfa Aesar), antimony shots (Sb, 99.8%; Alfa Aesar), bismuth pieces (Bi, 99.99%; Alfa Aesar), and tellurium pieces (Te, 99.999%; Alfa Aesar) were weighed according to the composition of Mg 3+d Bi 2 (d = 0, 0.05, 0.1, and 0.2), Mg 3.2 Bi 2-x Te x (x = 0, 0.002, and 0.004), Mg 3.2 Bi 1.998- x Sb x Te 0.002 (x = 0, 0.1, 0.3, 0.5, and 0.7), and Mg 3.2 Bi 2-x Sb x (x = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, and 0.7). The elements were loaded into a stainless-steel ball-milling jar in a glove box under an argon atmosphere with an oxygen and water level below 0.1 ppm. The materials were ball milled for 10 hours. The ball-milled powders were then loaded into a graphite die with an inner diameter of 12.7 mm in a glove box. The graphite die with the loaded powders was immediately hot pressed at 1073 K for 2 minutes. The hot-pressed disks are about 2-3 mm thick. The polycrystalline bismuth sample was prepared by ball milling the bismuth pieces (Bi, 99.99%; Alfa Aesar) for one minute and then hot- pressing at 513 K for 5 minutes.

  [0066] Phase composition: Phase compositions of the samples were characterized by X-ray diffraction (PANalytical X’pert PRO diffractometer). All of the samples prepared in this work are single-phase (Figure 18).

[0067] Thermoelectric property characterizations: Samples were cut into bar-shaped pieces with the dimensions of ~3.5 × ~3.5 × ~7 mm 3 . Thermoelectric transport properties of each sample were measured from 350 K to 2 K at high vacuum (below 1×10 -5 Torr) via the thermal transport option in a physical properties measurement system (PPMS) by Quantum Design. The four-probe lead configuration method was used to avoid the influence of lead contact resistance. The gold-plated copper leads were attached to the sample by using the silver-filled H20E epoxy. To ensure good mechanical bonding, samples with attached leads were heated on the hot-plate at 400 K for 10 minutes. The Hall coefficient R H was obtained from the magnetic field- dependent Hall resistivities via the electrical transport option in PPMS. The carrier concentration n H was calculated by n H ^1 eR H .

[0068] Thermoelectric cooling measurement: To characterize the thermoelectric cooling performance, a unicouple consisting of n-type Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and p-type Bi 0.5 Sb 1.5 Te 3 was prepared. To ensure low electrical contact resistance, it is necessary to electroplate Ni on both sides of Bi 0.5 Sb 1.5 Te 3 . In addition, Fe contact layers are necessary for Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 to realize a low contact resistance (Figure 17) [50] . The Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 powders were sandwiched between Fe powders and then hot pressed at 1073 K for 2 minutes. To realize the maximum cooling performance, the dimensions of thermoelectric legs needs to be considered according to the relation of ^ A l ^

^ 1 1 ^

^ ^ 2 ^ 2 (based on the thermoelectric properties at 350 K), where A is the a and l is the leg length. The Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 is ~4.01 mm thick and the two Fe contact layers together are ~0.82 mm thick, so the total length of the Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 leg is ~4.83 mm. The Bi 0.5 Sb 1.5 Te 3 leg is also ~4.83 mm. The cross-section area is ~3.52 × ~3.76 mm 2 for Mg 3.2 Bi 1.498 Sb 0.5 Te 0.002 and ~3.45× ~2.99 mm 2 for Bi 0.5 Sb 1.5 Te 3 . The thermoelectric legs were soldered to the copper blocks and the copper plate by using Pb 37 Sn 63 foils (with a melting point of ~430 K). In fact, Ni was also electroplated on the copper blocks and copper plate to minimize the surface oxidation during the soldering process. The cooling measurement was performed at vacuum (below 10 4 Torr).

[0069] Theoretical calculation: All of the first-principles calculations were performed by using the Quantum ESPRESSO package 1511 . A fully relativistic projector-augmented- wave (PAW) pseudopotential with Perdew-Burke-Ernzerhof (PBE) for exchange- correlation energy functional was employed to calculate the electronic Hamiltonian. Monkhorst-Pack k-point meshes with 8 c 8 c 6 and 80 c 80 c 60 configurations were employed in the self-consistent and Fermi-surface calculations respectively. The finite difference method was employed to calculate the effective mass tensor, of which the eigenvalues were evaluated.

[0070] In conclusion, in some embodiments herein, the optimized Mg 3.2 Bi2-based thermoelectric materials demonstrate a high ZT of ~0.9 at 350 K. Such a high performance is better than the state-of-the-art n-type Bi 2 Te 2. 7Se 0. 3. However, minimization of the contact resistance remains to be a challenge for the nanostructured Bi 2 Te 2.7 Seo . 3. In comparison, it has been found herein that Fe is a good contact material for n-type Mg 3.2 Bi 2 -based materials and it is very easy to fabricate the Fe/Mg3 .2 Bi 1. 498Sbo .5 Te 0.002 /Fe thermoelectric legs (Many other elements can also be used as contact materials such as Ni, or Cu).

[0071] Additionally, the Mg 3.2 Bi 2 -based materials are much cheaper than Bi 2 Te 2.7 Seo.3 due to the higher price of Te. Therefore, replacing Bi 2 Te 2.7 Se 0 .3 with the Mg 3.2 Bi 2 -based materials can significantly reduce the cost of thermoelectric module. Due to the excellent cooling performance of the n-type Mg 3 2 Bi 2 -based materials and their advantages over Bi 2 Te 2.7 Se 0. 3, the n-type Mg 3.2 Bi 2 -based materials paired with the existing p-type Bi 2 Te 3 -based materials can potentially change the landscape of cooling applications by using thermoelectric modules.

[0072] Herein the following references are incorporated in their entirety:

[1 ] A. F. Ioffe, Semiconductor thermoelements and thermoelectric cooling, Infosearch, London 1957; [2] H. Goldsmid, Thermoelectric refrigeration, Springer US, New York 1964; [3] F. J. DiSalvo, Science 1999, 285, 703; [4] L. E. Bell, Science 2008, 321, 1457; [5] Y. Z. Pei, H. Wang, G. J. Snyder, Adv. Mater. 2012, 24, 6125. [6] G. J. Tan, L. D. Zhao, M. G. Kanatzidis, Chem. Rev. 2016, 116, 12123. [7] X. Shi, L. D. Chen, C. Uher, Int. Mater. Rev. 2016, 61, 379. [8] J. He, T. M. Tritt, Science 2017, 357, eaak9997. [9] T. J. Zhu, Y T Liu, C. G. Fu, J. P. Heremans, J. G. Snyder, X. B. Zhao, Adv. Mater. 2017, 29, 1605884. [10] Y. Zhou, L. D. Zhao, Adv. Mater. 2017, 29, 1702676. [11 ] J. Mao, Z. Liu, J. Zhou, H. Zhu, Q. Zhang, G. Chen, Z. Ren, Adv. Phys. 2018, 67, 69. [12] H. J. Goldsmid, Introduction to thermoelectricity, Vol. 121 , Springer Science & Business Media, 2009. [13] R. Chasmar, R. Stratton, Int. J. Electron. 1959, 7, 52.[14] G. D. Mahan, J. Appl. Phys. 1989, 65, 1578. [15] J. O. Sofo,

G. D. Mahan, Phys. Rev. B 1994, 49, 4565. [16] H. J. Goldsmid, J. W. Sharp, J. Electron. Mater. 1999, 28, 869.[17] J. Schmitt, Z. M. Gibbs, G. J. Snyder, C. Felser, Mater. Horiz. 2015, 2, 68. [18] Z. M. Gibbs, H.-S. Kim, H. Wang, G. J. Snyder, Appl. Phys. Lett. 2015, 106, 022112. [19] J. H. Bahk, A. Shakouri, Appl. Phys. Lett. 2014, 105, 052106. [20] S. Wang, J. Yang, T. Toll, J. Yang, W. Zhang, X. Tang, Sci. Rep. 2015, 5, 10136. [21 ] J. H. Bahk, A. Shakouri, Phys. Rev. B 2016, 93, 165209. [22] H. Tamaki, H. K. Sato, T. Kanno, Adv. Mater. 2016, 28, 10182. [23] J. Zhang, L. Song, S. H. Pedersen,

H. Yin, L. T. Hung, B. B. Iversen, Nat. Commun. 2017, 8, 13901 . [24] J. Mao, J. Shuai, S. W. Song, Y. X. Wu, J. W. Zhou, Z. H. Liu, J. F. Sun, Y. Z. Pei, Q. Y. Zhang, D. J. Singh, G. Chen, Z. F. Ren, Proc. Natl. Acad. Sci. U.S.A. 2017, 114, 10548. [25]J. W. Zhang, L. R. Song, A. Mamakhel, M. R. V. Jorgensen, B. B. Iversen, Chem. Mater. 2017, 29, 5371. [26] J. Mao, Y. Wu, S. Song, Q. Zhu, J. Shuai, Z. Liu, Y. Pei, Z. F. Ren, ACS Energy Lett. 2017, 2, 2245. [27]J. Shuai, J. Mao, S. W. Song, Q. Zhu, J. F. Sun, Y. M. Wang, R. He, J. W. Zhou, G. Chen, D. J. Singh, Z. F. Ren, Energy Environ. Sci. 2017, 10, 799. [28] J. Mao, Y. X. Wu, S. W. Song, J. Shuai, Z. H. Liu, Y. Z. Pei, Z. F. Ren, Mater. Today Phys. 2017, 3, 1 .[29] J. Zhang, L. Song, M. Sist, K. Tolborg, B. B. Iversen, Nat. Commun. 2018, 9, 4716. [30] T. Kanno, H. Tamaki, H. K. Sato, S. D. Kang, S. Ohno, K. Imasato, J. J. Kuo, G. J. Snyder, Y. Miyazaki, Appl. Phys. Lett. 2018, 112, 033903. [31 ] S. Ohno, K. Imasato, S. Anand, H. Tamaki, S. D. Kang, P. Gorai, H. K. Sato, E. S. Toberer, T. Kanno, G. J. Snyder, Joule 2018, 2, 141. [32] K. Imasato, M. Wood, J. J. Kuo, G. J. Snyder, J. Mater. Chem. A 2018. [33] J. J. Kuo, S. D. Kang, K. Imasato, H. Tamaki, S. Ohno, T. Kanno, G. J. Snyder, Energy Environ. Sci. 2018, 11, 429. [34] X. Chen, H. Wu, J. Cui, Y. Xiao, Y. Zhang, J. He, Y. Chen, J. Cao, W. Cai, S. J. Pennycook, Z. Liu, L.-D. Zhao, J. Sui, Nano Energy 2018, 52, 246.[35] R. Shu, Y. Zhou, Q. Wang, Z. Han, Y. Zhu, Y. Liu, Y. Chen, M. Gu, W. Xu, Y. Wang, W. Zhang, L. Huang, W. Liu, Adv. Fun . Mater. 2019, 29, 1807235. [36] K. Imasato, S. D. Kang, G. J. Snyder, Energy Environ. Sd. 2019. [37]V. Ponnambalam, D. T. Morelli, J. Eledron. Mater. 2013, 42, 1307. [38] J. Xin, G. Li, G. Auffermann, H. Borrmann, W. Schnelle, J. Gooth, X. Zhao, T. Zhu, C. Felser, C. Fu, Mater. Today Phys. 2018, 7, 61 . [39]T. R. Chang, I. Pletikosic, T. Kong, G. Bian, A. Huang, J. Denlinger, S. K. Kushwaha, B. Sinkovic, H. T. Jeng, T. Valla, W. Xie, R. J. Cava, Adv. Sci. 2019, 6, 1800897. [40] R. W. Boydston, Phys. Rev. 1927, 30, 911 .[41 ]B. Chandrasekhar, J. Phys. Chem. Solids 1959, 11, 268.[42] C. F. Gallo, B. S. Chandrasekhar, P. H. Sutter, J. Appl. Phys. 1963, 34, 144. [43]H. Goldsmid, Phys. Status Solidi A 1970, 1, 7. [44] H. J. Goldsmid, "Bismuth: the thermoelectric material of the future?", presented at 25th International Conference on Thermoeledrics, 2006. [45] X. Zhang, L. Jin, X. Dai, G. Liu, J. Phys. Chem. Lett. 2017, 8, 4814. [46] Y. Z. Pei, A. D. LaLonde, H. Wang, G. J. Snyder, Energy Environ. Sci. 2012, 5, 7963. [47] B. Abeles, S. Meiboom, Phys. Rev. 1956, 101, 544.[48] A. L. Jain, Phys. Rev. 1959, 114, 1518.[49] W. Liu, H. Wang, L. Wang, X. Wang, G. Joshi, G. Chen, Z. Ren, J. Mater. Chem. A 2013, 1, 13093.[50] Q. Zhu, S. Song, H. Zhu, Z. Ren, J. Power Sources 2019, 414, 393. [51 ] P. Giannozzi, S. Baroni, N. Bonini, M. Calandra,

R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococcioni, I. Dabo, J. Phys. Conden. Matter 2009, 21, 395502.U.S. 7,465871 ; U.S. 9,011763; U.S. 8,293168.