Next Article in Journal
Structure-Function Studies of Glucose Oxidase in the Presence of Carbon Nanotubes and Bio-Graphene for the Development of Electrochemical Glucose Biosensors
Previous Article in Journal
Photoluminescence of Argan-Waste-Derived Carbon Nanodots Embedded in Polymer Matrices
Previous Article in Special Issue
Deformation Mechanism of Aluminum, Copper, and Gold in Nanoimprint Lithography Using Molecular Dynamics Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermoelectric Properties of Mg3(Bi,Sb)2 under Finite Temperatures and Pressures: A First-Principles Study

1
School of Science, Harbin Institute of Technology, Shenzhen 518055, China
2
The State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China
3
Beijing MaiGao MatCloud Technology Co., Ltd., Beijing 100190, China
4
Guangdong Aerospace Research Academy, Guangzhou 511458, China
5
Computer Network Information Center, Chinese Academy of Sciences, Beijing 100190, China
6
University of Chinese Academy of Sciences, Beijing 100190, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Nanomaterials 2024, 14(1), 84; https://doi.org/10.3390/nano14010084
Submission received: 15 November 2023 / Revised: 14 December 2023 / Accepted: 21 December 2023 / Published: 28 December 2023
(This article belongs to the Special Issue Mathematical and Computational Modeling for Nanohybrids)

Abstract

:
Mg3Bi2−vSbv (0 ≤ v ≤ 2) is a class of promising thermoelectric materials that have a high thermoelectric performance around room temperatures, whereas their thermoelectric properties under pressures and temperatures are still illusive. In this study, we examined the influence of pressure, temperature, and carrier concentration on the thermoelectric properties of Mg3Bi2−vSbv using first-principle calculations accompanied with Boltzmann transport equations method. There is a decrease in the lattice thermal conductivity of Mg3Sb2 (i.e., v = 2) with increasing pressure. For a general Mg3Bi2−vSbv system, power factors are more effectively improved by n-type doping where electrons are the primary carriers over holes in n-type doping, and can be further enhanced by applied pressure. The figure of merit (zT) exhibits a positive correlation with temperature. A high zT value of 1.53 can be achieved by synergistically tuning the temperature, pressure, and carrier concentration in Mg3Sb2. This study offers valuable insights into the tailoring and optimization of the thermoelectric properties of Mg3Bi2−vSbv.

1. Introduction

In recent decades, the pressing challenges of climate change and pollution have underscored the critical importance of alternative energy sources. The increasing demand for energy, especially electricity, has become a central focus in the development of our civilization. To address these challenges sustainably, there is growing necessity for innovative renewable energy sources that can harness the power of nature and convert it into electricity [1]. Thermoelectric materials, which have the ability to convert heat energy into electrical energy and vice versa, are considered a promising approach to improving energy conversion efficiency [2,3]. Materials with exceptional thermoelectric properties are often used in waste heat recovery and refrigeration applications [4,5].
In general, the Seebeck coefficient ( S ) is defined as the ratio of the potential difference under zero current conditions to the applied temperature difference. On the other hand, the electrical conductivity ( σ ) is used to describe the system’s electrical performance. We can use the power factor ( S 2 σ ) to accurately characterize the thermoelectric properties of the system [6,7]. To maintain the temperature gradient of the material, we also need to consider the material’s thermal conductivity, which is generally divided into lattice thermal conductivity ( κ l ) and electronic thermal conductivity ( κ e ). The former is related to the transfer of phonons, while the latter is proportional to the electrical conductivity. An excellent thermoelectric material needs to have a higher power factor and lower thermal conductivity [8,9], which ensures a higher thermoelectric conversion efficiency. Quantitatively, we typically use z T = S 2 σ T / ( κ l + κ e ) , to describe the thermoelectric performance of the system [10,11]. The higher the zT value, the greater the thermoelectric conversion efficiency. In recent years, emerging nanomaterials such as TiS3 [12], MoS2 [13], and CNTs [14] have been reported to exhibit excellent thermoelectric properties.
Certain theoretical studies propose that the thermoelectric properties of materials may be modified through the application of stress or pressure [15,16,17,18]. Under the influence of applied pressure, some systems can undergo changes in electrical conductivity, resulting in a transition between insulating and metallic states [19,20]. In general, the energy difference between the bottom of the conduction band (CBM) and the top of the valence band (VBM), known as the bandgap, is a crucial parameter that describes the electronic structure of a system. When pressure is applied, both the CBM and the VBM shift, resulting in a change in the bandgap. This change in the bandgap also results in variations in the electrical transport coefficients, such as the Seebeck coefficient and the electrical conductivity. In the case of intrinsic semiconductors, the correlation between electrical conductivity and the Seebeck coefficient is given by [21,22]:
σ = σ 0 exp ( E g 2 k T )
S = k | e | [ σ n σ p σ n + σ p E g 2 k T + ( r n + 5 2 ) σ n σ n + σ p ( r p + 5 2 ) σ p σ n + σ p 3 4 l n m p * m n * ]
where σ is the total electrical conductivity, σ n and σ p are the partial conductivity for electrons and holes, E g is the bandgap, k is Boltzmann’s constant, T is the temperature, e is the electron charge, r n ( r p ) and m n * ( m p * ) are the scatterging parameters and effective masses of the density of states of electrons (holes), respectively. A noticeable trend is the decrease in the Seebeck coefficient as the bandgap decreases, accompanied by an increase in electrical conductivity. Consequently, in materials where the bandgap is affected by the application of pressure, such as PbTe, PbSe, and Bi2Te3 [23,24,25], the power factors also demonstrate notable variations under pressure.
In recent years, there has been a notable surge in research interest surrounding Mg3Bi2−vSbv as a promising material for thermoelectric applications at room temperature. According to the study by Peng et al. [26], the Mg3Bi2−vSbv system demonstrates significantly low lattice thermal conductivity. Kanno, T. et al. [27] indicated that disorder in the system enhances carrier mobility and reduces the lattice thermal conductivity. Imasato et al.’s [28] study demonstrated that alloying in the Mg3Bi2 system alters its electronic transport properties, thereby influencing the thermoelectric properties. Pan et al.’s study [29] revealed that Mg3Bi1.25Sb0.75, with moderate doping, achieves a higher carrier mobility. Shi et al.’s research [30] revealed that the high carrier mobility in Mg3Sb2 is attributed to the purification of phases and the presence of coarse grains. Some experiments and computational studies suggest that Mg3Sb1.5Bi0.5, with n-type doping, exhibits excellent thermoelectric properties at certain temperatures [31,32,33,34]. Additionally, there are studies that describe how the thermoelectric performance of the system can be experimentally optimized [35,36,37]. The corresponding experimental reports also exist for p-type doping in Mg3Sb2 [38]. Compared to conventional thermoelectric materials like Bi2Te3 [39], Ag2Se [40], and GeTe [41], Mg3Bi2−vSbv presents the benefit of being cost-effective. Extensive research has been conducted into the different aspects of the Mg3Bi2−vSbv system. These investigations focused on exploring band topology [42,43,44,45], phonon dynamics [46,47,48,49,50], mechanical properties [51,52], and the topological thermoelectric properties of nodal-line semimetals [53,54,55,56]. According to previous studies [52], van der Waals (vdW) corrections, which are a specific type of intermolecular force that occur between neutral atoms or molecules, also play a significant role in influencing the properties of the system. The comprehension of the alterations in the thermoelectric characteristics of Mg3Bi2−vSbv under specific pressure conditions is of the utmost importance, given the potential of this material as a thermoelectric substance. Nevertheless, there is a dearth of research in this particular field.
In this study, we conducted computational calculations to determine the lattice thermal conductivity, Seebeck coefficient, electrical conductivity, and power factor of Mg3Bi2−vSbv at various temperatures and pressures. Specifically, in the context of lattice thermal conductivity, we utilized various methodologies to simulate the impacts of vdW corrections and subsequently compared the obtained outcomes. We have undertaken a comprehensive analysis to investigate the impact of pressure and temperature on the thermoelectric properties of the system.

2. Materials and Methods

2.1. Structure of Mg3Bi2−vSbv

The crystal structure of Mg3Bi2−vSbv is depicted in Figure 1. The crystal structures of all concentrations within the system exhibit the trigonal crystal system, characterized by a symmetry group of P 3 m - 1. In these systems, the three crystal axes, denoted by a and b, exhibit equal lengths, and the angles formed between them are α = β = 90°, with γ being 120°. Sb and Bi atoms are located in identical positions within the system. By enlarging the unit cell and subsequently replacing Bi atoms with Sb atoms, structures with varying Sb compositions can be obtained. It is apparent that, when considering the same Sb content, there are typically multiple configurations present. We utilize the LAsou [57,58] method to determine the lowest energy configuration, which serves as the representative structure for the given Sb content. In the LAsou method, a minimal number of first-principles calculations is necessary to determine the lowest energy configuration, given a specific Sb content.

2.2. Method of DFT Calculations

The Vienna Ab initio Simulation Package (VASP) [59,60,61,62] was employed for conducting our primary first-principles calculations, with a cutoff energy of 450 eV being set. To streamline the production of K-point input files, we utilized VASPKIT [63]. The K-points were generated using the Methfessel–Paxton [64] method, which is a technique utilized by VASP to generate reciprocal space K-point grids. We configured the grid density to a value of 0.03. The Perdew–Burke–Ernzerhof (PBE) [62] functional, which is a commonly used exchange-correlation functional, was utilized for the computational calculations. It was determined that these specific parameter configurations yielded convergent and resource-efficient calculations.
We employ two approaches, namely PBE-D3 [65] and vdW-DFq [66], to count the non-local vdW corrections. The computational method DFT-D3, developed by Grimme et al., is used to characterize the vdW dispersion energy-correction terms. The incorporation of empirical parameters in density functional theory has been proposed as a means to address its limitations, particularly in calculations involving heavy elements. This approach has demonstrated promising outcomes in terms of accuracy and reliability. Another method, known as vdW-DFq, has been developed by Peng et al., which demonstrates high accuracy in calculating the density and geometry of semihard materials. Both of these methods can be readily incorporated into VASP calculations by adjusting the relevant calculation parameters in VASP.

2.3. Method of Electrical Transport Properties

For the electrical transport properties, we utilized the relaxation time approximation (RTA) to solve the Boltzmann equation [67] with the conductivity and Seebeck coefficient expressed as follows:
σ α β ( μ , T ) = 1 V n k v n k α v n k β τ n k [ f μ ( ε n k ,   T ) ε n k ]
S α β ( μ , T ) = 1 e T V σ α β ( μ , T ) 1 n k v n k α v n k β τ n k ( μ ε n k ) [ f μ ( ε n k ,   T ) ε n k ]
where v n k is electron group velocity corresponding to band index n and the reciprocal coordinate k, T, μ, V, fu and e are the absolute temperature, Fermi level, volume of unite cell, the Fermi–Dirac distribution, and electron charge, respectively. ε n k is the band energy, and τ n k is the electronic relaxation time. The TransOpt [68] package implements such calculations. Regarding the treatment of the relaxation time, we used two methods. One assumes a constant relaxation time (CRTA). In this approximation, the expression for the Seebeck coefficient does not include the relaxation time. By comparing the experimentally measured conductivity of a system with a specific carrier concentration and the calculated value of σ / τ , the relaxation time can be determined. The determined relaxation time is then applied to the results for other carrier concentrations to obtain the corresponding conductivity and electronic thermal conductivity. The second method involves introducing constant-time electron–phonon coupling (CEPCA) in TransOpt. It considers the primary scattering mechanism of electron–phonon coupling in the electrical transport process, treating the electron–phonon coupling matrix as a constant. And, the relaxation time is approximated as follows:
τ n k 1 = 2 π k B E d e f 2 V G m k δ ( ε n k ε n k )
Here, G is the Young’s modulus, E d e f is the DP of the band edge state. We can input such a value to TransOpt and it can deal τ n k by itself.

2.4. Method of Thermal Properties

We utilized the finite displacement method offered by PHONOPY [69] to calculate the phonon spectrum, phonon density of states, and second-order force constants for the system. The supercell of the system was set to 3 × 3 × 2 . The lattice thermal conductivity for the system was calculated using the SHENGBTE-v1.5.0 software, and third-order force constants were computed in conjunction with VASP using the ThirdOrder program [70,71]. We employed the third-order program to expand the unit cell and generated a large number of perturbed configurations. Subsequently, we used VASP to calculate the energies of these configurations. By combining these energies, we obtained the third-order force constants. In the calculation of force constants, we introduced a parameter called the cutoff neighbor n. This parameter represents the distance of the farthest n-th neighbors and is used as the cutoff radius. Interactions beyond this radius will be disregarded. This parameter needs to be tested for convergence. The supercell size was maintained at 3 × 3 × 2 . Due to the extensive number of configurations to be calculated, only the Γ-point was employed to reduce the computational overhead. To investigate the impact of vdW corrections on the thermal properties of the system, we compared the results obtained using two vdW simulation methods, namely PBE-D3 and vdW-DFq, with the results obtained without any vdW contributions using the PBE functional.
Particularly, due to the complex third-order force constant calculations, it is necessary to compute a significant number of structural configurations, often in the thousands. Ensuring the completeness of the VASP calculations for these configurations, and handling the convergence of the near-neighbor cutoff, is a demanding task. It requires generating a large number of perturbed structural configurations with various cutoff parameters, leading to a substantial number of redundant structures. Calculating each of them individually would be a wasteful consumption of computational resources. To address this challenge, we developed and integrated a workflow, as shown in Figure S1, using the high-throughput materials cloud platform, MatCloud [72]. This workflow allows us to set supercell sizes and the cutoff neighbor parameters using the structure generator node. This node can automatically drive the ThirdOrder program to generate the perturbed structure needed for calculation. The parallel controller node orchestrates high-performance computing (HPC) resources to carry out the VASP calculations for all generated configurations.
It is important to note that this workflow assigns a unique ID to each generated configuration file based on the MD5 hash value of its structure. The MD5 algorithm has an extremely low probability of hash collisions, ensuring uniqueness. In other words, structures with the same hash value have identical file contents. The entire computational result is stored and indexed according to the hash values of the configuration files. This approach offers a significant advantage. When a specific near-neighbor cutoff parameter is set, we attempt to generate third-order force constants for all parameters within this cutoff after the calculations. However, in cases where it is found that third-order force constants involving set atomic near-neighbors are insufficient to achieve the convergence of lattice thermal conductivity, increasing the near-neighbor cutoff parameter allows us to calculate the additional perturbed structures. During post-processing, we retrieve previously calculated structures from the entire database to generate results, thereby avoiding the resource consumption caused by redundant calculations. Furthermore, all results calculated using the MatCloud-v2023 platform are securely uploaded to the cloud service, ensuring data integrity.

3. Results

3.1. Lattice Thermal Conductivity of Mg3Sb2

The lattice thermal conductivity of Mg3Sb2 was initially calculated without the application of external pressure. The results, obtained using various vdW corrections and directions, are presented in Figure 2a,b, respectively.
The values of lattice thermal conductivity obtained from PBE, PBE-D3, and vdW-DFq are quite similar. In general, vdW-DFq calculates slightly higher values for lattice thermal conductivity than PBE-D3, which, in turn, are slightly higher than those obtained from PBE. Since the lattice thermal conductivity is primarily governed by phonon transport, it typically increases with stronger atomic bonding. Therefore, the vdW corrections adds additional contributions, resulting in a slight increase in lattice thermal conductivity, which is understandable. Additionally, as the temperature rises, these differences tend to diminish. They are all very close to the experimental data reported in the literature [38]. Moving on to the lattice thermal conductivity in different directions, the variations are relatively small. For the Mg3Sb2 system, the values of the thermal conductivity tensor components (xx, yy, zz) exhibit minor differences, indicating limited anisotropy. The total thermal conductivity is typically the average of these three directions, which is used for further discussions on thermal conductivity. Finally, as the temperature increases, there is a decreasing trend in thermal conductivity. This is primarily attributed to the elevated temperature, which leads to increased phonon scattering rates and consequently reduces the lattice thermal conductivity.
We further examined the trends in lattice thermal conductivity under different pressures, and the results are illustrated in Figure 2c and listed in Table 1. Without external pressure, the thermal conductivity is at its maximum. The application of external pressure, whether compressive or tensile, leads to a decrease in thermal conductivity. Moreover, with increasing external pressure, there is an overall decreasing trend in thermal conductivity. However, this trend is not consistent throughout. When comparing the thermal conductivity at 3 GPa and 5 GPa, it increases at 5 GPa compared to 3 GPa. At a pressure above 7 GPa, the lattice thermal conductivity approaches zero.
It is noteworthy that, when the pressure reaches 10 GPa, there is a slight increase in thermal conductivity. Considering that the phonon spectrum of the crystal starts exhibiting significant imaginary frequencies at pressures exceeding 7 GPa, indicating structural dynamical instability, it is not necessary to further consider higher pressures.
We plotted the phonon spectra and phonon density of states of Mg3Sb2 under various pressures in Figure 3. The phonon spectra of Mg3Sb2 can be roughly divided into three regions. Below 3 THz, the primary contribution comes from the acoustic branches, mainly driven by the vibrations of Sb atoms, as observed in the phonon density of states. The range from 3 THz to 5 THz is dominated by the vibrations of magnesium atoms. In the absence of external pressure, there is no distinct boundary between these two parts. However, under tensile conditions, there is a tendency for these two parts to merge. Conversely, under compression, they gradually separate as the pressure increases. This phenomenon is one of the reasons for the reduction in thermal conductivity. The third part consists of optical branches with frequencies exceeding 5 THz. These optical phonons exhibit a noticeable gap compared to the acoustic phonons below. They are primarily generated by the vibrations of magnesium atoms. With an increase in pressure, the peak of the density of states increases, and phonon frequencies also increase. Furthermore, the gaps between different phonon branches increase as pressure increases.
To gain a deeper understanding of the contributions of different phonon modes to lattice thermal conductivity, we plotted the cumulative thermal conductivity as well as the three-phonon scattering rates at 300 K without pressure in Figure 4. The cumulative thermal conductivity is a graph where the horizontal axis represents the vibrational frequency, and it indicates the contribution of phonons with frequencies up to a certain value to the lattice thermal conductivity [73]. Naturally, this should be an increasing curve. For the results at 0 GPa, the curve can be roughly divided into three sections. There is an increasing process below 5 THz with a significant increase from 0 to 3 THz and a relatively smaller increase from 3 THz to 5 THz. Based on the phonon spectrum, these phonons correspond to the acoustic branches and a small fraction of the optical branches. In general, the vast majority of the contribution to thermal conductivity comes from the acoustic branches. In the range of 5–6 THz, the thermal conductivity increases level off. This corresponds to phonon modes with a state density of almost zero, indicating that there is no increase in thermal conductivity in this range. Afterwards, there is a slight increase, suggesting that optical phonons also contribute to the thermal conductivity of Mg3Sb2, increasing it by approximately 20%.
It is important to note that the anharmonic scattering process encompasses both absorption and emission processes when examining the scattering rates. The former involves two phonons combining to form a single phonon, while the latter is a process where one phonon splits into two. It can be observed that the phonon scattering rates undergo an initial increase followed by a decrease, peaking at around 4 THz. This peak corresponds to the initial phase of the cumulative thermal conductivity increase, where the rate of increase gradually decreases. Below 4 THz, the phonon scattering processes are dominated by phonon absorption, while in the range after 4 THz, the phonon emission processes predominate. Meanwhile, the gap between 5 THz and 6 THz correspond to the aforementioned portion of the phonon spectrum with a state density of nearly zero. This indicates that there are almost no scattering processes occurring in this frequency range.
Next, we examined the changes with increased pressure. We plotted these results in Figure 5. In the frequency range of 0–5 THz, both systems under added pressure, whether −2 GPa or 3 GPa, exhibit significantly higher phonon scattering rates compared to the scenario without any added pressure. This also results in significantly lower increases in the lattice thermal conductivity within this frequency range, compared to the scenario without pressure. The comparison of the phonon spectra suggests that the introduction of 3 GPa of pressure results in a gap at the 4 THz position. This causes the thermal conductivity at 3 GPa to enter the flat region earlier. On the other hand, introducing a pressure of −2 GPa significantly increases the scattering rate for phonon emission processes near 4 THz. This, in turn, reduces the contribution of these phonons to thermal conductivity.
Below 4 THz, the lattice thermal conductivity and scattering rate of the −2 GPa system are higher than those of the 3 GPa system. This is because the specific heat of the −2 GPa system is slightly higher than that of the 3 GPa system. Moving onto the latter part, which is the cumulative thermal conductivity beyond 5 THz, the 3 GPa system exhibits slightly higher scattering rates and, consequently, lower thermal conductivity. In any case, despite variations in the specific scattering mechanisms, the thermal conductivity ultimately converges to a similar value, with minimal disparity between them. For the 5 GPa scenario, the trends are identical to those of the 3 GPa system, with the only difference being that the former has relatively lower scattering rates, which results in higher lattice thermal conductivity.
Furthermore, we examined the scenario with an increased pressure of 7 GPa which is shown in Figure 6. From the phonon spectrum, it is evident that increased pressure causes a significant gap to open in the phonon spectrum near 4 THz. Although there are limited changes in the spectral lines below 4 THz, these alterations are significant enough to affect the overall thermal conductivity. Based solely on the phonon spectrum, it can be anticipated that, at 7 GPa, the thermal conductivity will be relatively low. When considering the scattering rates and cumulative thermal conductivity, it is observed that the scattering rates at 7 GPa are significantly higher, by several orders of magnitude, compared to those in the other systems. The excessively high scattering rates indicate that the phonons of the acoustic branch contribute very little to the lattice thermal conductivity in this scenario. It should be noted that, at 7 GPa, the Mg3Sb2 system starts to show some imaginary frequencies. As the pressure increases, these imaginary frequencies become more significant. This pressure may reach a critical point at which certain computational conditions are no longer met. The exact mechanisms involved in this behavior require further investigation.
In brief, applying pressure, regardless of its polarity, leads to a reduction in the lattice thermal conductivity of the Mg3Sb2 system. This observation suggests that increasing the pressure could potentially enhance the thermoelectric performance of the system to some extent. Moreover, when the pressure surpasses 7 GPa, the system exhibits the discernible presence of imaginary phonon frequencies, making the calculation of thermal conductivity unreliable beyond this threshold.

3.2. Thermoelectricity Properties of Ternary Mg3Bi2−vSbv

In terms of electronic structure, we utilized the PBE method to calculate the band structures of ternary Mg3Bi2−vSbv. The results obtained without applying pressure are shown in Figure 7. When comparing the band structures for different v values, we can observe that they are initially quite similar. During the transition from Mg3Bi2 to Mg3Sb2, the conduction band gradually shifts upward at the Γ point, eventually reaching a degeneracy with the preceding band. The bandgap also gradually increases. With the exception of the v = 0 and v = 0.5 structures, the other three structures all exhibit indirect bandgaps, with gap widths below 0.2 eV.
To evaluate the impact of increased pressure on the bandgaps of the structures, we subjected the system to various pressures and calculated the resulting bandgaps. For the purpose of comparison, we plotted the results as the difference between the CBM and VBM in Figure 7d. From the figure, it is evident that the band gaps for all systems initially increase and then decrease as the pressure varies. Moreover, this change, whether it is an increase or decrease in pressure, follows a linear trend. Although the v = 0 and v = 0.5 structures exhibit metallic properties under no pressure, they still develop bandgaps and demonstrate semiconductor characteristics when pressure is applied. The peak values of the band gaps differ for structures with a varying Sb content, but the overall trend indicates that the distance between the CBM and VBM increases with a higher Sb content. Empirically, narrow-gap semiconductors of this type tend to exhibit favorable thermoelectric properties.
We used TransOpt-v2.0 to calculate the electrical transport properties of the Mg3Bi2−vSbv system. This software can simultaneously output results based on both the CRTA and the CEPCA approaches. When comparing with experimental data, we observed a significant discrepancy. Specifically, the electrical conductivity obtained through the CEPCA was higher by two to three orders of magnitude compared to the experimental results [38]. This discrepancy indicates that considering the electron–phonon coupling matrix as a constant does not accurately represent the electron scattering scenarios in the Mg3Bi2−vSbv system. For the sake of simplicity, we initially disregard the relaxation time. Figure 8 illustrates the changes in the power factor, the absolute value of the Seebeck coefficient, and electrical conductivity as a function of pressure for the Mg3Bi2−vSbv system under n-type doping, while maintaining a constant carrier concentration of 0.03 ×   10 20   cm 3 . Combining the bandgap information from Figure 7d, it is evident that, for the Mg3BiSb (v = 1), Mg3Bi0.5Sb1.5 (v = 1.5), and Mg3Sb2 (v = 2) structures, as per Equations (1) and (2), the bandgap exhibits an initial increase followed by a decrease with pressure. Consequently, there is an initial decrease followed by an increase in electrical conductivity, while the Seebeck coefficient exhibits an initial increase followed by a decrease. This leads to the power factor reaching its maximum value at moderate pressures. However, the observed trend is not as evident for Mg3Bi2 (v = 0) and Mg3Bi1.5Sb0.5 (v = 0.5), which may be attributed to the initially small bandgap, leading to less pronounced changes. Further exploration is necessary to elucidate a more detailed mechanism.
We plotted the variation of the maximum power factor for both p-type and n-type carriers with pressure in Figure 9 and Figure 10, respectively.
For p-type Mg3Bi2−vSbv, an overall trend is that the peak power factor of all systems increases as v increases, except for the structure with v = 0.5. With increasing pressure, the peak values of structures with v = 0, 1.0, and 1.5 increase, while the other two structures exhibit varying trends of increase and decrease, similar to the Seebeck coefficient. It is evident that, despite the relatively low Seebeck coefficient for v = 0.5, the higher electrical conductivity results in a reasonably significant power factor. This suggests that, at higher doping levels, the system with v = 0.5 may exhibit favorable thermoelectric performance. For the Mg3Bi2 structure, a significant drop in electrical conductivity is observed at a pressure of 7 GPa. This is because the peak power factor at this point shifts the carrier concentration from around 5.0 ×   10 20   cm 3 to 1.5 ×   10 20   cm 3 , indicating a decrease in the required doping concentration for optimal performance.
For n-type Mg3Bi2−vSbv, the power factor increases with the increase in Sb content. Additionally, as pressure increases, the peak power factors of all systems also increase. Furthermore, in n-type semiconductors, the peak power factors are generally higher than those in p-type semiconductors. This indicates that n-type doping tends to result in higher thermoelectric performance when the system is doped. However, the Seebeck coefficient and electrical conductivity do not exhibit a distinct trend with increasing pressure. This behavior is related to the different positions of the optimal carrier concentration.
In order to obtain an accurate calculation of the zT value for the system, it is necessary to take into account the relaxation time. During the calculation, the relaxation time is considered to be an unknown constant and is determined by comparing it with experimental values. We adopt the simplifying assumption that the relaxation time remains constant regardless of changes in doping concentration and pressure. When comparing the experimental results obtained at a doping concentration of 0.39 ×   10 20   cm 3 and a temperature of 300 K, the Seebeck coefficient calculated using the CRTA method in TransOpt is 122 μV/K. This value closely aligns with the experimental value of 114 μV/K [38]. Through this comparative analysis, it is observed that the relaxation time is estimated to be approximately 4.55 femtoseconds. We will utilize this period of relaxation to engage in a discussion regarding the computation of the zT value for Mg3Sb2.
The zT values of Mg3Sb2 were plotted as a function of pressure in Figure 8c and Figure 9d, incorporating the thermal conductivity previously calculated for Mg3Sb2. Some data are listed in Table 2. Figure 8c illustrates the p-type semiconductor, whereas Figure 9d depicts the n-type doping scenario. The zT values reported in this study were determined based on the optimal power factor. Higher temperatures lead to higher zT values at the same pressure. For p-type semiconductors, the zT value is initially at its lowest without the application of additional pressure. However, increasing the pressure up to a certain threshold enhances the peak value of zT. Beyond this threshold, further increases in pressure lead to a decrease in the zT value. At a pressure of 3 GPa and a temperature of 800 K, the maximum achieved zT value is 0.55, accompanied by a carrier concentration of 1.99 ×   10 20   cm 3 . For n-type semiconductors, the zT exhibits an increasing trend with rising pressure, reaching its maximum value at 3 GPa, after which it starts to decline. At a temperature of 800 K, the highest achieved zT value is 1.53, accompanied by a carrier concentration of 1.99 ×   10 20   cm 3 .
In summary, the thermoelectric performance of n-type doping in the Mg3Bi2−vSbv system is generally superior to that of p-type doping. Under pressure modulation, the power factor exhibits fluctuations. For the Mg3Sb2 system, the maximum zT value is achieved at a pressure of 3 GPa and a temperature of 800 K, reaching a maximum of 1.53. This result of the unrecorded high zT value of Mg3Sb2 might suggest that the zT value of Mg3Bi2−vSbv systems can be further improved by synergistically tuning the temperature, pressure, and carrier concentration.

4. Conclusions

We conducted a comprehensive investigation to examine the influence of pressure on the thermoelectric properties of Mg3Bi2−vSbv. We used two computational methods, PBE-D3 and vdW-DFq, to include non-local vdW corrections in calculating the lattice thermal conductivity of Mg3Sb2. The inclusion of vdW effects generally led to slightly larger outcomes compared to simulations that did not take vdW corrections into account. However, it is important to note that the observed difference was not statistically significant.
We calculated the lattice thermal conductivity of Mg3Sb2 at different pressures and temperatures and conducted an analysis. The thermal conductivity of the lattice in this structure decreases as the temperature increases. Additionally, the value decreases as the pressure increases, eventually approaching zero at 7 GPa. We conducted a thorough analysis of the contributions made by different phonon modes to the overall lattice thermal conductivity. We then discussed the effect of pressure on the electronic transport properties of Mg3Bi2−vSbv. Both Mg3Sb2 and Mg3Bi1.5Sb0.5 showed significantly higher power factors compared to the other three systems when used for n-type doping. No distinct pattern in their behavior was observed when different levels of pressure were applied. However, the other three systems showed improvements in power factor.
The power factor increased with a higher Sb content and pressure for n-type doping, surpassing that of p-type doping. The zT values of Mg3Sb2 were calculated at different pressures and temperatures. Regardless of the type of doping, whether it is n-type or p-type, the maximum values of the zT were achieved at a temperature of 800 K and a pressure of 3 GPa. A dopant concentration of 1.99 × 10 20   cm 3 was found to be necessary to achieve the maximum zT values. The highest achieved zT value for n-type doping was 1.53, while for p-type doping, it reached 0.55. Our research uncovered the performance characteristics of Mg3Bi2−vSbv under varying external pressure conditions, offering valuable insights into the potential applications of this system in the realm of thermoelectric materials and additional revenue to further improve zT by synergistically tuning the temperature, pressure, and carrier concentration.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano14010084/s1, Figure S1: Workflow for calculating third-order constants on MatCloud.

Author Contributions

Q.P. and X.-J.C. conceived the project; X.M. performed the first-principle-based calculation and analysis. X.M. and Q.P. share equal contributions. X.Y. (Xiaoyu Yang) provides software support. X.Y. (Xiaoze Yuan) provides structural configuration. All authors analyzed data, discussed the results, wrote the draft, and reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Shenzhen Science and Technology Program with grant number KQTD20200820113045081, National Natural Science Foundation of China (Grant No. 12272378), High-level Innovation Research Institute Program of Guangdong Province (Grant No. 2020B0909010003), and LiYing Program of the Institute of Mechanics, Chinese Academy of Sciences (Grant No. E1Z1011001).

Data Availability Statement

Data is contained within the article (and Supplementary Materials).

Acknowledgments

The authors also express their gratitude to Beijing MaiGao MatCloud Technology Co., Ltd. for providing computational resources and the computing platform.

Conflicts of Interest

Authors Qing Peng, Xinjie Ma and Xiaoyu Yang are employed by the company Beijing MaiGao MatCloud Technology Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Shi, X.-L.; Zou, J.; Chen, Z.-G. Advanced Thermoelectric Design: From Materials and Structures to Devices. Chem. Rev. 2020, 120, 7399–7515. [Google Scholar] [CrossRef] [PubMed]
  2. Tan, G.; Zhao, L.-D.; Kanatzidis, M.G. Rationally Designing High-Performance Bulk Thermoelectric Materials. Chem. Rev. 2016, 116, 12123–12149. [Google Scholar] [CrossRef] [PubMed]
  3. Shi, X.; Chen, L. Thermoelectric materials step up. Nat. Mater. 2016, 15, 691–692. [Google Scholar] [CrossRef] [PubMed]
  4. DiSalvo, F.J. Thermoelectric Cooling and Power Generation. Science 1999, 285, 703–706. [Google Scholar] [CrossRef] [PubMed]
  5. Bell, L.E. Cooling, Heating, Generating Power, and Recovering Waste Heat with Thermoelectric Systems. Science 2008, 321, 1457–1461. [Google Scholar] [CrossRef] [PubMed]
  6. Qiu, P.; Shi, X.; Chen, L. Cu-based thermoelectric materials. Energy Storage Mater. 2016, 3, 85–97. [Google Scholar] [CrossRef]
  7. Mehdizadeh Dehkordi, A.; Zebarjadi, M.; He, J.; Tritt, T.M. Thermoelectric power factor: Enhancement mechanisms and strategies for higher performance thermoelectric materials. Mater. Sci. Eng. R Rep. 2015, 97, 1–22. [Google Scholar] [CrossRef]
  8. Moshwan, R.; Yang, L.; Zou, J.; Chen, Z.-G. Eco-Friendly SnTe Thermoelectric Materials: Progress and Future Challenges. Adv. Funct. Mater. 2017, 27, 1703278. [Google Scholar] [CrossRef]
  9. Zhang, X.; Zhao, L.-D. Thermoelectric materials: Energy conversion between heat and electricity. J. Mater. 2015, 1, 92–105. [Google Scholar] [CrossRef]
  10. Liu, W.; Yan, X.; Chen, G.; Ren, Z. Recent advances in thermoelectric nanocomposites. Nano Energy 2012, 1, 42–56. [Google Scholar] [CrossRef]
  11. Minnich, A.J.; Dresselhaus, M.S.; Ren, Z.F.; Chen, G. Bulk nanostructured thermoelectric materials: Current research and future prospects. Energy Environ. Sci. 2009, 2, 466–479. [Google Scholar] [CrossRef]
  12. Tripathi, N.; Pavelyev, V.; Sharma, P.; Kumar, S.; Rymzhina, A.; Mishra, P. Review of titanium trisulfide (TiS3): A novel material for next generation electronic and optical devices. Mater. Sci. Semicond. Process. 2021, 127, 105699. [Google Scholar] [CrossRef]
  13. Kim, S.; Lee, C.; Lim, Y.S.; Shim, J.-H. Investigation for Thermoelectric Properties of the MoS2 Monolayer–Graphene Heterostructure: Density Functional Theory Calculations and Electrical Transport Measurements. ACS Omega 2021, 6, 278–283. [Google Scholar] [CrossRef] [PubMed]
  14. Kumanek, B.; Stando, G.; Stando, P.; Matuszek, K.; Milowska, K.Z.; Krzywiecki, M.; Gryglas-Borysiewicz, M.; Ogorzałek, Z.; Payne, M.C.; MacFarlane, D.; et al. Enhancing thermoelectric properties of single-walled carbon nanotubes using halide compounds at room temperature and above. Sci. Rep. 2021, 11, 8649. [Google Scholar] [CrossRef]
  15. Gaul, A.; Peng, Q.; Singh, D.J.; Ramanath, G.; Borca-Tasciuc, T. Pressure-induced insulator-to-metal transitions for enhancing thermoelectric power factor in bismuth telluride-based alloys. Phys. Chem. Chem. Phys. 2017, 19, 12784–12793. [Google Scholar] [CrossRef] [PubMed]
  16. Yang, G.; Cui, H.; Ma, D.; He, C. The elastic and thermoelectric properties of the Zintl compound Ca5Al2Sb6 under high pressure. J. Appl. Phys. 2014, 116, 223709. [Google Scholar] [CrossRef]
  17. Hao, S.; Dravid, V.P.; Kanatzidis, M.G.; Wolverton, C. Research Update: Prediction of high figure of merit plateau in SnS and solid solution of (Pb,Sn)S. APL Mater. 2016, 4, 104505. [Google Scholar] [CrossRef]
  18. Hao, S.; Shi, F.; Dravid, V.P.; Kanatzidis, M.G.; Wolverton, C. Computational Prediction of High Thermoelectric Performance in Hole Doped Layered GeSe. Chem. Mater. 2016, 28, 3218–3226. [Google Scholar] [CrossRef]
  19. Ma, Y.; Eremets, M.; Oganov, A.R.; Xie, Y.; Trojan, I.; Medvedev, S.; Lyakhov, A.O.; Valle, M.; Prakapenka, V. Transparent dense sodium. Nature 2009, 458, 182–185. [Google Scholar] [CrossRef]
  20. Matsuoka, T.; Shimizu, K. Direct observation of a pressure-induced metal-to-semiconductor transition in lithium. Nature 2009, 458, 186–189. [Google Scholar] [CrossRef]
  21. Seeger, K. Semiconductor Physics; Springer: Wien, Austria, 1973. [Google Scholar]
  22. Morozova, N.V.; Korobeinikov, I.V.; Ovsyannikov, S.V. Strategies and challenges of high-pressure methods applied to thermoelectric materials. J. Appl. Phys. 2019, 125, 220901. [Google Scholar] [CrossRef]
  23. Besson, J.M.; Paul, W.; Calawa, A.R. Tuning of PbSe Lasers by Hydrostatic Pressure from 8 to 22 μ. Phys. Rev. 1968, 173, 699–713. [Google Scholar] [CrossRef]
  24. Nabi, Z.; Abbar, B.; Méçabih, S.; Khalfi, A.; Amrane, N. Pressure dependence of band gaps in PbS, PbSe and PbTe. Comput. Mater. Sci. 2000, 18, 127–131. [Google Scholar] [CrossRef]
  25. Li, C.-Y.; Ruoff, A.L.; Spencer, C.W. Effect of Pressure on the Energy Gap of Bi2Te3. J. Appl. Phys. 1961, 32, 1733–1735. [Google Scholar] [CrossRef]
  26. Peng, Q.; Yuan, X.; Zhao, S.; Chen, X.-J. Lattice Thermal Conductivity of Mg3(Bi,Sb)2 Nanocomposites: A First-Principles Study. Nanomaterials 2023, 13, 2938. [Google Scholar] [CrossRef] [PubMed]
  27. Kanno, T.; Tamaki, H.; Yoshiya, M.; Uchiyama, H.; Maki, S.; Takata, M.; Miyazaki, Y. High-Density Frenkel Defects as Origin of N-Type Thermoelectric Performance and Low Thermal Conductivity in Mg3Sb2-Based Materials. Adv. Funct. Mater. 2021, 31, 2008469. [Google Scholar] [CrossRef]
  28. Imasato, K.; Kang, S.D.; Ohno, S.; Snyder, G.J. Band engineering in Mg3Sb2 by alloying with Mg3Bi2 for enhanced thermoelectric performance. Mater. Horiz. 2018, 5, 59–64. [Google Scholar] [CrossRef]
  29. Pan, Y.; Yao, M.; Hong, X.; Zhu, Y.; Fan, F.; Imasato, K.; He, Y.; Hess, C.; Fink, J.; Yang, J.; et al. Mg3(Bi,Sb)2 single crystals towards high thermoelectric performance. Energy Environ. Sci. 2020, 13, 1717–1724. [Google Scholar] [CrossRef]
  30. Shi, X.; Sun, C.; Bu, Z.; Zhang, X.; Wu, Y.; Lin, S.; Li, W.; Faghaninia, A.; Jain, A.; Pei, Y. Revelation of Inherently High Mobility Enables Mg3Sb2 as a Sustainable Alternative to n-Bi2Te3 Thermoelectrics. Adv. Sci. 2019, 6, 1802286. [Google Scholar] [CrossRef]
  31. Zhang, F.; Chen, C.; Yao, H.; Bai, F.; Yin, L.; Li, X.; Li, S.; Xue, W.; Wang, Y.; Cao, F.; et al. High-Performance N-type Mg3Sb2 towards Thermoelectric Application near Room Temperature. Adv. Funct. Mater. 2020, 30, 1906143. [Google Scholar] [CrossRef]
  32. Zhang, J.; Song, L.; Pedersen, S.H.; Yin, H.; Hung, L.T.; Iversen, B.B. Discovery of high-performance low-cost n-type Mg3Sb2-based thermoelectric materials with multi-valley conduction bands. Nat. Commun. 2017, 8, 13901. [Google Scholar] [CrossRef] [PubMed]
  33. Imasato, K.; Kang, S.D.; Snyder, G.J. Exceptional thermoelectric performance in Mg3Sb0.6Bi1.4 for low-grade waste heat recovery. Energy Environ. Sci. 2019, 12, 965–971. [Google Scholar] [CrossRef]
  34. Tamaki, H.; Sato, H.K.; Kanno, T. Isotropic Conduction Network and Defect Chemistry in Mg3+δSb2-Based Layered Zintl Compounds with High Thermoelectric Performance. Adv. Mater. 2016, 28, 10182–10187. [Google Scholar] [CrossRef] [PubMed]
  35. Imasato, K.; Wood, M.; Kuo, J.J.; Snyder, G.J. Improved stability and high thermoelectric performance through cation site doping in n-type La-doped Mg3Sb1.5Bi0.5. J. Mater. Chem. A 2018, 6, 19941–19946. [Google Scholar] [CrossRef]
  36. Wood, M.; Kuo, J.J.; Imasato, K.; Snyder, G.J. Improvement of Low-Temperature zT in a Mg3Sb2–Mg3Bi2 Solid Solution via Mg-Vapor Annealing. Adv. Mater. 2019, 31, 1902337. [Google Scholar] [CrossRef]
  37. Zhang, J.; Song, L.; Iversen, B.B. Insights into the design of thermoelectric Mg3Sb2 and its analogs by combining theory and experiment. NPJ Comput. Mater. 2019, 5, 76. [Google Scholar] [CrossRef]
  38. Shuai, J.; Wang, Y.; Kim, H.S.; Liu, Z.; Sun, J.; Chen, S.; Sui, J.; Ren, Z. Thermoelectric properties of Na-doped Zintl compound: Mg3−xNaxSb2. Acta Mater. 2015, 93, 187–193. [Google Scholar] [CrossRef]
  39. Ao, D.-W.; Liu, W.-D.; Chen, Y.-X.; Wei, M.; Jabar, B.; Li, F.; Shi, X.-L.; Zheng, Z.-H.; Liang, G.-X.; Zhang, X.-H.; et al. Novel Thermal Diffusion Temperature Engineering Leading to High Thermoelectric Performance in Bi2Te3-Based Flexible Thin-Films. Adv. Sci. 2022, 9, 2103547. [Google Scholar] [CrossRef]
  40. Yang, Q.; Yang, S.; Qiu, P.; Peng, L.; Wei, T.-R.; Zhang, Z.; Shi, X.; Chen, L. Flexible thermoelectrics based on ductile semiconductors. Science 2022, 377, 854–858. [Google Scholar] [CrossRef]
  41. Bayikadi, K.S.; Wu, C.T.; Chen, L.-C.; Chen, K.-H.; Chou, F.-C.; Sankar, R. Synergistic optimization of thermoelectric performance of Sb doped GeTe with a strained domain and domain boundaries. J. Mater. Chem. A 2020, 8, 5332–5341. [Google Scholar] [CrossRef]
  42. Yu, R.; Fang, Z.; Dai, X.; Weng, H. Topological nodal line semimetals predicted from first-principles calculations. Front. Phys. 2017, 12, 127202. [Google Scholar] [CrossRef]
  43. Teshome, T.; Datta, A. Topological Phase Transition in Sb2Mg3 Assisted by Strain. ACS Omega 2019, 4, 8701–8706. [Google Scholar] [CrossRef] [PubMed]
  44. Li, S.; Yu, Z.-M.; Liu, Y.; Guan, S.; Wang, S.-S.; Zhang, X.; Yao, Y.; Yang, S.A. Type-II nodal loops: Theory and material realization. Phys. Rev. B 2017, 96, 081106. [Google Scholar] [CrossRef]
  45. Chen, W.; Lu, H.-Z.; Zilberberg, O. Weak Localization and Antilocalization in Nodal-Line Semimetals: Dimensionality and Topological Effects. Phys. Rev. Lett. 2019, 122, 196603. [Google Scholar] [CrossRef] [PubMed]
  46. Zhu, Y.; Liu, J.; Wei, B.; Xu, S.; Song, Y.; Wang, X.; Xia, T.-L.; Chen, J.; Snyder, G.J.; Hong, J. Giant phonon anharmonicity driven by the asymmetric lone pairs in Mg3Bi2. Mater. Today Phys. 2022, 27, 100791. [Google Scholar] [CrossRef]
  47. Xin, J.; Li, G.; Auffermann, G.; Borrmann, H.; Schnelle, W.; Gooth, J.; Zhao, X.; Zhu, T.; Felser, C.; Fu, C. Growth and transport properties of Mg3X2 (X = Sb, Bi) single crystals. Mater. Today Phys. 2018, 7, 61–68. [Google Scholar] [CrossRef]
  48. Agne, M.T.; Imasato, K.; Anand, S.; Lee, K.; Bux, S.K.; Zevalkink, A.; Rettie, A.J.E.; Chung, D.Y.; Kanatzidis, M.G.; Snyder, G.J. Heat capacity of Mg3Sb2, Mg3Bi2, and their alloys at high temperature. Mater. Today Phys. 2018, 6, 83–88. [Google Scholar] [CrossRef]
  49. Ding, J.; Lanigan-Atkins, T.; Calderón-Cueva, M.; Banerjee, A.; Abernathy, D.L.; Said, A.; Zevalkink, A.; Delaire, O. Soft anharmonic phonons and ultralow thermal conductivity in Mg3(Sb,Bi)2 thermoelectrics. Sci. Adv. 2021, 7, eabg1449. [Google Scholar] [CrossRef]
  50. Peng, W.; Petretto, G.; Rignanese, G.-M.; Hautier, G.; Zevalkink, A. An Unlikely Route to Low Lattice Thermal Conductivity: Small Atoms in a Simple Layered Structure. Joule 2018, 2, 1879–1893. [Google Scholar] [CrossRef]
  51. Peng, Q.; Zhao, S.; Yuan, X.; Chen, X.-J. Elasticity of Mg3Bi2-xSbx. Materials 2022, 15, 7161. [Google Scholar] [CrossRef]
  52. Peng, Q.; Ma, X.; Yang, X.; Zhao, S.; Yuan, X.; Chen, X. Assessing Effects of van der Waals Corrections on Elasticity of Mg3Bi2−xSbx in DFT Calculations. Materials 2023, 16, 6482. [Google Scholar] [CrossRef]
  53. Gooth, J.; Schierning, G.; Felser, C.; Nielsch, K. Quantum materials for thermoelectricity. MRS Bull. 2018, 43, 187–192. [Google Scholar] [CrossRef]
  54. Li, M.; Chen, G. Thermal transport for probing quantum materials. MRS Bull. 2020, 45, 348–356. [Google Scholar] [CrossRef]
  55. Barati, S.; Abedinpour, S.H. Thermoelectric response of nodal-line semimetals: Probing the Fermi surface topology. Phys. Rev. B 2020, 102, 125139. [Google Scholar] [CrossRef]
  56. Hosoi, M.; Tateishi, I.; Matsuura, H.; Ogata, M. Thin films of topological nodal line semimetals as a candidate for efficient thermoelectric converters. Phys. Rev. B 2022, 105, 085406. [Google Scholar] [CrossRef]
  57. Yuan, X.; Zhou, Y.; Peng, Q.; Yang, Y.; Li, Y.; Wen, X. Active learning to overcome exponential-wall problem for effective structure prediction of chemical-disordered materials. NPJ Comput. Mater. 2023, 9, 12. [Google Scholar] [CrossRef]
  58. Peng, Q.; Yuan, X.; Zhao, S.; Zhou, Y.; Wen, X.; Chen, X. Active-learning search for unitcell structures: A case study on Mg3Bi2−xSbx. Comput. Mater. Sci. 2023, 226, 112260. [Google Scholar] [CrossRef]
  59. Kohn, W.; Sham, L.J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133–A1138. [Google Scholar] [CrossRef]
  60. Kresse, G.; Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  61. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef]
  62. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  63. Wang, V.; Xu, N.; Liu, J.-C.; Tang, G.; Geng, W.-T. VASPKIT: A user-friendly interface facilitating high-throughput computing and analysis using VASP code. Comput. Phys. Commun. 2021, 267, 108033. [Google Scholar] [CrossRef]
  64. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  65. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef]
  66. Peng, Q.; Wang, G.; Liu, G.-R.; De, S. Van der Waals Density Functional Theory vdW-DFq for Semihard Materials. Crystals 2019, 9, 243. [Google Scholar] [CrossRef]
  67. Allen, P.B. Boltzmann Theory and Resistivity of Metals. In Quantum Theory of Real Materials; Chelikowsky, J.R., Louie, S.G., Eds.; The Kluwer International Series in Engineering and Computer Science; Springer: Boston, MA, USA, 1996; pp. 219–250. [Google Scholar] [CrossRef]
  68. Li, X.; Zhang, Z.; Xi, J.; Singh, D.J.; Sheng, Y.; Yang, J.; Zhang, W. TransOpt. A code to solve electrical transport properties of semiconductors in constant electron–phonon coupling approximation. Comput. Mater. Sci. 2021, 186, 110074. [Google Scholar] [CrossRef]
  69. Togo, A.; Tanaka, I. First principles phonon calculations in materials science. Scr. Mater. 2015, 108, 1–5. [Google Scholar] [CrossRef]
  70. Li, W.; Carrete, J.A.; Katcho, N.; Mingo, N. ShengBTE: A solver of the Boltzmann transport equation for phonons. Comput. Phys. Commun. 2014, 185, 1747–1758. [Google Scholar] [CrossRef]
  71. Han, Z.; Yang, X.; Li, W.; Feng, T.; Ruan, X. FourPhonon: An extension module to ShengBTE for computing four-phonon scattering rates and thermal conductivity. Comput. Phys. Commun. 2022, 270, 108179. [Google Scholar] [CrossRef]
  72. Yang, X.; Wang, Z.; Zhao, X.; Song, J.; Zhang, M.; Liu, H. MatCloud: A high-throughput computational infrastructure for integrated management of materials simulation, data and resources. Comput. Mater. Sci. 2018, 146, 319–333. [Google Scholar] [CrossRef]
  73. Li, W.; Mingo, N. Ultralow lattice thermal conductivity of the fully filled skutterudite YbFe4Sb12 due to the flat avoided-crossing filler modes. Phys. Rev. B 2015, 91, 144304. [Google Scholar] [CrossRef]
Figure 1. The atomistic structures of the Mg3Bi2−vSbv: the 5-atom unit cell of (a) Mg3Bi2 and (b) Mg3Sb2, and (ce) the 10-atom cells of Mg3Bi2−vSbv with v in 0.5, 1.0, and 1.5, respectively.
Figure 1. The atomistic structures of the Mg3Bi2−vSbv: the 5-atom unit cell of (a) Mg3Bi2 and (b) Mg3Sb2, and (ce) the 10-atom cells of Mg3Bi2−vSbv with v in 0.5, 1.0, and 1.5, respectively.
Nanomaterials 14 00084 g001
Figure 2. (a) Lattice thermal conductivity of Mg3Sb2 at varying temperatures in various vdW and experimental data [38]. (b) Lattice thermal conductivity of Mg3Sb2 at varying temperatures in various directions. (c) Lattice thermal conductivity under different pressures and at different temperatures.
Figure 2. (a) Lattice thermal conductivity of Mg3Sb2 at varying temperatures in various vdW and experimental data [38]. (b) Lattice thermal conductivity of Mg3Sb2 at varying temperatures in various directions. (c) Lattice thermal conductivity under different pressures and at different temperatures.
Nanomaterials 14 00084 g002
Figure 3. Phonon spectra and phonon density of states of Mg3Sb2 at various pressures.
Figure 3. Phonon spectra and phonon density of states of Mg3Sb2 at various pressures.
Nanomaterials 14 00084 g003
Figure 4. (a) Cumulative lattice thermal conductivity of Mg3Sb2 at 300 K without pressure. (b) Anharmonic scattering rate of Mg3Sb2 at 300 K without pressure.
Figure 4. (a) Cumulative lattice thermal conductivity of Mg3Sb2 at 300 K without pressure. (b) Anharmonic scattering rate of Mg3Sb2 at 300 K without pressure.
Nanomaterials 14 00084 g004
Figure 5. (a) Cumulative lattice thermal conductivity of Mg3Sb2 at 300 K under various pressures. (bd) Anharmonic scattering rate of Mg3Sb2 at 300 K under various pressures.
Figure 5. (a) Cumulative lattice thermal conductivity of Mg3Sb2 at 300 K under various pressures. (bd) Anharmonic scattering rate of Mg3Sb2 at 300 K under various pressures.
Nanomaterials 14 00084 g005
Figure 6. (a) Cumulative lattice thermal conductivity of Mg3Sb2 at 300 K without pressure. (b) Anharmonic scattering rate of Mg3Sb2 at 300 K under 7 GPa.
Figure 6. (a) Cumulative lattice thermal conductivity of Mg3Sb2 at 300 K without pressure. (b) Anharmonic scattering rate of Mg3Sb2 at 300 K under 7 GPa.
Nanomaterials 14 00084 g006
Figure 7. (a) Band structure of Mg3Bi2. (b) Band structure of Mg3Sb2. (c) Band structure of Mg3Bi2−vSbv for v = 0.5, v = 1.0, and v = 1.5. (d) Bandgap (CBM-VBM) of Mg3Bi2−vSbv under various pressures.
Figure 7. (a) Band structure of Mg3Bi2. (b) Band structure of Mg3Sb2. (c) Band structure of Mg3Bi2−vSbv for v = 0.5, v = 1.0, and v = 1.5. (d) Bandgap (CBM-VBM) of Mg3Bi2−vSbv under various pressures.
Nanomaterials 14 00084 g007
Figure 8. (a) The changes in power factor as a function of pressure for the Mg3Bi2−vSbv system under n-type doping, maintaining a constant carrier concentration of 0.03 ×   10 20   cm 3 . (b) The absolute value of the Seebeck coefficient as a function of pressure for the Mg3Bi2−vSbv system under n-type doping, maintaining a constant carrier concentration of 0.03 ×   10 20   cm 3 . (c) The electrical conductivity as a function of pressure for the Mg3Bi2−vSbv system under n-type doping, maintaining a constant carrier concentration of 0.03 ×   10 20   cm 3 .
Figure 8. (a) The changes in power factor as a function of pressure for the Mg3Bi2−vSbv system under n-type doping, maintaining a constant carrier concentration of 0.03 ×   10 20   cm 3 . (b) The absolute value of the Seebeck coefficient as a function of pressure for the Mg3Bi2−vSbv system under n-type doping, maintaining a constant carrier concentration of 0.03 ×   10 20   cm 3 . (c) The electrical conductivity as a function of pressure for the Mg3Bi2−vSbv system under n-type doping, maintaining a constant carrier concentration of 0.03 ×   10 20   cm 3 .
Nanomaterials 14 00084 g008
Figure 9. (a) The maximum PF/τ of p-type-doped Mg3Bi2−vSbv under various pressures. (b) Seeback coefficient of p-type-doped Mg3Bi2−vSbv under various pressures. (c) σ /τ of p-type-doped Mg3Bi2−vSbv under various pressures. (d) Max zT of p-type doped Mg3Sb2 under various pressures.
Figure 9. (a) The maximum PF/τ of p-type-doped Mg3Bi2−vSbv under various pressures. (b) Seeback coefficient of p-type-doped Mg3Bi2−vSbv under various pressures. (c) σ /τ of p-type-doped Mg3Bi2−vSbv under various pressures. (d) Max zT of p-type doped Mg3Sb2 under various pressures.
Nanomaterials 14 00084 g009
Figure 10. (a) The maximum PF/τ of n-type-doped Mg3Bi2−vSbv under various pressures. (b) Seeback coefficient of n-type-doped Mg3Bi2−vSbv under various pressures. (c) σ /τ of n-type-doped Mg3Bi2−vSbv under various pressures. (d) Max zT of n-type-doped Mg3Sb2 under various pressures.
Figure 10. (a) The maximum PF/τ of n-type-doped Mg3Bi2−vSbv under various pressures. (b) Seeback coefficient of n-type-doped Mg3Bi2−vSbv under various pressures. (c) σ /τ of n-type-doped Mg3Bi2−vSbv under various pressures. (d) Max zT of n-type-doped Mg3Sb2 under various pressures.
Nanomaterials 14 00084 g010
Table 1. This table shows the lattice thermal conductivity of Mg3Sb2 at different temperatures and different pressures.
Table 1. This table shows the lattice thermal conductivity of Mg3Sb2 at different temperatures and different pressures.
Temperature (K)Pressure (GPa)Lattice Thermal Conductivity (W/m·K)
300−2.00.573
0.01.251
1.00.988
5.00.519
500−2.00.258
0.00.765
1.00.605
5.00.322
800−2.00.165
0.00.483
1.00.383
5.00.204
Table 2. This table list max zT of Mg3Sb2 at different temperatures under 3 GPa.
Table 2. This table list max zT of Mg3Sb2 at different temperatures under 3 GPa.
Main Carrier TypeTemperature (K)Concentration
(×1020 cm−3)
zT
Holes (p-type)3000.850.17
5000.890.43
8001.990.55
Electrons (n-type)3002.000.21
5002.000.74
8001.991.53
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Peng, Q.; Ma, X.; Yang, X.; Yuan, X.; Chen, X.-J. Thermoelectric Properties of Mg3(Bi,Sb)2 under Finite Temperatures and Pressures: A First-Principles Study. Nanomaterials 2024, 14, 84. https://doi.org/10.3390/nano14010084

AMA Style

Peng Q, Ma X, Yang X, Yuan X, Chen X-J. Thermoelectric Properties of Mg3(Bi,Sb)2 under Finite Temperatures and Pressures: A First-Principles Study. Nanomaterials. 2024; 14(1):84. https://doi.org/10.3390/nano14010084

Chicago/Turabian Style

Peng, Qing, Xinjie Ma, Xiaoyu Yang, Xiaoze Yuan, and Xiao-Jia Chen. 2024. "Thermoelectric Properties of Mg3(Bi,Sb)2 under Finite Temperatures and Pressures: A First-Principles Study" Nanomaterials 14, no. 1: 84. https://doi.org/10.3390/nano14010084

APA Style

Peng, Q., Ma, X., Yang, X., Yuan, X., & Chen, X. -J. (2024). Thermoelectric Properties of Mg3(Bi,Sb)2 under Finite Temperatures and Pressures: A First-Principles Study. Nanomaterials, 14(1), 84. https://doi.org/10.3390/nano14010084

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop