
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39237544
55519
10.1038/s41598-024-55519-2
Article
Advancing thermoelectric potential: strontium telluride under compression strain analyzed with HSE hybrid functional and Wannier interpolation
Patel Hiren S.
Dabhi Vishnu A.
Vora Aditya M. voraam@gmail.com

https://ror.org/017f2w007 grid.411877.c 0000 0001 2152 424X Department of Physics, University School of Sciences, Gujarat University, Navrangpura, Ahmedabad, Gujarat 380 009 India
5 9 2024
5 9 2024
2024
14 207002 1 2024
24 2 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
In the present era, the energy sector is undergoing an intense transformation, which encourages numerous research efforts aimed at reducing and reusing energy waste. One of the main areas of focus is thermoelectric energy, where telluride compounds have attracted researchers due to their remarkable ability to convert thermal energy into electrical energy. We focused this study on finding out how well strontium telluride (SrTe) can be used to generate thermoelectric power by testing it under up to 10% compression strain. We have used advanced computational approaches to increase the accuracy of our results, specifically the HSE hybrid functional with the Wannier interpolation method. This method is primarily employed to analyze electronic properties; however, our research extends its utility to investigate thermoelectric characteristics. Our findings provide accurate predictions for both electronic and thermoelectric properties. The above method has successfully achieved a significant improvement of 58% in the electronic band gap value, resulting in a value of 2.83 eV, which closely matches the experimental results. Furthermore, the Figure of Merit 0.95 is obtained, which is close to the ideal range. Both the band gap value and the thermoelectric figure of merit decrease when the compression strain is increased. These findings emphasize the importance of using SrTe under specific conditions. The findings of this work provide motivation for future researchers to investigate the environmental changes in the thermoelectric potential of SrTe.

Keywords

Strontium telluride
Compression strain
Thermoelectric properties
Wannier interpolation
Hybrid functional
Subject terms

Materials science
Condensed-matter physics
Theory and computation
Physics
Condensed-matter physics
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

In today’s world, we are constantly working hard to meet increasing energy needs while also minimizing energy waste. These ongoing efforts, characterized by their approach, have opened up opportunities in the field of energy exploration. Thermoelectric energy has become an area in the past century within the diverse energy sector. Researchers are continuously pushing boundaries in this field by experimenting with materials. Among these materials, telluride compounds have gained attention, especially ‘SrTe’, for its ability to harness thermal energy and convert it into electricity1–3.

According to the existing body of literature, SrTe is naturally found in a NaCl-type crystal structure. However, it undergoes a phase transition, adopting a CsCl-type structure, when subjected to a pressure of 14.7 GPa4. Shi et al. calculated the structural stabilities, electronic, elastic, and optical properties of all possible phases of SrTe and found that only the NaCl phase undergoes a phase transition to the CsCl phase at 10.9 GPa5. High-pressure X-ray diffraction experiments by Zimmer et al. revealed a phase transition in SrTe from the NaCl phase to the CsCl phase at 12 GPa6. From all the above reviews, we can infer that SrTe has a phase transition between 10 and 15 GPa. The bulk modulus of SrTe was calculated to be 36 GPa using the FP-LAPW method by Khenata et al.7. By calculating the elastic and thermodynamic properties of SrTe at various temperatures, it was found that it is brittle at high temperatures8. These reviews indicate that SrTe has a low resistance to pressure, indicating that SrTe deforms at high pressure so that its mechanical properties may change. Rached et al. used the full-potential linear muffin-tin orbital method (FP-LMTO) along with a plane-wave basis (PLW) to find that SrTe has a band gap of 1.55 eV, a lattice constant of 6.64 Å, and a bulk modulus of 39.2 GPa9. SrTe has an indirect band gap of 2.76 eV at the Γ-X high-symmetry point, according to the GW approximation report by Zhu et al.10. Simons et al. confirmed the major participation of ‘p’ orbital’s electrons of ‘Te’ in the upper valence region, while ‘s’ orbital’s electrons are in the lower conduction region11. The lattice thermal conductivity of SrTe is a key parameter for its thermoelectric applications. It has been reported12 to be between 9 and 10 W m−1 K−1. It has been observed that the Seebeck coefficient exhibits higher values in p-type SrTe (100) and SrTe (110) as compared to the corresponding n-type doping concentrations13. The enhancement of thermoelectric properties can be achieved through the incorporation of ‘Mn’ and ‘Na’ doping into SrTe14. Tan et al. reported that hole-doping of PbTe–8%SrTe greatly boosted its power factor, with maximal values over 30 μW cm−1 K−2, and achieved a thermoelectric figure of merit (ZT) of 2.5 at 923 K15. Biswas et al. experimentally revealed that the highest thermoelectric figure of merit (ZT) observed for a p-type PbTe-SrTe system was ~ 1.3 at 715 K, and the p-type PbTe-SrTe system has the potential to impact power generation applications in the temperature range of 500–800 K16. Hexagonal-SrTe and hexagonal-SrS monolayers exhibit excellent thermoelectric properties, with a ZT > 0.8 for temperatures above 600 K17. Most of the above research has been obtained by traditional methods, and we have tried them with advanced computational methods to improve the accuracy of the final results.

From deep-sea exploration to space missions, thermoelectric material compounds are increasingly being used, and the evaluation of their performance under pressure conditions is significant. Specifically, this study examines the power generation capabilities of strontium telluride (SrTe) under applied compression strains ranging from 0 to 10%. In order to improve the accuracy of the results, the HSE hybrid functional and Wannier interpolation have been applied. It has been our goal to produce accurate predictions for both electronic and thermoelectric properties based on our findings. In this study, we observed a significant change in thermoelectric power generation under compression strain. Researchers in this field may find the results and method of this study to be valuable reference points in the future.

In the present study, we used SrTe with a rock salt structure and Fm3¯m space group. As shown in Fig. 1, Sr2+ is bonded to six equivalent Te2− atoms to form a mixture of corner and edge-sharing SrTe6 octahedra. Here, all the Sr–Te bond lengths are 3.36 Å. Also, Te2− is bonded to six equivalent Sr2+ atoms to form a mixture of corner and edge-sharing TeSr6 octahedra as well. Where, ‘Sr’ and ‘Te’ atoms have been located ‘4a’ (0, 0, 0) and ‘4b’ (0, 0, 0.5) as per Wyckoff positions respectively18–22.Figure 1 Atomic arrangements of SrTe, (a) conventional cell of SrTe with rock-salt formation and (b) primitive cell of SrTe having space group Fm3¯m.

Computational details

In this study, we employed Density Functional Theory (DFT) as implemented in the Quantum ESPRESSO (QE) software package for the electronic structure calculations23,24. The choice of DFT allows for accurate and efficient modeling of the electronic properties of materials to some extent25. QE, being a widely recognized and extensively validated computational tool, provided a robust framework for our investigations. The electronic wave functions were expanded using a plane wave basis set, ensuring a comprehensive representation of the electronic states. We used the Optimized Norm-Conserving Vanderbilt Pseudopotential (ONCV-PP) to account for the scalar relativistic effects26. This method is known for being accurate at describing relativistic effects while keeping computational efficiency high. We used the Perdew–Burke–Ernzerhof (PBE) semi-local exchange–correlation (XC) functional to completely optimize the cell volume so that our calculations would be accurate and converge27. The optimization process involved adjusting the cell volume until the system reached its minimum energy configuration, ensuring the accuracy of our structural predictions. For the numerical accuracy of our calculations, a kinetic energy cutoff of 120 Ry was employed. The optimization of the cell parameters was carried out using the BFGS algorithm, a quasi-Newton optimization algorithm known for its efficiency in minimizing energy with respect to the structural parameters. In addition, stringent convergence criteria were applied, with an energy convergence threshold set at 10−5 Ry/atom and a force convergence threshold set at 10−4 Ry/Bohr. These criteria ensured that both the total energy and atomic forces reached the level of precision necessary for reliable results28.

For refining the plane wave functions and determining the band gap value, we employed the HSE hybrid functional and Wannier interpolation29,30. This additional step allowed us to improve the accuracy of our electronic structure calculations and obtain more reliable predictions for the electronic band structure and band gap. For the present computation, we configured a three-dimensional mesh with a 2 × 2 × 2 q-sampling for the Fock operator. The electron energy convergence threshold was established at 10−8 Ry/electron, and the K-mesh grid was refined up to 0.15 Å31. Specifically, we selectively projected the ‘s’ orbital of ‘Sr’ and the ‘p’ orbital of ‘Te’ due to their significant involvement near the Fermi region. Thermodynamic properties were computed separately utilizing the ‘Gibbs 2’ code32,33. Furthermore, thermoelectric properties were assessed using both the wave functions obtained from ‘Quantum ESPRESSO (QE)’ and ‘Wannier interpolations’34. We utilized BoltzTraP-2 for comprehensive thermoelectric calculations, enhancing the depth of our investigation and ensuring accurate assessments of the material’s electronic and transport characteristics35. We have repeated the above steps under various compression strain conditions.

This combined methodology, characterized by its use of advanced computational tools, convergence criteria, and hybrid functional approaches, ensured the accuracy and reliability of our results in investigating the electronic and thermoelectric properties of the materials under consideration.

Results and discussion

The lattice parameter of the studied material was determined through an E-V fitting curve and energy minimization approach (Fig. 3a), resulting in a value of 6.7189 Å (Table 1). This finding demonstrates excellent agreement with both previously available theoretical predictions and experimental data, affirming the reliability and accuracy of our computational methodology.Table 1 Lattice parameter calculations for SrTe in angstroms, in comparison with experimental and theoretical references.

	Present work	Other work	Experimental	
Lattice parameter (Å)	6.7189	6.7336, 6.7217,37,38, 6.767	6.666,36	

Electronic properties

We have observed 1.78 eV indirect energy band gap between Γ and X high symmetry points in the first Brillouin zone for SrTe using electronic band structure calculation (Fig. 2). Interestingly, this band gap exhibited a decreasing trend with applied compression strain (Table 2). The observed behavior suggests that SrTe behaves as a semiconductor with a moderate band gap, and the increase in conductivity under applied compression strain highlights its potential for enhanced electronic transport properties.Figure 2 Electronic band structure of SrTe at different compression strains: (a) 0%, (b) 5%, and (c) 10%, depicted using selected high symmetry points. Blue band lines were obtained using standard Quantum ESPRESSO (QE) calculation, while red bands were derived using the HSE hybrid functional and Wannier interpolation, providing enhanced accuracy. (EF indicates Fermi level).

Table 2 Comparison of the change in energy band gap (eV) with applied compression strain (%), in relation to existing theoretical and experimental data.

Compression strain/band gap (eV)	Present work	Other work	Experimental	
Without HSE + Wannier	With HSE + Wannier	
0%	1.78	2.83	1.8036, 1.637, 1.737, 1.559	2.939	
5%	1.23	2.34	–	–	
10%	0.70	1.86	–	–	

It is known that, standard DFT functionals, like the Generalized Gradient Approximation (GGA) or the Local Density Approximation (LDA), may not adequately capture the electronic correlations and the self-interaction errors inherent in some systems40,41. These shortcomings can lead to inaccurate predictions of electronic properties, such as band gaps and electronic structures. Hybrid functional incorporates a fraction of exact exchange from Hartree–Fock theory, which is known to provide a more accurate description of electronic correlations42. By combining the exchange term from Hartree–Fock with the correlation part from standard DFT, hybrid functional aim to strike a balance between accuracy and computational efficiency. The HSE functional is specifically formulated to include a portion of screened Hartree–Fock exchange in addition to the standard DFT exchange–correlation functional. The inclusion of this non-local exchange term helps to improve the description of charge transfer, band gaps, and electronic structures in systems with strong correlation effects. HSE has been shown to provide more accurate results for a wide range of materials, including semiconductors, insulators, and transition metal complexes30. While, Wannier interpolation in DFT calculations allows for the transformation of electronic states from a plane-wave basis (used in typical DFT calculations) to localized Wannier functions. This transformation aids in constructing a more interpretable and computationally efficient description of the electronic structure, enabling the study of various material properties and behaviors34. These two methods collectively improved the accuracy of the electronic calculation at the cost of computational time. As far as p-DOS is concerned, the ONCV-PP demonstrates notable accuracy in predicting electronic and optical properties; it exhibits a limitation in its inability to calculate the projected Density of States (p-DOS), as its computational focus is primarily directed towards the valence part of the electronic structure43,44. As indicated in the literature review, we have conducted the aforementioned calculations exclusively for the ‘5p’ and ‘5s’ orbitals of the elements ‘Te’ and ‘Sr’ (Fig. 2). We achieved a band gap closer to experimental values (Table 2), affirming the reliability of our novel computational methodology for subsequent calculations. This advancement positions us to obtain more precise estimations of thermoelectric potential.

Thermodynamic properties

The change in energy versus volume of SrTe is shown in Fig. 3 under ± 10% applied strain, and its bulk modulus obtained using the Birch–Murnaghan fitting curve is found to be 32.51 GPa. For which the following equations have been used:1 PV=32B0V0V7/3-V0V5/31+34B0′-4V0V2/3-1

2 EV=E0+916V0B0V0V2/3-13B0′+V0V2/3-126-4V0V2/3

3 B0=-V∂P∂VP=0andB0′=∂B∂PP=0

where, P, B0, V0, V, B0′, E, and E0 are pressure, bulk modulus, reference volume, deformed volume, internal energy of the structure, and reference energy of the structure respectively.Figure 3 (a) Energy–volume (E–V) plot of SrTe with a second-order Birch–Murnaghan fitting curve, and (b) graphical representation depicting the variation in volume of SrTe in relation to temperature and compression strain.

As shown in Fig. 3b, the volume of SrTe increases with temperature and decreases under compression strain, which follows the rules of the general principle of physics. After applying this data in the Gibbs2 tool, we found the relation of bulk modulus (B) and Debye temperature (θD) with temperature for 0%, 5%, and 10% strains, and we found that both are increasing with the strain while decreasing with increasing temperature (Fig. 4a,b)32. This is a general trend, as intermolecular distances increase with increasing temperature, thereby weakening attractive forces between atoms such as London dispersion forces or hydrogen bonds, which are responsible for resisting compression in the material, resulting in a decrease in bulk modulus. When pressure is applied to a solid, it affects the inter-atomic distances, bond strengths, and overall crystal structure, which influences the bulk modulus as well as the vibrational modes of atoms in the crystal lattice leads to increasing trends in Debye temperature. While the Gruneisen parameter (ϒ) shows a reverse trend (Table 3)43. 4 Γ=B9RTθD3∫0θD/Tx4exex-12dx∂V∂TP

Figure 4 Temperature-dependent variations in (a) bulk modulus, (b) Debye temperature (θD), (c) Grüneisen parameter, and (d) lattice thermal conductivity of SrTe subjected to 0%, 5%, and 10% compression strains, providing comprehensive insights into its thermodynamic properties.

Table 3 Bulk modulus, Debye temperature, Gruneisen parameter, and lattice thermal conductivity values at 300 K under compression strains of 0%, 5%, and 10%: a comparative analysis with theoretical and experimental studies.

Compression strain	Present work	Other work	Experimental	
0%	5%	10%	0%	0%	
Bulk modulus (GPa)	30.69	62.06	98.01	31.8036, 31.9747, 3237,38, 367	39.536	
Debye Temperature (K)	192.34	258.65	326.36	17237, 20912	–	
Gruneisen parameter	1.93	1.69	1.53	1.4812	–	
Lattice thermal conductivity (W m−1 K−1)	8.10	19.01	33.59	10.6412, 10.548	9.6612	

At lower temperatures in semiconductor materials, electrons are typically in a bound state, emphasizing the significance of lattice thermal conductivity (KL) over electron thermal conductivity (Ke). As temperature rises, lattice vibrations become more dispersed, leading to a reduction in the mean free path and a subsequent decrease in lattice thermal conductivity. Conversely, at higher temperatures, the emergence of electron–hole pairs contributes to an increase in electron thermal conductivity. The compression strain-induced increase in phonon group velocity results in an elevated lattice thermal conductivity, as illustrated in the accompanying graph. Here the Slack’s equation is employed to calculate KL in this context12,45,46.5 KL=AM¯θD3δγ2Tn2/3

6 A=2.43×10-81-0.514γ+0.228γ2

In the above equations, we denote M¯ as the average atomic mass, θD as the Debye temperature, δ3 as the volume of the unit cell per atom, T as the temperature, n as the number of atoms per unit cell, and ϒ as the Gruneisen parameter.

Thermoelectric properties

The thermoelectric properties of SrTe were assessed using the semi-classical Boltzmann transport theory and the rigid band approach, which were implemented in the BoltzTraP algorithm44,49. The theoretical maximum efficiency of a thermoelectric material is determined by the hot end temperature (Th), the cold end temperature (Tc), and the material’s ZT value50. This relationship is stated by the following equation:7 ηmax=Th-TcTh×1+ZT-11+ZT+TcTh

Here, we denote T as the average temperature, and ZT is the figure of merit of the material, a dimensionless quantity51. Here, ZT depends on Seebeck’s coefficient (S), electrical conductivity (σ), and total thermal conductivity (K=κe+KL),8 ZT=S2σKT

The figure of merit for a high-quality thermoelectric material should ideally be equal to or greater than 1. The heat conductivity must be sufficiently low, but the electrical conductivity must be sufficiently high to meet this requirement. The value of the figure of merit (ZT) is influenced by temperature and charge carrier concentration. Therefore, it is important to choose the temperature and carrier concentration values carefully in order to achieve the optimum ZT value52. It is crucial to keep in mind that in the rigid band shift model, the chemical potential controls the compound’s carrier concentration. During n-type doping, the Fermi level increases, indicating a positive EF. Conversely, p-type doping causes the Fermi level to decrease, indicating a negative EF53. The Seebeck coefficient and both electrical and thermal conductivity exhibit an inverse relationship with carrier concentration. Therefore, it is not a straightforward assumption that increasing the carrier concentration will necessarily lead to a rise in the value of ZT. The establishment of an ideal value of ZT for a certain material is always true under ideal conditions. This is because when the carrier concentration increases, both the electrical conductivity and electron thermal conductivity increase, but the Seebeck coefficient decreases51.

This study was mostly about looking at SrTe at temperatures between 300 and 1200 K, with charge carrier concentrations of ± 1018 cm−3, ± 1020 cm−3, and ± 1022 cm−3, and compression strains of 0%, 5%, and 10%. Figs. 5b and 6b,e show the unexpected behavior of the semiconducting material ‘SrTe’ at high charge carrier concentrations in which the electrical conductivity remains unchanged with temperature. As we know, at high carrier concentrations, most of the available states in the conduction band are filled even at low temperatures. This saturation limits the increase in conductivity with temperature, as there are few empty states for excited carriers to occupy. This leads to an electrical conductivity that is relatively independent of the temperature. The investigation of figures also reveals that negative charge carriers exhibit greater electrical conductivity than positive charge carriers. This phenomenon stems from the dominance of electrons, the majority charge carriers in n-type semiconductors, which possess higher intrinsic mobility and lower effective mass than holes, the majority carriers in p-type semiconductors. Consequently, electrons experience less resistance when moving through the material under an applied electric field, thereby enhancing current flow and conductivity46. Figures 5c and 6c,f illustrate that the electron thermal conductivity of semiconductors exhibits a positive correlation with temperature due to the heightened density and mobility of charge carriers at elevated temperatures. It is noted that higher concentrations of negative charge carriers are correlated with elevated levels of thermal conductivity and electrical conductivity.Figure 5 Graph illustrating variation in (a) Seebeck co-efficient, (b) electrical conductivity, (c) thermal conductivity, and (d) figure of merit (ZT) with respect to temperature and doping concentration at 0% compression strain. (Line graph indicating value obtained from normal QE calculation, while square in graph indicating values derived from advanced computing using hybrid and Wannier interpolation).

Figure 6 Graphs illustrating the variation in (a) Seebeck coefficient, (b) electrical conductivity, and (c) Thermal conductivity with respect to temperature and doping concentration under 5% compression strain. Additional graphs depicting the variation in (d) Seebeck coefficient, (e) electrical conductivity, and (f) thermal conductivity with respect to temperature and doping concentration under 10% compression strain. (Line graphs represent values obtained from normal QE calculation, while squares on the graphs indicate values derived from advanced computing using hybrid and Wannier interpolation).

In semiconductors, the Seebeck coefficient exhibits a significant increase with temperature, owing to the growing influence of thermally excited carriers. The Seebeck coefficient in heavily doped semiconductors can display temperature-independent properties as a result of the displacement of the Fermi level. This is prominently evident in the depicted Figs. 5a and 6a,d. Furthermore, the ZT plot reveals that SrTe exhibits superior thermoelectric performance with a negative charge carrier concentration compared to a positive one (Figs. 5d and 7a,b).Figure 7 Figures depicting the variation in (a) ZT under 5% compression strain and (b) ZT under 10% compression strain concerning temperature. Line graphs represent values obtained from normal QE calculation, while squares on the graphs indicate values derived from advanced computing using hybrid and Wannier interpolation.

If you look closely at Fig. 5, you can see that there is no difference between the thermoelectric properties that come from the normal QE calculations and those that come from the advanced (HSE + Wannier) calculations when there is no compression strain. However, for 5% and 10% compression strain, the results diverge (Fig. 6). As expected, at higher charge carrier concentrations, the differences become insignificant. However, at high temperatures with a ± 1018 cm−3 charge carrier concentration, the normal thermoelectric calculations yield anomalous results. In contrast, the thermoelectric properties obtained with the advanced (HSE + Wannier) calculations remain stable across all temperature ranges and charge carrier concentrations. Our data also revealed a consistent trend: a charge carrier concentration of − 1020 cm−3 consistently yielded the highest ZT values across all compression strain ratios. This finding is further solidified by the ZT values listed in Table 4 at 1200 K for each strain ratio.Table 4 Figure of Merit (ZT) at 1200 K and − 1020 cm−3 charge carrier concentration under varying compression strains.

Compression strain	0%	5%	10%	
Figure of merit (ZT)	0.95	0.80	0.70	

Conclusion

This study highlights the synergistic power of combining Wannier functions and HSE hybrid functional within the Quantum Espresso (QE) framework for accurate and efficient characterization of electronic and thermoelectric properties in materials. Wannier functions offer a computationally efficient route for analyzing properties on a reduced k-point grid, while HSE hybrid functional significantly improves the accuracy of energy band gaps and wave functions. Our results demonstrate that this combined approach provides a favorable environment for calculating and interpreting reliable electronic and thermoelectric properties, especially under applied strain. Notably, this is the first attempt at this advanced technique for computing stable thermoelectric properties under strain, and the obtained results show better agreement than conventional QE calculations. We achieved a lattice parameter of 6.72 Å, with a minimal deviation (0.1%) from the experimental value. Our calculations revealed a decrease in the electronic band gap with increasing compression, and by employing the HSE hybrid functional and Wannier interpolation, we successfully reduced the band gap error to 2%, with a value of 2.83 eV, closer to the experimental data. The lattice thermal conductivity of SrTe was found to be 8.10 W m−1 K−1, and interestingly, the thermoelectric power generation ability was shown to decrease under compression strain. Notably, the highest figure of Merit (ZT) of 0.95 was attained at 1200 K with a negative charge carrier concentration of 1020 cm−3. However, here the decrease in thermoelectric figure of merit under the effect of compression strain raises questions about the use of telluride compounds as thermoelectric materials for applications under high pressure conditions such as submarines and space shuttles. Overall these findings showcase the potential of this combined approach for advanced material characterization and pave the way for further exploration of complex material properties with improved accuracy and computational efficiency.

Acknowledgements

The authors utilized the computer facility established by the Department of Science and Technology, Government of India, New Delhi, India, as a component of the DST-FIST initiative. We would like to extend our appreciation to the University Grant Commission of India, situated in New Delhi, for the financial assistance it offered under DRS-SAP-II.

Author contributions

H.S.P., V.A.D., A.M.V. studied the conception and designed the final approval of the version to be published. H.S.P. was associated with the calculations. H.S.P., V.A.D., A.M.V. contributed to the improvement of the scientific content of the article. All authors reviewed the manuscript.

Data availability

No datasets were generated or analysed during the current study.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Shi X-L Chen Z-G Quasi-one-dimensional bulk thermoelectrics Joule 2023 7 6 1108 1110 10.1016/j.joule.2023.05.008
Shi, X.-L. & Chen, Z.-G. Quasi-one-dimensional bulk thermoelectrics. Joule 7(6), 1108–1110. 10.1016/j.joule.2023.05.008 (2023).10.1016/j.joule.2023.05.008
2. Shi X A solvothermal synthetic environmental design for high-performance SnSe-based thermoelectric materials Adv. Energy Mater. 2022 12 20 2200670 10.1002/aenm.202200670
Shi, X. et al. A solvothermal synthetic environmental design for high-performance SnSe-based thermoelectric materials. Adv. Energy Mater. 12(20), 2200670. 10.1002/aenm.202200670 (2022).10.1002/aenm.202200670
3. Shi X Advances in flexible inorganic thermoelectrics EcoEnergy 2023 1 2 296 343 10.1002/ece2.17
Shi, X. et al. Advances in flexible inorganic thermoelectrics. EcoEnergy 1(2), 296–343. 10.1002/ece2.17 (2023).10.1002/ece2.17
4. Kunduru L Roshan R Yedukondalu N Sainath M Pressure driven topological semi metallic phase in SrTe Nucleation Atmos. Aerosols 1966 1 2018 10.1063/1.5038708
Kunduru, L., Roshan, R., Yedukondalu, N. & Sainath, M. Pressure driven topological semi metallic phase in SrTe. Nucleation Atmos. Aerosols 1, 2018. 10.1063/1.5038708 (1966).10.1063/1.5038708
5. Shi L Duan Y Yang X Qin L Structural stabilities, electronic, elastic and optical properties of SrTe under pressure: A first-principles study Phys. B Condens. Matter 2011 406 2 181 186 10.1016/j.physb.2010.10.038
Shi, L., Duan, Y., Yang, X. & Qin, L. Structural stabilities, electronic, elastic and optical properties of SrTe under pressure: A first-principles study. Phys. B Condens. Matter 406(2), 181–186. 10.1016/j.physb.2010.10.038 (2011).10.1016/j.physb.2010.10.038
6. Zimmer H-G Winzen H Syassen K High-pressure phase transitions in CaTe and SrTe Phys. Rev. 1985 32 6 4066 4070 10.1103/physrevb.32.4066
Zimmer, H.-G., Winzen, H. & Syassen, K. High-pressure phase transitions in CaTe and SrTe. Phys. Rev. 32(6), 4066–4070. 10.1103/physrevb.32.4066 (1985).10.1103/physrevb.32.4066
7. Khenata R First-principle study of structural, electronic and elastic properties of SrS, SrSe and SrTe under pressure Phys. B Condens. Matter 2003 339 4 208 215 10.1016/j.physb.2003.07.003
Khenata, R. et al. First-principle study of structural, electronic and elastic properties of SrS, SrSe and SrTe under pressure. Phys. B Condens. Matter 339(4), 208–215. 10.1016/j.physb.2003.07.003 (2003).10.1016/j.physb.2003.07.003
8. Saoud FS Rabah K Bouhemadou A Plenet JC Henini M Djabou REH Structural stabilities and elastic thermodynamic properties of SrTe compound and SrTe1−xCax alloy under high pressure J. Electron. Mater. 2016 46 2 766 774 10.1007/s11664-016-4830-9
Saoud, F. S. et al. Structural stabilities and elastic thermodynamic properties of SrTe compound and SrTe1 − xCax alloy under high pressure. J. Electron. Mater. 46(2), 766–774. 10.1007/s11664-016-4830-9 (2016).10.1007/s11664-016-4830-9
9. Rached D Rabah M Benkhettou N Soudini B Abid H Electronic and structural properties of strontium chalcogenides SrS, SrSe and SrTe Phys. Status Solidi 2004 241 11 2529 2537 10.1002/pssb.200402053
Rached, D., Rabah, M., Benkhettou, N., Soudini, B. & Abid, H. Electronic and structural properties of strontium chalcogenides SrS, SrSe and SrTe. Phys. Status Solidi 241(11), 2529–2537. 10.1002/pssb.200402053 (2004).10.1002/pssb.200402053
10. Zhu XL Zhang H Cheng X Theoretical study of the band gap of SrSe and SrTe in GW approximation Mod. Phys. Lett. B 2011 25 23 1905 1914 10.1142/s0217984911027121
Zhu, X. L., Zhang, H. & Cheng, X. Theoretical study of the band gap of SrSe and SrTe in GW approximation. Mod. Phys. Lett. B 25(23), 1905–1914. 10.1142/s0217984911027121 (2011).10.1142/s0217984911027121
11. Simons J Hempelmann J Fries KS Müller PC Dronskowski R Steinberg S Bonding diversity in rock salt-type tellurides: Examining the interdependence between chemical bonding and materials properties RSC Adv. 2021 11 34 20679 20686 10.1039/d1ra02999a 35479374
Simons, J. et al. Bonding diversity in rock salt-type tellurides: Examining the interdependence between chemical bonding and materials properties. RSC Adv. 11(34), 20679–20686. 10.1039/d1ra02999a (2021).35479374 10.1039/d1ra02999a
12. Roshan R Anomalous lattice thermal conductivity in rocksalt IIA–VIA compounds ACS Appl. Energy Mater. 2021 5 1 882 896 10.1021/acsaem.1c03310
Roshan, R. et al. Anomalous lattice thermal conductivity in rocksalt IIA–VIA compounds. ACS Appl. Energy Mater. 5(1), 882–896. 10.1021/acsaem.1c03310 (2021).10.1021/acsaem.1c03310
13. Li J Wang Y Zhang G Chen D Ren F First-principles investigation of the electronic structures and Seebeck coefficients of PbTe/SrTe interfaces J. Appl. Phys. 2019 125 3 035107 10.1063/1.5053710
Li, J., Wang, Y., Zhang, G., Chen, D. & Ren, F. First-principles investigation of the electronic structures and Seebeck coefficients of PbTe/SrTe interfaces. J. Appl. Phys. 125(3), 035107. 10.1063/1.5053710 (2019).10.1063/1.5053710
14. Luo J Enhanced average thermoelectric figure of merit of the PbTe–SrTe–MnTe alloy ACS Appl. Mater. Interfaces 2017 9 10 8729 8736 10.1021/acsami.6b16060 28256136
Luo, J. et al. Enhanced average thermoelectric figure of merit of the PbTe–SrTe–MnTe alloy. ACS Appl. Mater. Interfaces 9(10), 8729–8736. 10.1021/acsami.6b16060 (2017).28256136 10.1021/acsami.6b16060
15. Tan G Non-equilibrium processing leads to record high thermoelectric figure of merit in PbTe–SrTe Nat. Commun. 2016 7 1 12167 10.1038/ncomms12167 27456303
Tan, G. et al. Non-equilibrium processing leads to record high thermoelectric figure of merit in PbTe–SrTe. Nat. Commun. 7(1), 12167. 10.1038/ncomms12167 (2016).27456303 10.1038/ncomms12167
16. Biswas K Investigation of the thermoelectric properties of the PbTe-SrTe system MRS Proc. 2010 1267 187 192 10.1557/proc-1267-dd06-05
Biswas, K. et al. Investigation of the thermoelectric properties of the PbTe-SrTe system. MRS Proc. 1267, 187–192. 10.1557/proc-1267-dd06-05 (2010).10.1557/proc-1267-dd06-05
17. Rajput K Roy DR Structure, stability, electronic and thermoelectric properties of strontium chalcogenides Phys. E Low-Dimens. Syst. Nanostruct. 2020 119 113965 10.1016/j.physe.2020.113965
Rajput, K. & Roy, D. R. Structure, stability, electronic and thermoelectric properties of strontium chalcogenides. Phys. E Low-Dimens. Syst. Nanostruct. 119, 113965. 10.1016/j.physe.2020.113965 (2020).10.1016/j.physe.2020.113965
18. Ganose AM Jain A Robocrystallographer: Automated crystal structure text descriptions and analysis MRS Commun. 2019 9 3 874 881 10.1557/mrc.2019.94
Ganose, A. M. & Jain, A. Robocrystallographer: Automated crystal structure text descriptions and analysis. MRS Commun. 9(3), 874–881. 10.1557/mrc.2019.94 (2019).10.1557/mrc.2019.94
19. Ong SP The materials application programming interface (API): A simple, flexible and efficient API for materials data based on REpresentational State Transfer (REST) principles Comput. Mater. Sci. 2015 97 209 215 10.1016/j.commatsci.2014.10.037
Ong, S. P. et al. The materials application programming interface (API): A simple, flexible and efficient API for materials data based on REpresentational State Transfer (REST) principles. Comput. Mater. Sci. 97, 209–215. 10.1016/j.commatsci.2014.10.037 (2015).10.1016/j.commatsci.2014.10.037
20. Jain A Commentary: The Materials Project: A materials genome approach to accelerating materials innovation APL Mater. 2013 1 1 11002 10.1063/1.4812323
Jain, A. et al. Commentary: The Materials Project: A materials genome approach to accelerating materials innovation. APL Mater. 1(1), 11002. 10.1063/1.4812323 (2013).10.1063/1.4812323
21. T. M. Project. Materials Data on SrTe by Materials Project. United States. 10.17188/1194783 (2020).
22. Momma K Izumi F VESTA: A three-dimensional visualization system for electronic and structural analysis J. Appl. Crystallogr. 2008 41 3 653 658 10.1107/S0021889808012016
Momma, K. & Izumi, F. VESTA: A three-dimensional visualization system for electronic and structural analysis. J. Appl. Crystallogr. 41(3), 653–658. 10.1107/S0021889808012016 (2008).10.1107/S0021889808012016
23. Giannozzi P QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials J. Phys. Condens. Matter 2009 21 39 395502 10.1088/0953-8984/21/39/395502 21832390
Giannozzi, P. et al. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 21(39), 395502. 10.1088/0953-8984/21/39/395502 (2009).21832390 10.1088/0953-8984/21/39/395502
24. March, U. & Espresso, Q. The Quantum ESPRESSO Software Distribution, vol. 395502 (2018).
25. Patel HS Dabhi VA Vora AM Singh DK Das S Materny A To study the structural and electronic properties of TiBeO 3 using density functional theory Advances in Spectroscopy: Molecules to Materials 2019 Springer 389 395
Patel, H. S., Dabhi, V. A. & Vora, A. M. To study the structural and electronic properties of TiBeO 3 using density functional theory. In Advances in Spectroscopy: Molecules to Materials Vol. 236 (eds Singh, D. K. et al.) 389–395 (Springer, 2019).
26. Schlipf M Gygi F Optimization algorithm for the generation of ONCV pseudopotentials Comput. Phys. Commun. 2015 196 36 44 10.1016/j.cpc.2015.05.011
Schlipf, M. & Gygi, F. Optimization algorithm for the generation of ONCV pseudopotentials. Comput. Phys. Commun. 196, 36–44. 10.1016/j.cpc.2015.05.011 (2015).10.1016/j.cpc.2015.05.011
27. Perdew JP Wang Y Accurate and simple analytic representation of the electron-gas correlation energy Phys. Rev. B 1992 45 23 13244 10.1103/PhysRevB.45.13244
Perdew, J. P. & Wang, Y. Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B 45(23), 13244 (1992).10.1103/PhysRevB.45.13244
28. Patel HS Dabhi VA Vora AM Elastic constants of beryllium oxide: A first-principles investigation AIP Conf. Proc. 2020 2224 1 30006 10.1063/5.0000485
Patel, H. S., Dabhi, V. A. & Vora, A. M. Elastic constants of beryllium oxide: A first-principles investigation. AIP Conf. Proc. 2224(1), 30006 (2020).10.1063/5.0000485
29. Enkovaara J Electronic structure calculations with GPAW: A real-space implementation of the projector augmented-wave method J. Phys. Condens. Matter 2010 22 25 253202 10.1088/0953-8984/22/25/253202 21393795
Enkovaara, J. et al. Electronic structure calculations with GPAW: A real-space implementation of the projector augmented-wave method. J. Phys. Condens. Matter 22(25), 253202. 10.1088/0953-8984/22/25/253202 (2010).21393795 10.1088/0953-8984/22/25/253202
30. Heyd J Scuseria GE Ernzerhof M Hybrid functionals based on a screened Coulomb potential J. Chem. Phys. 2003 118 18 8207 8215 10.1063/1.1564060
Heyd, J., Scuseria, G. E. & Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 118(18), 8207–8215. 10.1063/1.1564060 (2003).10.1063/1.1564060
31. Patel HS Dabhi VA Vora AM Adverse effect of K-Mesh shifting in several crystal Systems: An analytical study Mater. Today Proc. 2022 10.1016/j.matpr.2022.02.599 35846849
Patel, H. S., Dabhi, V. A. & Vora, A. M. Adverse effect of K-Mesh shifting in several crystal Systems: An analytical study. Mater. Today Proc.10.1016/j.matpr.2022.02.599 (2022).35846849 10.1016/j.matpr.2022.02.599
32. Otero-de-la-Roza A Luaña V Gibbs2: A new version of the quasi-harmonic model code. I. Robust treatment of the static data Comput. Phys. Commun. 2011 182 8 1708 1720 10.1016/j.cpc.2011.04.016
Otero-de-la-Roza, A. & Luaña, V. Gibbs2: A new version of the quasi-harmonic model code. I. Robust treatment of the static data. Comput. Phys. Commun. 182(8), 1708–1720. 10.1016/j.cpc.2011.04.016 (2011).10.1016/j.cpc.2011.04.016
33. Otero-de-la-Roza A Abbasi-Pérez D Luaña V Gibbs2: A new version of the quasiharmonic model code. II. Models for solid-state thermodynamics, features and implementation Comput. Phys. Commun. 2011 182 10 2232 2248 10.1016/j.cpc.2011.05.009
Otero-de-la-Roza, A., Abbasi-Pérez, D. & Luaña, V. Gibbs2: A new version of the quasiharmonic model code. II. Models for solid-state thermodynamics, features and implementation. Comput. Phys. Commun. 182(10), 2232–2248. 10.1016/j.cpc.2011.05.009 (2011).10.1016/j.cpc.2011.05.009
34. Mostofi AA An updated version of wannier90: A tool for obtaining maximally-localised Wannier functions Comput. Phys. Commun. 2014 185 8 2309 2310 10.1016/j.cpc.2014.05.003
Mostofi, A. A. et al. An updated version of wannier90: A tool for obtaining maximally-localised Wannier functions. Comput. Phys. Commun. 185(8), 2309–2310. 10.1016/j.cpc.2014.05.003 (2014).10.1016/j.cpc.2014.05.003
35. Madsen GKH Carrete J Verstraete MJ BoltzTraP2, a program for interpolating band structures and calculating semi-classical transport coefficients Comput. Phys. Commun. 2018 231 140 145 10.1016/j.cpc.2018.05.010
Madsen, G. K. H., Carrete, J. & Verstraete, M. J. BoltzTraP2, a program for interpolating band structures and calculating semi-classical transport coefficients. Comput. Phys. Commun. 231, 140–145. 10.1016/j.cpc.2018.05.010 (2018).10.1016/j.cpc.2018.05.010
36. Labidi S Meradji H Labidi M Ghemid S Drablia S El Haj Hassan F First principles calculations of structural, electronic and thermodynamic properties of SrS, SrSe, SrTe compounds and SrS1−xSex alloy Phys. Procedia 2009 2 3 1205 1212 10.1016/j.phpro.2009.11.083
Labidi, S. et al. First principles calculations of structural, electronic and thermodynamic properties of SrS, SrSe, SrTe compounds and SrS1−xSex alloy. Phys. Procedia 2(3), 1205–1212. 10.1016/j.phpro.2009.11.083 (2009).10.1016/j.phpro.2009.11.083
37. Ren Y First principles study of electronic, phonon and elastic properties of rock-salt-phase MTe (M = Mg, Ca, Sr, Ba) Comput. Condens. Matter 2017 11 3 69 76 10.1016/j.cocom.2017.05.004
Ren, Y. First principles study of electronic, phonon and elastic properties of rock-salt-phase MTe (M = Mg, Ca, Sr, Ba). Comput. Condens. Matter 11(3), 69–76. 10.1016/j.cocom.2017.05.004 (2017).10.1016/j.cocom.2017.05.004
38. Souadkia M Bennecer B Kalarasse F Mellouki A Ab initio calculation of vibrational and thermodynamic properties of SrX (S, Se, Te) in the B1 (NaCl) and B2 (CsCl) structures Comput. Mater. Sci. 2011 50 5 1701 1710 10.1016/j.commatsci.2010.12.031
Souadkia, M., Bennecer, B., Kalarasse, F. & Mellouki, A. Ab initio calculation of vibrational and thermodynamic properties of SrX (S, Se, Te) in the B1 (NaCl) and B2 (CsCl) structures. Comput. Mater. Sci. 50(5), 1701–1710. 10.1016/j.commatsci.2010.12.031 (2011).10.1016/j.commatsci.2010.12.031
39. Bahloul B Ab initio calculations of structural, electronic, optical and thermodynamic properties of alkaline earth tellurides BaxSr1−XTe J. Phys. Chem. Solids 2014 75 3 307 314 10.1016/j.jpcs.2013.11.014
Bahloul, B. et al. Ab initio calculations of structural, electronic, optical and thermodynamic properties of alkaline earth tellurides BaxSr1−XTe. J. Phys. Chem. Solids 75(3), 307–314. 10.1016/j.jpcs.2013.11.014 (2014).10.1016/j.jpcs.2013.11.014
40. Bamba CO Accurate electronic, transport, and related properties of wurtzite beryllium oxide (w-BeO) J. Mod. Phys. 2017 08 12 1938 1949 10.4236/jmp.2017.812116
Bamba, C. O. et al. Accurate electronic, transport, and related properties of wurtzite beryllium oxide (w-BeO). J. Mod. Phys. 08(12), 1938–1949. 10.4236/jmp.2017.812116 (2017).10.4236/jmp.2017.812116
41. Lakel S Elhamra F Almi K Structural phase transition, electronic, and mechanical properties of beryllium oxide: Temperature and pressure-induced effects Phys. Status Solidi Basic Res. 2018 255 4 2 9 10.1002/pssb.201700524
Lakel, S., Elhamra, F. & Almi, K. Structural phase transition, electronic, and mechanical properties of beryllium oxide: Temperature and pressure-induced effects. Phys. Status Solidi Basic Res. 255(4), 2–9. 10.1002/pssb.201700524 (2018).10.1002/pssb.201700524
42. Mikajlo EA Dorsett HE Ford MJ Trends in the band structures of the group-I and -II oxides J. Chem. Phys. 2004 120 22 10799 10806 10.1063/1.1738635 15268107
Mikajlo, E. A., Dorsett, H. E. & Ford, M. J. Trends in the band structures of the group-I and -II oxides. J. Chem. Phys. 120(22), 10799–10806. 10.1063/1.1738635 (2004).15268107 10.1063/1.1738635
43. Patel HS Dabhi VA Vora AM First principles investigation of electronic and elastic properties of strontium oxide and strontium peroxide phases: Insights into structural stability and phase transitions Phys. Scr. 2023 98 12 125930 10.1088/1402-4896/ad0802
Patel, H. S., Dabhi, V. A. & Vora, A. M. First principles investigation of electronic and elastic properties of strontium oxide and strontium peroxide phases: Insights into structural stability and phase transitions. Phys. Scr. 98(12), 125930. 10.1088/1402-4896/ad0802 (2023).10.1088/1402-4896/ad0802
44. Patel HS Dabhi VA Vora AM Investigating electronic, optical, thermodynamic, and thermoelectric properties of SrO and SrO2 phases: A density functional theory approach ACS Omega 2023 8 45 43008 43023 10.1021/acsomega.3c06221 38024681
Patel, H. S., Dabhi, V. A. & Vora, A. M. Investigating electronic, optical, thermodynamic, and thermoelectric properties of SrO and SrO2 phases: A density functional theory approach. ACS Omega 8(45), 43008–43023. 10.1021/acsomega.3c06221 (2023).38024681 10.1021/acsomega.3c06221
45. Wei J Modified embedded atom method potential for modeling the thermodynamic properties of high thermal conductivity beryllium oxide ACS Omega 2019 4 4 6339 6346 10.1021/acsomega.9b00174 31459774
Wei, J. et al. Modified embedded atom method potential for modeling the thermodynamic properties of high thermal conductivity beryllium oxide. ACS Omega 4(4), 6339–6346. 10.1021/acsomega.9b00174 (2019).31459774 10.1021/acsomega.9b00174
46. Alipour H Hamedani A Alahyarizadeh G First-principles calculations to investigate the thermal response of the ZrC(1–x)N x ceramics at extreme conditions High Temp. Mater. Process. 2023 42 1 20220241 10.1515/htmp-2022-0241
Alipour, H., Hamedani, A. & Alahyarizadeh, G. First-principles calculations to investigate the thermal response of the ZrC(1–x)N x ceramics at extreme conditions. High Temp. Mater. Process. 42(1), 20220241. 10.1515/htmp-2022-0241 (2023).10.1515/htmp-2022-0241
47. Strauch D Semiconductors, New Data and Updates for IIa–VI Compounds with Mg and Sr (Structural Properties, Thermal and Thermodynamic Properties, and Lattice Properties) 2017 Springer
Strauch, D. Semiconductors, New Data and Updates for IIa–VI Compounds with Mg and Sr (Structural Properties, Thermal and Thermodynamic Properties, and Lattice Properties) (Springer, 2017). 10.1007/978-3-662-53620-9.
48. Xia Y Hodges JM Kanatzidis MG Lattice thermal transport in group II-alloyed PbTe Appl. Phys. Lett. 2018 112 18 193121 10.1063/1.5002587
Xia, Y., Hodges, J. M. & Kanatzidis, M. G. Lattice thermal transport in group II-alloyed PbTe. Appl. Phys. Lett. 112(18), 193121. 10.1063/1.5002587 (2018).10.1063/1.5002587
49. Madsen GKH Singh DJ BoltzTraP. A code for calculating band-structure dependent quantities Comput. Phys. Commun. 2006 175 1 67 71 10.1016/j.cpc.2006.03.007
Madsen, G. K. H. & Singh, D. J. BoltzTraP. A code for calculating band-structure dependent quantities. Comput. Phys. Commun. 175(1), 67–71. 10.1016/j.cpc.2006.03.007 (2006).10.1016/j.cpc.2006.03.007
50. Nemir D Beck J On the significance of the thermoelectric figure of merit Z J. Electron. Mater. 2010 39 9 1897 1901 10.1007/s11664-009-1060-4
Nemir, D. & Beck, J. On the significance of the thermoelectric figure of merit Z. J. Electron. Mater. 39(9), 1897–1901. 10.1007/s11664-009-1060-4 (2010).10.1007/s11664-009-1060-4
51. Kang, C., Wang, H., Bahk, J.-H. & Kim, W. Thermoelectric Materials and Devices. ResearchGate (The Royal Society of Chemistry, 2015) [Online]. https://www.researchgate.net/publication/272055575_Thermoelectric_Materials_and_Devices
52. Salik L Structural, elastic, electronic, magnetic, optical, and thermoelectric properties of the diamond-like quaternary semiconductor CuMn2InSe4 J. Supercond. Nov. Magn. 2019 33 4 1091 1102 10.1007/s10948-019-05331-1
Salik, L. et al. Structural, elastic, electronic, magnetic, optical, and thermoelectric properties of the diamond-like quaternary semiconductor CuMn2InSe4. J. Supercond. Nov. Magn. 33(4), 1091–1102. 10.1007/s10948-019-05331-1 (2019).10.1007/s10948-019-05331-1
53. Zou D Xie S Liu YY Lin J Li JY Electronic structure and thermoelectric properties of half-Heusler Zr0.5Hf0.5NiSn by first-principles calculations J. Appl. Phys. 2013 10.1063/1.4804939
Zou, D., Xie, S., Liu, Y. Y., Lin, J. & Li, J. Y. Electronic structure and thermoelectric properties of half-Heusler Zr0.5Hf0.5NiSn by first-principles calculations. J. Appl. Phys.10.1063/1.4804939 (2013).10.1063/1.4804939
