Next Article in Journal
Effect of Mischmetal Addition on Microstructure and Mechanical Properties of As-Cast and As-Rolled Mg–Sn–Ca Alloys
Next Article in Special Issue
Screening of Complex Layered Chalcogenide Structures as High-Performance Thermoelectrics by High-Throughput Calculations
Previous Article in Journal
Particle Swarm Method for Optimization of ATIG Welding Process to Joint Mild Steel to 316L Stainless Steel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Throughput Exploration of Half-Heusler Phases for Thermoelectric Applications

by
Kaja Bilińska
and
Maciej J. Winiarski
*
Institute of Low Temperature and Structure Research, Polish Academy of Sciences, Okólna 2, 50-370 Wrocław, Poland
*
Author to whom correspondence should be addressed.
Crystals 2023, 13(9), 1378; https://doi.org/10.3390/cryst13091378
Submission received: 30 August 2023 / Revised: 13 September 2023 / Accepted: 15 September 2023 / Published: 17 September 2023

Abstract

:
As a result of the high-throughput ab initiocalculations, the set of 34 stable and novel half-Heusler phases was revealed. The electronic structure and the elastic, transport, and thermoelectric properties of these systems were carefully investigated, providing some promising candidates for thermoelectric materials. The complementary nature of the research is enhanced by the deformation potential theory applied for the relaxation time of carriers (for power factor, PF) and the Slack formula for the lattice thermal conductivity (for figure of merit, ZT). Moreover, two exchange-correlation parametrizations were used (GGA and MBJGGA), and a complete investigation was provided for both p- and n-type carriers. The distribution of the maximum PF and ZT for optimal doping at 300 K in all systems was disclosed. Some chemical trends in electronic and transport properties were discussed. The results suggest TaFeAs, TaFeSb, VFeAs, and TiRuAs as potentially valuable thermoelectric materials. TaFeAs revealed the highest values of both PF and ZT at 300 K (PF p = 1.67 mW/K 2 m, ZT p = 0.024, PF n = 2.01 mW/K 2 m, and ZT p = 0.025). The findings presented in this work encourage further studies on the novel phases, TaFeAs in particular.

1. Introduction

The wide application of thermoelectric (TE) materials is desirable due to the accumulating heat losses of various origins (e.g., energy and automotive industries). Despite numerous comprehensive studies [1,2,3,4,5,6,7], the subject of TE performance among half-Heusler (hH) phases is still an interesting area of research. Many alloys and whole families of compounds (e.g., tellurides) were not investigated up to now [8,9,10,11]. The experimental reports for known systems are generally focused on one regime of carriers [12,13,14,15,16,17,18,19,20,21,22,23,24]. In many theoretical studies, the predictions of thermoelectric properies are based on the standard exchange-correlation functional (XCF) [25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53], which may be insufficient for an accurate description of electronic structures of semiconducting materials. Recent studies showed new directions for a search of novel hH-based systems with complex crystal structures [54,55].
In this study, the continuation of previous investigations of electronic structures in over 150 hH phases with eighteen valence electrons [11] is presented. The set of 34 novel and stable (according to the Born’s criteria and thermodynamics) potentially valuable TE materials, previously predicted, is a subject of careful examination of TE properties by ab initio calculations. The results were obtained with two XCF parametrizations. Electron- and hole-like carriers are considered. The relaxation time of the carriers and lattice thermal conductivity are approximated in physical units, which provides a limited number of reasonable candidate systems for further efforts in the experimental synthesis of novel hH phases.

2. Computational Details

Electronic structure calculations with spin-orbit coupling included, based on the density functional theory (DFT), were performed with the use of VASP [56,57,58,59]. The generalized gradient approximation (GGA [60]) and modified Becke–Johnson GGA (MBJGGA [61]) approaches were selected for XCF. The cut-off energy for a plane-wave basis was set to 500 eV. The meshes of k-points for electronic structure and transport calculations were set to 20 3 and 50 3 points in the Brillouin zone, respectively. The equilibrium lattice parameters and elastic constants were only obtained with the GGA calculations. The relaxation time of carriers was approximated with the deformation potential theory [62]. The lattice thermal conductivity calculations were performed following the Slack formula [63,64,65]. Transport analysis were carried out with the Boltzmann transport equation, employing the BoltzTraP2 package [66]. The diagram of the complete calculation scheme is shown in Figure 1. The open quantum materials database was employed for studies of thermodynamic stability [67].

3. Results and Discussion

A primary set of potentially valuable hH alloys was created based on the valence electron count (VEC) of eighteen (due to the expected stability of such systems [68]). The elastic stability and thermodynamic stability of these phases were carefully examined in the previous study [11]. It is worth recalling that the criteria of mechanical stability (related to elastic constants, i.e., born conditions) are easy to meet in cubic hH materials. The thermodynamic stability of the phases considered here was evaluated based on analyses of formation energies (convex hull and hull distance, Δ E H D ). The phases with zero hull distance, which means that their ground state crystal structures are expected to be cubic, were selected for further investigations of electronic structure. One may consider that other systems with formation energies close to the convex hull may also be stable or metastable in some nonequilibrium conditions. Only systems far from the most favorable composition should decompose into another phase or a mixture of various phases. Some additional hH compounds with very small hull distances (below 50 meV with respect to the lower crystal symmetry Pnma), e.g., TiNiPb, NbRuAs, and VRhGe, are not considered in this work. Although the results of ab initio calculations are usually reasonable, the synthesis of novel materials may be challenging even in the case of selected 34 phases. It is worth recalling that the thermal properties of similar materials may also be successfully studied with the use of molecular dynamics [69,70].
Only the semiconducting hH phases are selected here for further investigations. The band structures of the systems considered are depicted in Figure 2. One may find some general trends in electronic structure of hH materials, e.g., the band gaps ( E g ) are narrow for relatively heavy elements (when comparing TiCoAs and TiCoBi, VCoGe and VIrGe, NbCoGe, and NbIrGe). Indirect band gaps are found in most systems (with the MBJGGA approach in particular). Direct E g Γ Γ is predicted with GGA for ScPdBi, ZrRuTe, HfPtGe, and TaRuAs, whereas the MBJGGA approach yields indirect E g Γ X in these compounds. Only HfIrSb exhibits a direct E g within both XCFs. The use of more than one XCF and the effects of spin-orbit coupling on the electronic structures of hH systems were extensively discussed in previous reports (e.g., for YPdAs [11] and ScPtSb [8]). As presented in Table 1, the effect of MBJ potential on band gaps may be surprising, i.e., the direct-indirect transitions of E g result in narrower band gaps in some cases. Nevertheless, the general comparison between GGA- and MBJGGA-derived E g , depicted in Figure 3a, clearly indicates wider band gaps yielded by this approach.
The thermoelectric performance of materials is characterized by two quantities: the power factor (PF) and the figure of merit (ZT), which are defined as follows:
P F = S 2 σ ( τ ) Z T = P F κ e + κ L T .
where S is Seebeck coefficient, σ ( τ ) is electrical conductivity ( τ is relaxation time of carriers), κ e is thermal conductivity of carriers, κ L is lattice thermal conductivity, and T is temperature. The values of the Seebeck coefficient, σ , and κ e are calculated from first principles with the Boltzmann transport equation. The theoretical estimations of τ and κ L require careful discussion. According to the deformation potential theory [62], the relaxation time of carries may be computed with the formula given as follows:
τ β = 2 2 π C β ħ 3 ( k B T m ) 3 / 2 E β 2
where C β is the elastic constant in the β direction, m is the effective mass, and E β is the deformation potential in the β direction. The effective masses of electron- and hole-like carries, determined in the vicinity of the valence band maximum (VBM) and conduction band minimum (CBM), are gathered in Table 1. The relatively low values of m are desirable for TE performance due to the long relaxation time of carriers. The strong variations between the p- and n-type m occurred in some cases (ZrPtGe, TiCoAs), which promotes particular carrier regimes. The use of two XCF also may lead to similar discrepancies in m . The relaxation time of carriers is also strongly dominated by the deformation potential and elastic constant (not shown here), which generally results in good accordance between τ G G A and τ M B J for a particular carrier type. Investigations of the deformation potential E β M B J for VBM and CBM revealed remarkably high values for some compounds, i.e, VFeSb: 101.00 and 121.66 eV; NbFeAs: 113.89 and 102.74 eV; and TiFeTe: 101.86 and 119.86 eV, respectively. These results are clearly larger than those for similar systems (e.g., HfCoAs: 31.11; and 36.20 eV, for VBM and CBM). One may consider some limitations of the MBJ functional in the case of non-equillibrium conditions and prefer the GGA deformation potentials for VFeSb (36.17 and 36.70 eV), NbFeAs (39.31 and 39.20 eV), and TiFeTe (39.29 and 41.40 eV, for VBM and CBM, respectively). This issue is crucial for the further discussion of TE performance, which could be strongly overestimated due to the extraordinarily long relaxation time of carriers.
The total thermal conductivity of semiconducting materials is dominated by the lattice contribution ( κ L ), which may be approximated with the Slack’s formula [63]. This approach also employs elastic constants of materials. The DFT-derived data are presented in Table 2. The elastic constants C 12 are inversely proportional to the unit cell volumes of particular compounds. As depicted in Figure 3b, the obtained values of k L exhibit a pronounced dependence on the cubic lattice parameter. Compounds with relatively small lattice parameters and high values of C 12 are characterized by low lattice thermal conductivity. Furthermore, there is no clear interplay between k L and the mass of the particular constituent elements. For some systems, the big summary ionic mass of the phase leads to the high k L (e.g., HfCoAs and HfCoBi), and for others, the relation is the opposite (e.g., HfNiSn and HfPdSn). One may only consider some general trends, i.e., the presence of heavy X and Z elements leads to increased lattice thermal conductivity in X Y Z hH systems, whereas for the heavy Y ions an opposite relation is seen.
Figure 4 and Figure 5 depict distributions of the PF and ZT at 300 K. The values of the carrier concentrations considered are limited to the range from ≈10 19 to ≈10 21 , which seems to be reasonable according to the experimental results for hH phases [71]. The optimal values of carrier concentrations, maximalizing PF and ZT, are presented in the bottom panels of Figure 4 and Figure 5, respectively. As seen if Figure 4, the low values of PF (below ≈0.5 mW/K 2 m) are found for most systems for both p- and n-type regimes. The phases with the highest PF are TaFeAs, TaFeSb, VFeAs, and TiRuAs. As discussed above, the MBJGGA results for other materials may be biased by the artificially long relaxation time of carries. Namely, the high values of PF for ZrRuTe may be questionable because the corresponding GGA results are significantly smaller. The careful analysis of MBJGGA and GGA band structures for this material, seen in Figure 2, suggests only minor change in E g and similar shapes of VBM and CBM. The relatively high values of the PF in compounds gathered in Table 3 are obtained with both XCF approaches, which may confirm an optimal electronic structure for TE performance in these phases.
Although the values of the ZT presented in Figure 5 are proportional to the corresponding PF, the distribution of the ZT is slightly different from that of the PF. Namely, a relatively low lattice thermal conductivity of TiRhAs leads to the ZT comparable to those of VFeAs and TaFeAs. A similar effect is also found in MBJGGA results for VIrGe. The TE properties of compounds with the highest values of the PF and ZT, gathered in Table 3, are lower than those of the best hH TE systems (e.g., NbFeSb [36]). It is worth noting that the lattice thermal conductivity of hH materials may be poorly modeled with Slack’s formula, whereas the values of the PF predicted here suggest very good TE performance of selected phases.
Regarding available experimental data, Downie et al. disclosed the n-type PF of ≈0.3 mW/mK 2 in HfNiSn at 330 K with the possibility of further enhancement for various dopings [72]. The accordance between the experimental and n-type PF of 0.34 mW/mK 2 calculated within the MBJGGA approach for HfNiSn is very good. Although experimental studies reported extremely high ZT at temperatures above 800 K for various HfNiSn materials [73], the room temperature ZT in this system is negligible. While a good accordance between the experimental results and our predictions was obtained for TaFeSb [74], the systems based on this compound are also considered as promising only at high T. In the case of TaCoSn [75], the low experimental ZT at 300 K in the n-type regime is also predicted in the present work. However, the p-type TE performance of TaCoSn seems to be an interesting area for further research. For some compounds studied here, the available experimental reports were focused on properties other than TE performance, e.g., the electrical transport in ScPdBi [76]. Syntheses of most novel materials predicted in this work were not reported up to now. The generally good accordance between the theoretical and available experimental data for hH phases indicates the DFT-based methods combined with the transport modeling and approximations for the relaxation time of carriers and lattice thermal conductivity as a reasonable tool for materials science. Furthermore, an enhancement of TE in particular phases can be achieved with experimental methods, which are not considered in theoretical calculations, e.g., the modifications of a crystal microstructure cause a strong decrease in lattice thermal conductivity and higher ZT, as reported for VFeSb [77]. It is worth recalling that NbFeSb, which has been extensively studied in recent years, exhibits the PF and ZT of ≈33 μ W/cm K 2 and the ZT ≈ 0.4 at 300 K, respectively [16]. Furthermore, some of the properties of the two systems considered here were reported in previous theoretical studies, i.e., HfPdSn (half-Heusler compounds having full-Heusler counterparts regarded in machine learning models [21]) and ZrPdSn (electronic properties calculated with ab initio methods [78]).
Among the novel hH systems regarded in this work, TaFeAs is the most promising TE material. Other hH arsenides were recently predicted to show interesting properties [79,80,81,82]. Interestingly, the band gaps of materials selected here are not optimal according to the ‘10 k B T rule’ [83]. It is worth recalling that the band structures of promising TE materials among hH alloys, i.e., LaPtSb [84], exhibit characteristic shapes of VBM or CBM, which results in very low effective masses of carriers. The Fe-based arsenides exhibit such features in both p- and n-type regimes. The best n-type TE performance of TaFeAs may be explained to some extent by the significant flattening of conduction bands in the Γ X L direction without the flattening of CBM (high S and low m ). Furthermore, the Fe-bearing phases are also expected to be promising TE materials in the p-type regime, which is related to the relatively high carrier relaxation time and low lattice thermal conductivity in such systems.

4. Conclusions

Based on the high-throughput ab initio calculations, including the two XCFs regarded, the deformation potential theory applied for carrier relaxation time, and Slack’s formula for lattice thermal conductivity, the potential TE performance of 34 novel hH semiconductors was carefuly examined. TaFeAs, TaFeSb, VFeAs, and TiRuAs attract particular attention due to their relatively high PF and ZT at 300 K. The desired set of electronic properties for TE materials is difficult to obtain in a single hH phase. Such a complex examination of multiple parameters for each potentially valuable TE material reveals some interesting trends for predicting novel systems. The hH arsenides are expected to be a promising family of novel hH phases for further experimental research.

Author Contributions

K.B.: Methodology, investigation, visualization, and writing; M.J.W.: conceptualization, methodology, resources, writing, and supervision. All authors have read and agreed to the published version of the manuscript.

Funding

Calculations for this work were performed at the Wroclaw Center for Networking and Supercomputing (Project No. 158).

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationship that could have appeared to influence the work reported in this paper.

References

  1. Wang, G.; Wang, D. Electronic structure and thermoelectric properties of Pb-based half-Heusler compounds: ABPb (A = Hf, Zr; B = Ni, Pd). J. Alloys Compd. 2016, 682, 375–380. [Google Scholar] [CrossRef]
  2. Touia, A.; Benyahia, K.; Tekin, A. First-principles calculations of structural, electronic, optical, and thermoelectric properties of LuNiBi and LuNiSb half-heusler. J. Supercond. Nov. Magn. 2021, 34, 2689–2698. [Google Scholar] [CrossRef]
  3. Xing, G.; Sun, J.; Li, Y.; Fan, X.; Zheng, W.; Singh, D.J. Electronic fitness function for screening semiconductors as thermoelectric materials. Phys. Rev. Mater. 2017, 1, 065405. [Google Scholar] [CrossRef]
  4. Abid, O.M.; Menouer, S.; Yakoubi, A.; Khachai, H.; Omran, S.B.; Murtaza, G.; Prakash, D.; Kgebata, R.; Verma, K.D. Structural, electronic, elastic, thermoelectric and thermodynamic properties of the NbMSb half heusler (M = Fe, Ru, Os) compounds with first principle calculations. Superlattices Microstruct. 2016, 93, 171–185. [Google Scholar] [CrossRef]
  5. Touia, A.; Benkhaled, M.; Khobzaoui, C.; Fodil, M. Optical and Thermodynamic Properties of Half-Heusler Compound TaIrSn: Using Modified Becke–Johnson (mBJ). J. Supercond. Nov. Magn. 2021, 34, 2865–2877. [Google Scholar] [CrossRef]
  6. Candan, A.; Kushwaha, A.K. A first-principles study of the structural, electronic, optical, and vibrational properties for paramagnetic half-Heusler compound TiIrBi by GGA and GGA + mBJ functional. Mater. Today Commun. 2021, 27, 102246. [Google Scholar] [CrossRef]
  7. Chibani, S.; Arbouche, O.; Zemouli, M.; Benallou, Y.; Amara, K.; Chami, N.; Ameri, M.; El Keurti, M. First-principles investigation of structural, mechanical, electronic, and thermoelectric properties of Half-Heusler compounds RuVX (X = As, P, and Sb). Comput. Condens. Matter 2018, 16, e00312. [Google Scholar] [CrossRef]
  8. Winiarski, M.J.; Bilińska, K.; Ciesielski, K.; Kaczorowski, D. Thermoelectric performance of p-type half-Heusler alloys ScMSb (M = Ni, Pd, Pt) by ab initio calculations. J. Alloys Compd. 2018, 762, 901–905. [Google Scholar] [CrossRef]
  9. Winiarski, M.J.; Bilińska, K. High thermoelectric power factors of p-type half-Heusler alloys YNiSb, LuNiSb, YPdSb, and LuPdSb. Intermetallics 2019, 108, 55–60. [Google Scholar] [CrossRef]
  10. Winiarski, M.J.; Bilinska, K. Power Factors of p-type Half-Heusler Alloys ScNiBi, YNiBi, and LuNiBi by ab initio Calculations. Acta Phys. Pol. A 2020, 138, 533–538. [Google Scholar] [CrossRef]
  11. Bilińska, K.; Winiarski, M.J. Search for semiconducting materials among 18-electron half-Heusler alloys. Solid State Commun. 2023, 365, 115133. [Google Scholar] [CrossRef]
  12. Mallick, M.M.; Vitta, S. Thermoelectric properties of ultra-low thermal conductivity half-Heusler alloy. AIP Conf. Proc. 2016, 1731, 110027. [Google Scholar]
  13. Kimura, Y.; Zama, A. Thermoelectric properties of p-type half-Heusler compound HfPtSn and improvement for high-performance by Ir and Co additions. Appl. Phys. Lett. 2006, 89, 172110. [Google Scholar] [CrossRef]
  14. Sekimoto, T.; Kurosaki, K.; Muta, H.; Yamasaka, S. Thermoelectric and thermophysical properties of TiCoSb, ZrCoSb, HfCoSb prepared by SPS. In Proceedings of the ICT 2005. 24th International Conference on Thermoelectrics, Clemson, SC, USA, 19–23 June 2005; pp. 347–350. [Google Scholar]
  15. Serrano-Sánchez, F.; Luo, T.; Yu, J.; Xie, W.; Le, C.; Auffermann, G.; Weidenkaff, A.; Zhu, T.; Zhao, X.; Alonso, J.A.; et al. Thermoelectric properties of n-type half-Heusler NbCoSn with heavy-element Pt substitution. Phys. Chem. Chem. Phys. 2020, 8, 14822–14828. [Google Scholar] [CrossRef]
  16. Joshi, G.; He, R.; Engber, M.; Samsonidze, G.; Pantha, T.; Dahal, E.; Yang, J.; Lan, Y.; Kozinsky, B.; Ren, Z. NbFeSb-based p-type half-Heuslers for power generation applications. J. Am. Chem. Soc. 2019, 12, 4070–4076. [Google Scholar] [CrossRef]
  17. Zakutayev, A.; Zhang, X.; Nagaraja, A.; Yu, L.; Lany, S.; Mason, T.O.; Ginley, D.S.; Zunger, A. Theoretical prediction and experimental realization of new stable inorganic materials using the inverse design approach. J. Am. Chem. Soc. 2013, 135, 10048–10054. [Google Scholar] [CrossRef]
  18. Kawaharada, Y.; Kurosaki, K.; Muta, H.; Uno, M.; Yamanaka, S. High temperature thermoelectric properties of CoTiSb half-Heusler compounds. J. Alloys Compd. 2004, 384, 308–311. [Google Scholar] [CrossRef]
  19. Asaad, M.; Buckman, J.; Smith, R.I.; Bos, J.W.G. Thermoelectric properties and high-temperature stability of the Ti1-xVxCoSb1-xSnx half-Heusler alloys. RSC Adv. 2016, 6, 5611–56517. [Google Scholar] [CrossRef]
  20. Hooshm Zaferani, S.; Darebaghi, A.; Hong, S.J.; Vashaee, D.; Ghomashchi, R. Experimental realization of heavily p-doped half-Heusler CoVSn compound. Energies 2020, 13, 1459. [Google Scholar] [CrossRef]
  21. Gzyl, A.S.; Oliynyk, A.O.; Mar, A. Half-heusler structures with full-heusler counterparts: Machine-learning predictions and experimental validation. Cryst. Growth Des. 2020, 20, 6469–6477. [Google Scholar] [CrossRef]
  22. Romaka, V.V.; Romaka, L. Experimental and theoretical investigation of the Y–Ni–Sb and Tm–Ni–Sb systems. J. Alloys Compd. 2021, 855, 157334. [Google Scholar] [CrossRef]
  23. Oestreich, J.; Probst, U.; Richardt, F.; Bucher, E. Thermoelectrical properties of the compounds ScM {sup VIII} Sb and YM {sup VIII} Sb (M {sup VIII} = Ni, Pd, Pt). J. Phys. Condens. Matter 2003, 15, 635. [Google Scholar] [CrossRef]
  24. Dai, C.K.; Song, Q.F.; Xie, L.; Liu, R.H.; Bai, S.Q.; Chen, L.D. Improving thermoelectric properties of ZrPtSn-based half-Heusler compound by Sb doping. Rare Met. 2021, 40, 2838. [Google Scholar] [CrossRef]
  25. Adetunji, B.I.; Adebambo, P.O.; Bamgbose, M.K.; Musari, A.A.; Adebayo, G.A. Predicting the elastic, phonon and thermodynamic properties of cubic HfNiX (X = Ge and Sn) Half Heulser alloys: A DFT study. Eur. Phys. J. B 2019, 92, 231. [Google Scholar] [CrossRef]
  26. Bendahma, F.; Mana, M.; Terkhi, S.; Cherid, S.; Bestani, B.; Bentata, S. Investigation of high figure of merit in semiconductor XHfGe (X = Ni and Pd) half-Heusler alloys: Ab-initio study. Comput. Condens. Matter 2019, 21, e00407. [Google Scholar] [CrossRef]
  27. Vikram, J.K.; Enamullah Alam, A. Bismuth based half-Heusler alloys with giant thermoelectric figures of merit. J. Mater. Chem. A 2017, 13, 6131. [Google Scholar] [CrossRef]
  28. Kaur, K.; Kumar, R.; Rai, D.P. A promising thermoelectric response of HfRhSb half Heusler compound at high temperature: A first principle study. J. Alloys Compd. 2018, 763, 1018. [Google Scholar] [CrossRef]
  29. Kaur, K.; Rai, D.P.; Thapa, R.K.; Srivastava, S. Structural, electronic, mechanical, and thermoelectric properties of a novel half Heusler compound HfPtPb. J. Appl. Phys. 2017, 122, 045110. [Google Scholar] [CrossRef]
  30. Rani, B.; Wani, A.F.; Sharopov, U.B.; Patra, L.; Singh, J.; Ali, A.M.; El-Rehim, A.F.A.; Khandy, S.A.; Dhiman, S.; Kaur, K. Electronic Structure-, Phonon Spectrum-, and Effective Mass-Related Thermoelectric Properties of PdXSn (X = Zr, Hf) Half Heuslers. Molecules 2022, 27, 6567. [Google Scholar] [CrossRef]
  31. Bamgbose, M.K. First-principles study of electronic structure and thermoelectric properties of p-type XIrSb (X = Ti, Zr and Hf) half-Heusler compounds. Mater. Sci. Semicond. Process 2021, 129, 105792. [Google Scholar] [CrossRef]
  32. Popoola, A.I.; Odusote, Y.A. The properties of NbRhGe as high temperature thermoelectric material. J. Appl. Phys. 2019, 11, 51–56. [Google Scholar]
  33. Osafile, O.E.; Nenuwe, O.N. Lattice dynamics and thermodynamic responses of XNbSn half-heusler semiconductors: A first-principles approach. J. Nig. Soc. Phys. Sci. 2021, 3, 121–130. [Google Scholar] [CrossRef]
  34. Kaur, K.; Kumar, R. On the possibility of thermoelectricity in half Heusler XRuSb (X = V, Nb, Ta) materials: A first principles prospective. J. Phys. Chem. Solids 2017, 110, 108–115. [Google Scholar] [CrossRef]
  35. Fang, T.; Zheng, S.; Zhou, T.; Yan, L.; Zhang, P. Computational prediction of high thermoelectric performance in p-type half-Heusler compounds with low band effective mass. Phys. Chem. Chem. Phys. 2017, 19, 4411–4417. [Google Scholar] [CrossRef] [PubMed]
  36. Naydenov, G.A.; Hasnip, P.J.; Lazarov, V.K.; Probert, M.I.J. Huge power factor in p-type half-Heusler alloys NbFeSb and TaFeSb. J. Phys. Mater. 2019, 2, 035002. [Google Scholar] [CrossRef]
  37. Wang, L.L.; Miao, L.; Wang, Z.Y.; Wei, W.; Xiong, R.; Liu, H.J.; Shi, J.; Tang, X.F. Thermoelectric performance of half-Heusler compounds TiNiSn and TiCoSb. J. Appl. Phys. 2009, 105, 013709. [Google Scholar] [CrossRef]
  38. Kaur, K.; Kumar, R. High temperature thermoelectric performance of p-type TaRhSn half Heusler compound: A computational assessment. Ceram. Int. 2017, 43, 15160. [Google Scholar] [CrossRef]
  39. Kaur, K.; Kumar, R. Giant thermoelectric performance of novel TaIrSn Half Heusler compound. Phys. Lett. A 2017, 381, 3760. [Google Scholar] [CrossRef]
  40. Hong, D.; Zeng, W.; Xin, Z.; Liu, F.S.; Tang, B.; Liu, Q.J. First-principles calculations of structural, mechanical and electronic properties of TiNi-X (X = C, Si, Ge, Sn, Pb) alloys. Int. J. Mod. Phys. 2019, 33, 1950167. [Google Scholar] [CrossRef]
  41. Kaur, K. TiPdSn: A half Heusler compound with high thermoelectric performance. Europhys. Lett. 2017, 117, 47002. [Google Scholar] [CrossRef]
  42. Zheng, W.; Lu, Y.; Li, Y.; Wang, J.; Hou, Z.; Shao, X. Structural and thermoelectric properties of Zr-doped TiPdSn half-Heusler compound by first-principles calculations. Chem. Phys. Lett. 2020, 741, 137055. [Google Scholar] [CrossRef]
  43. Adebambo, P.O.; Agbaoye, R.O.; Bamgbose, M.K.; Ayedun, F.; Solola, G.T.; Adebayo, G.A. Assessing the structural, electronic, elastic and thermoelectric properties of PtTiSn and PdLaBi transition metal alloys from the first-principles prospective. Mater. Sci. Semicond. 2021, 129, 105796. [Google Scholar] [CrossRef]
  44. Ma, H.; Yang, C.L.; Wang, M.S.; Ma, X.G.; Yi, Y.G. Effect of M elements (M = Ti, Zr, and Hf) on thermoelectric performance of the half-Heusler compounds MCoBi. J. Phys. D Appl. Phys. 2019, 52, 255501. [Google Scholar] [CrossRef]
  45. Umukoro, J.O.; Omagbemi, O.G.; Osafile, O.E. Effect of Spin-Orbit Coupling (SOC) On the Electronic and Thermoelectric Properties of Ticobi Half Heusler Alloy. J. Phys. Chem. Res. 2022, 4, 4. [Google Scholar]
  46. Yang, J.; Li, H.; Wu, T.; Zhang, W.; Chen, L.; Yang, J. Evaluation of half-Heusler compounds as thermoelectric materials based on the calculated electrical transport properties. Adv. Funct. Mater. 2008, 18, 2880–2888. [Google Scholar] [CrossRef]
  47. Kaur, K.; Kumar, R. Ti based half Heusler compounds: A new on the screen with robustic thermoelectric performance. J. Alloys Compd. 2017, 727, 1171–1177. [Google Scholar] [CrossRef]
  48. Cherifi, F.; Mostefa, Z.; Boukra, A.; Meghoufel, Z.F.; Bouattou, M.; Kadi Allah, F.; Terki, F. Thermoelectric Transport Parameters of p-Type RuVAs and RuNbAs Heusler Alloys. Phys. Status Solidi B 2020, 257, 2000271. [Google Scholar] [CrossRef]
  49. Mostari, F.; Rahman, M.A.; Khatun, R. First principles study on the structural, elastic, electronic and optical properties of cubic ‘half-Heusler’alloy RuVAs under pressure. Int. J. Mat. Math. Sci. 2020, 2, 51–63. [Google Scholar]
  50. Sharma, S.K.; Ahmed, S.S. Transport properties of RuV-based half-Heusler semiconductors for thermoelectric applications: A computational study. J. Phys. Condens. Mat. 2020, 32, 405501. [Google Scholar]
  51. Sarwan, M.; Shukoor, A.; Singh, S. A first principle study of structural, elastic, electronic and thermodynamic properties of Half-Heusler compounds; YNiPn (Pn = As, Sb, and Bi). Solid State Sci. 2021, 112, 106507. [Google Scholar] [CrossRef]
  52. Sharma, S.; Kumar, P. Tuning the thermoelectric properties of YNiBi half-Heusler alloy. MRX 2018, 5, 046528. [Google Scholar] [CrossRef]
  53. Kong, F.; Hu, Y.; Hou, H.; Liu, Y.; Wang, B.; Wang, L. Thermoelectric and thermodynamic properties of half-Heulser alloy YPdSb from first principles calculations. J. Solid State Chem. 2012, 196, 511–517. [Google Scholar] [CrossRef]
  54. Nenuwe, N.O.; Agbawe, N.O. Ab initio predictions of thermoelectric, mechanical, and phonon characteristics of FeTiSe half-Heusler compound. Curr. Appl. Phys. 2023, 53, 132–141. [Google Scholar] [CrossRef]
  55. Wang, C.; Zhou, X.; Cong, D.; Tang, G.; Yang, J. A novel valence-balanced double half-Heusler Ti2Zr2Hf2NbVFe5Ni3Sb8 alloy by high entropy engineering. Mater. Today Phys. 2023, 36, 101172. [Google Scholar] [CrossRef]
  56. Kresse, G.; Hafner, J. Ab initio molecular dynamics for open-shell transition metals. Phys. Rev. B 1993, 48, 13115. [Google Scholar] [CrossRef]
  57. Kresse, G.; Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 1994, 49, 14251. [Google Scholar] [CrossRef]
  58. 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]
  59. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169. [Google Scholar] [CrossRef]
  60. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865. [Google Scholar] [CrossRef]
  61. Tran, F.; Blaha, P. Accurate band gaps of semiconductors and insulators with a semilocal exchange-correlation potential. Phys. Rev. Lett. 2009, 102, 226401. [Google Scholar] [CrossRef]
  62. Bardeen, J.; Shockley, W.J.P.R. Deformation potentials and mobilities in non-polar crystals. Phys. Rev. 1950, 80, 72. [Google Scholar] [CrossRef]
  63. Slack, G.A. Nonmetallic crystals with high thermal conductivity. J. Phys. Chem. Solids 1973, 34, 321–335. [Google Scholar] [CrossRef]
  64. Yang, K.; Wan, R.; Zhang, Z.; Lei, Y.; Tian, G. First-principle investigation on the thermoelectric and electronic properties of HfCoX (X = As, Sb, Bi) half-Heusler compounds. J. Solid State Chem. 2022, 312, 123386. [Google Scholar] [CrossRef]
  65. Sun, H.L.; Yang, C.L.; Wang, M.S.; Ma, X.G. Remarkably high thermoelectric efficiencies of the half-Heusler compounds BXGa (X = Be, Mg, and Ca). ACS Appl. Mater. Interfaces 2020, 12, 5838–5846. [Google Scholar] [CrossRef]
  66. Madsen, G.K.; Carrete, J.; Verstraete, M.J. BoltzTraP2, a program for interpolating band structures and calculating semi-classical transport coefficients. Comput. Phys. Commun. 2018, 231, 140–145. [Google Scholar] [CrossRef]
  67. Kirklin, S.; Saal, J.E.; Meredig, B.; Thompson, A.; Doak, J.W.; Aykol, M.; Rühl, S.; Wolverton, C. The Open Quantum Materials Database (OQMD): Assessing the accuracy of DFT formation energies. NPJ Comput. Mater. 2015, 1, 15010. [Google Scholar] [CrossRef]
  68. Gautier, R.; Zhang, X.; Hu, L.; Yu, L.; Lin, Y.; Sunde, T.O.; Chon, D.; Paeppelmeirer, K.; Zunger, A. Prediction and accelerated laboratory discovery of previously unknown 18-electron ABX compounds. Nat. Chem. 2015, 7, 308–316. [Google Scholar] [CrossRef]
  69. Lundgren, C.; Kakanakova-Georgieva, A.; Gueorguiev, G.K. A perspective on thermal stability and mechanical properties of 2D Indium Bismide from ab initio molecular dynamics. Nanotechnology 2022, 33, 335706. [Google Scholar] [CrossRef]
  70. Filho, M.A.M.; Hsiao, C.-L.; Batista dos Santos, R.; Hultman, L.; Birch, J.; Gueorguiev, G.K. Self-Induced Core–Shell InAlN Nanorods: Formation and Stability Unraveled by Ab Initio Simulations. ACS Nanosci. Au 2023, 3, 84–93. [Google Scholar] [CrossRef]
  71. Xia, K.; Liu, Y.; An, S.; Snyder, G.J.; Xin, J.; Yu, J.; Zhao, X.; Zhu, T. Enhanced thermoelectric performance in 18-electron Nb0.8CoSb half-heusler compound with intrinsic Nb vacancies. Adv. Funct. Mater. 2018, 28, 1705845. [Google Scholar] [CrossRef]
  72. Downie, R.A.; MacLaren, D.A.; Bos, J.W. Thermoelectric performance of multiphase XNiSn (X = Ti, Zr, Hf) half-Heusler alloys. J. Mater. Chem. A 2014, 2, 6107–6114. [Google Scholar] [CrossRef]
  73. Gürth, M.; Rogl, G.; Romaka, V.V.; Grytsiv, A.; Bauer, E.; Rogl, P. Thermoelectric high ZT half-Heusler alloys Ti1-x-yZrxHfyNiSn (0 < x < 1; 0 < x < 1). Acta Mater. 2016, 104, 210–222. [Google Scholar]
  74. Zhu, H.; Mao, J.; Li, Y.; Sun, J.; Wang, Y.; Zhu, Q.; Li, G.; Song, Q.; Zhou, J.; Fu, Y.; et al. Discovery of TaFeSb-based half-Heuslers with high thermoelectric performance. Nat. Commun. 2019, 10, 270. [Google Scholar] [CrossRef] [PubMed]
  75. Li, S.; Zhu, H.; Mao, J.; Feng, Z.; Li, X.; Chen, C.; Cao, F.; Liu, X.; Singh, D.J.; Ren, Z.; et al. n-Type TaCoSn-Based Half-Heuslers as Promising Thermoelectric Materials. ACS Appl. Mater. Interfaces 2019, 11, 41321–41329. [Google Scholar] [CrossRef]
  76. Ramarao, S.D.; Pawbake, A.; Singh, A.K.; Núñez-Regueiro, M.; Méasson, M.A.; Peter, S.C. Electrical transport properties of half-heusler ScPdBi single crystals under extreme conditions. J. Alloys Compd. 2020, 848, 156632. [Google Scholar] [CrossRef]
  77. Chauhan, N.S.; Miyazaki, Y. Low lattice thermal conductivity and microstructural evolution in VFeSb half-Heusler alloys. Materialia 2022, 22, 101430. [Google Scholar] [CrossRef]
  78. Çiftci, Y.Ö.; Çoban, C. Ab initio calculations on the structural, mechanical, electronic, dynamic, and optical properties of semiconductor half-heusler compound ZrPdSn. Z. Naturforsch. A 2016, 71, 135–143. [Google Scholar] [CrossRef]
  79. Gupta, Y.; Sinha, M.M.; Verma, S.S. First-principles investigation on the electronic, mechanical and lattice dynamical properties of novel AlNiX (X = As and Sb) half-Heusler alloys. Mater. Today Commun. 2021, 26, 101885. [Google Scholar] [CrossRef]
  80. Mokhtari, H.; Boumia, L.; Mokhtari, M.; Dahmane, F.; Mansour, D.; Khenata, R. Mechanical Stability, Electronic, and Magnetic Properties of XZrAs (X = Cr, Mn, V) Half-Heusler Compounds. J. Supercond. Nov. Magn. 2023, 36, 1217–1244. [Google Scholar] [CrossRef]
  81. Bamgbose, M.K. Electronic structure and thermoelectric properties of HfRhZ (Z = As, Sb and Bi) half-Heusler compounds. Appl. Phys. A 2020, 126, 1–8. [Google Scholar] [CrossRef]
  82. Razzaq, S.; Ismail, K.; Murtaza, G.; Raza, H.H. Theoretical Study of Half-Heusler CsXAs (X = Ca, Sr, and Ba) from First Principle Calculations. J. Supercond. Nov. Magn. 2022, 35, 3291–3299. [Google Scholar] [CrossRef]
  83. Mahan, G.D. Figure of merit for thermoelectrics. J. Appl. Phys. 1989, 65, 1578–1583. [Google Scholar] [CrossRef]
  84. Xue, Q.Y.; Liu, H.J.; Fan, D.D.; Cheng, L.; Zhao, B.Y.; Shi, J. LaPtSb: A half-Heusler compound with high thermoelectric performance. Phys. Chem. Chem. Phys. 2016, 18, 17912–17916. [Google Scholar] [CrossRef] [PubMed]
Figure 1. A scheme of the calculation process.
Figure 1. A scheme of the calculation process.
Crystals 13 01378 g001
Figure 2. Band structures of hH alloys obtained with GGA and MBJGGA, marked with black and red colors, respectively.
Figure 2. Band structures of hH alloys obtained with GGA and MBJGGA, marked with black and red colors, respectively.
Crystals 13 01378 g002
Figure 3. Band gaps and lattice thermal conductivity distributions. In (a), the relation between band gaps calculated with the GGA and MBJGGA parametrizations is depicted, taking into consideration Z-ion (XYZ) constituents. In (b), the lattice thermal conductivity as a function of the equilibrium lattice parameter is presented, with the elastic constant C 12 and volume of cell marked with the color and size of the markers, respectively.
Figure 3. Band gaps and lattice thermal conductivity distributions. In (a), the relation between band gaps calculated with the GGA and MBJGGA parametrizations is depicted, taking into consideration Z-ion (XYZ) constituents. In (b), the lattice thermal conductivity as a function of the equilibrium lattice parameter is presented, with the elastic constant C 12 and volume of cell marked with the color and size of the markers, respectively.
Crystals 13 01378 g003
Figure 4. A distribution of PF at 300 K for the novel hH systems regarded with XCF parametrizations denoted with color (black—GGA, red—MBJ) and carrier types with different markers’ shapes for p-type, i.e., holes (circle), and n-type, i.e., electrons (full dot). In the bottom part of the graph, the corresponding values of the carrier concentration for each combination (compound, XCF, and type of carrier) are presented.
Figure 4. A distribution of PF at 300 K for the novel hH systems regarded with XCF parametrizations denoted with color (black—GGA, red—MBJ) and carrier types with different markers’ shapes for p-type, i.e., holes (circle), and n-type, i.e., electrons (full dot). In the bottom part of the graph, the corresponding values of the carrier concentration for each combination (compound, XCF, and type of carrier) are presented.
Crystals 13 01378 g004
Figure 5. A distribution of ZT at 300 K for the hH systems regarded with XCF parametrizations denoted with color (black—GGA, red—MBJ) and carrier types with different markers’ shapes for p-type (circle), and n-type (full dot). In the bottom part of the graph, the corresponding values of the carrier concentrations for each combination (compound, XCF, type of carrier) are presented.
Figure 5. A distribution of ZT at 300 K for the hH systems regarded with XCF parametrizations denoted with color (black—GGA, red—MBJ) and carrier types with different markers’ shapes for p-type (circle), and n-type (full dot). In the bottom part of the graph, the corresponding values of the carrier concentrations for each combination (compound, XCF, type of carrier) are presented.
Crystals 13 01378 g005
Table 1. A comparison of band gaps E X C F (eV), effective mass m X C F p / n (m e ), and relaxation time values τ X C F p / n (fs) calculated for hH phases.
Table 1. A comparison of band gaps E X C F (eV), effective mass m X C F p / n (m e ), and relaxation time values τ X C F p / n (fs) calculated for hH phases.
CompoundE GGA E MBJ m GGA p m GGA n m MBJ p m MBJ n τ GGA p τ GGA n τ MBJ p τ MBJ n
ScPdBi0.07050.12110.344.800.190.3121.50.346.321.5
HfNiSn0.32390.27050.650.750.590.7818.513.118.510.8
HfPdSn0.37910.33210.680.60.770.6311.211.68.410.0
ZrPdSn0.43730.41680.520.820.510.9613.96.313.24.7
NbRuSb0.34750.4830.320.210.360.2231.664.86.119.2
TiNiGe0.63070.58121.371.201.231.055.65.95.96.5
VIrGe0.26850.64020.610.540.460.3814.018.120.632.6
VFeSb0.34360.65690.650.240.490.4114.261.52.82.5
NbFeAs0.57310.69870.321.430.330.3139.04.24.46.0
NbIrGe0.60310.72140.820.380.820.368.629.38.030.3
ZrPtPb0.66750.72940.480.190.470.3714.853.414.018.0
NbIrSn0.62510.72940.810.390.810.377.823.57.424.1
HfIrSb0.65820.75380.920.590.241.549.76.164.11.3
TiCoBi0.87860.7790.540.560.541.1314.211.813.33.8
TaRuAs0.36570.78180.360.230.440.2433.728.623.368.9
TiPtSn0.6650.80630.530.530.650.7713.011.69.16.0
TiRhBi0.65740.81380.390.580.551.1818.59.110.32.8
VFeAs0.36140.86250.390.330.390.3630.639.627.135.6
TaFeSb0.81350.87080.380.330.390.445.351.448.637.3
TiPtGe0.72320.90880.440.420.480.4719.619.016.015.2
HfCoBi0.97630.92210.540.590.580.4018.813.815.222.3
TiRhAs0.76840.95710.320.410.320.4231.620.229.819.4
TiIrSb0.68130.97160.350.280.360.2732.441.728.842.4
TaFeAs0.8790.98010.340.290.353.0152.565.945.61.9
VCoGe0.68320.99260.790.460.820.4911.022.59.320.1
HfPtGe0.92561.02360.560.440.251.3617.68.453.03.4
TaCoSn1.00831.03290.70.441.320.6515.727.45.614.0
ZrPtGe1.01191.08070.479.060.473.8319.80.218.30.6
NbCoGe1.08871.12720.580.340.600.3517.233.614.030.1
TiFeTe0.9841.1840.560.380.570.512.620.21.81.6
TaCoGe1.16011.18810.690.450.780.4620.032.214.226.7
TiCoAs1.30011.2410.701.700.831.5714.02.910.33.0
ZrRuTe0.9321.2470.510.340.240.2811.49.837.629.1
HfCoAs1.28611.36071.671.102.843.034.36.62.11.4
Table 2. Equilibrium lattice parameters (Å), elastic constants C 12 (GPa), unit cell volumes (Å 3 ), and values of lattice thermal conductivity (W/mK).
Table 2. Equilibrium lattice parameters (Å), elastic constants C 12 (GPa), unit cell volumes (Å 3 ), and values of lattice thermal conductivity (W/mK).
Compounda C 12 Volume κ L
ScPdBi6.52565.9369.4572.20
HfNiSn6.11173.7257.05113.52
HfPdSn6.36093.0964.3199.54
ZrPdSn6.39686.5565.4191.21
NbRuSb6.187125.2259.2176.27
TiNiGe5.66895.4045.5261.53
VIrGe5.818169.8849.2338.92
VFeSb5.78898.8348.4845.43
NbFeAs5.689115.0246.0363.29
NbIrGe6.010156.9454.2762.45
ZrPtPb6.50890.2768.9195.69
NbIrSn6.230132.1460.4576.60
HfIrSb6.333112.1963.50112.06
TiCoBi6.03375.0054.9058.22
TaRuAs5.972158.4953.2552.81
TiPtSn6.231108.2560.4865.36
TiRhBi6.28093.9561.9250.51
VFeAs5.496126.6241.5044.79
TaFeSb5.96098.4552.9397.50
TiPtGe5.991126.3553.7651.50
HfCoBi6.18865.6359.24102.54
TiRhAs5.889128.9651.0643.12
TiIrSb6.165119.5858.5873.10
TaFeAs5.692128.4346.1069.27
VCoGe5.512121.0741.8747.04
HfPtGe6.171124.5258.7572.45
TaCoSn5.96296.3652.9888.66
ZrPtGe6.200119.1159.5866.56
NbCoGe5.698116.3346.2569.23
TiFeTe5.86467.4350.4160.65
TaCoGe5.715125.2346.7479.43
TiCoAs5.605104.9244.0257.09
ZrRuTe6.29882.0262.4586.20
HfCoAs5.78399.4848.3576.60
Table 3. Compounds with the best predicted thermoelectric performance, based on PF and ZT at 300 K (see also Figure 4 and Figure 5). With indexes p / n , types of carriers are denoted. XCF used for presented results were GGA and MBJGGA (MBJ).
Table 3. Compounds with the best predicted thermoelectric performance, based on PF and ZT at 300 K (see also Figure 4 and Figure 5). With indexes p / n , types of carriers are denoted. XCF used for presented results were GGA and MBJGGA (MBJ).
Comp.PF GGA PF MBJ ZT GGA ZT MBJ
TaFeAs1.67 p , 2.01 n 1.48 p , 1.17 n 0.024 p , 0.025 n 0.021 p , 0.015 n
TaFeSb1.49 p , 1.15 n 1.63 p , 1.06 n 0.015 p , 0.011 n 0.017 p , 0.010 n
VFeAs1.10 p , 0.95 n 0.98 p , 1.25 n 0.024 p , 0.017 n 0.021 p , 0.025 n
TiRhAs0.95 p , 0.69 n 0.80 p , 0.96 n 0.018 p , 0.014 n 0.014 p , 0.021 n
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

Bilińska, K.; Winiarski, M.J. High-Throughput Exploration of Half-Heusler Phases for Thermoelectric Applications. Crystals 2023, 13, 1378. https://doi.org/10.3390/cryst13091378

AMA Style

Bilińska K, Winiarski MJ. High-Throughput Exploration of Half-Heusler Phases for Thermoelectric Applications. Crystals. 2023; 13(9):1378. https://doi.org/10.3390/cryst13091378

Chicago/Turabian Style

Bilińska, Kaja, and Maciej J. Winiarski. 2023. "High-Throughput Exploration of Half-Heusler Phases for Thermoelectric Applications" Crystals 13, no. 9: 1378. https://doi.org/10.3390/cryst13091378

APA Style

Bilińska, K., & Winiarski, M. J. (2023). High-Throughput Exploration of Half-Heusler Phases for Thermoelectric Applications. Crystals, 13(9), 1378. https://doi.org/10.3390/cryst13091378

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