# Electronic and Transport Properties of Strained and Unstrained Ge2Sb2Te5: A DFT Investigation

Jing Tian, Weiliang Ma, Pascal Boulet, Marie-Christine Record

► To cite this version:

Jing Tian, Weiliang Ma, Pascal Boulet, Marie-Christine Record. Electronic and Transport Properties of Strained and Unstrained Ge2Sb2Te5: A DFT Investigation. Materials, 2023, 16 (14), pp.5015. 10.3390/ma16145015. [https://arxiv.org/abs/2304.01507](https://arxiv.org/abs/2304.01507)

HAL Id: hal-04527085

https://amu.hal.science/hal-04527085

Submitted on 24 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_1.png" /></figure>

HAL
open science

# Electronic and Transport Properties of Strained and Unstrained Ge2Sb2Te5: A DFT Investigation

Jing Tian, Weiliang Ma, Pascal Boulet, Marie-Christine Record

<figure><img src="image_2.png" /></figure>

## To cite this version:

Jing Tian, Weiliang Ma, Pascal Boulet, Marie-Christine Record. Electronic and Transport Properties of Strained and Unstrained Ge2Sb2Te5: A DFT Investigation. Materials, 2023, 16 (14), pp.5015. 10.3390/ma16145015. [https://arxiv.org/abs/2305.01507](https://arxiv.org/abs/2305.01507)

HAL Id: hal-04527085

https://amu.hal.science/hal-04527085v1

Submitted on 24 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Electronic and Transport Properties of Strained and Unstrained Ge₂Sb₂Te₅: A DFT Investigation

Jing Tian <fnref n="1" />,<fnref n="2" />, Weiliang Ma <fnref n="1" />,<fnref n="2" />, Pascal Boulet <fnref n="1" />,* and Marie-Christine Record <fnref n="2" /> ID

**Abstract:** In recent years, layered chalcogenides have attracted interest for their appealing thermoelectric properties. We investigated the Ge₂Sb₂Te₅ compound in two different stacking sequences, named stacking 1 (S1) and stacking 2 (S2), wherein the Ge and Sb atomic positions can be interchanged in the structure. The compound unit cell, comprising nine atoms, is made of two layers separated by a gap. We show, using the quantum theory of atoms in molecules, that the bonding across the layers has characteristics of transit region bonding, though with a close resemblance to closed-shell bonding. Both S1 and S2 are shown to bear a similar small gap. The full determination of their thermoelectric properties, including the Seebeck coefficient, electrical conductivity and electronic and lattice thermal conductivities, was carried out by solving the Boltzmann transport equation. We show that stacking 1 exhibits a larger Seebeck coefficient and smaller electrical conductivity than stacking 2, which is related to their small electronic gap difference, and that S1 is more suitable for thermoelectric application than S2. Moreover, under certain conditions of temperature and doping level, it could be possible to use S1-Ge₂Sb₂Te₅ as both a *p* and *n* leg in a thermoelectric converter. Under biaxial, tensile and compressive strains, we observe that the thermoelectric properties are improved for both S1 and S2. Furthermore, the increase in the power factor of S1 in the cross-plane direction, namely perpendicular to the gap between the layers, shows that strains can counteract the electronic transport hindrance due to the gap.

**Keywords:** layered chalcogenides; thermoelectrics; strains; density functional theory; Boltzmann transport

## 1. Introduction

The rise in electrical energy needs over the decades calls for conducting deep investigations to improve storage and conversion systems. The discovery of efficient materials is one of the cornerstones that should allow for meeting the challenge of energy demand. Among the means to produce electrical power, both photovoltaic and thermoelectric conversions have regained strong interest. But both suffer from low conversion efficiency so far. Heat-to-electricity conversion has been known for more than two centuries, as this phenomenon was discovered, but remained unexplained at that time, in 1794 by Volta [1] and then rediscovered in 1821 by Seebeck [2], though erroneously interpreted by him as a manifestation of a magnetic field. Ioffe's work in the early 1960s [3] has done much to popularize the use of semiconductors as potential thermoelectric materials, and thermoelectric modules have found a niche usage in spacecrafts [4]. More recently, thermoelectric materials have been implemented for heat-to-electricity conversion in cars' exhaust-gas systems [5]. However, thermoelectricity has hardly found its way to large-scale deployment due to the lack of semiconductor materials with a figure of merit ZT greater than 2. One reason for this is that the figure of merit combines properties that counteract one another, i.e., the Seebeck coefficient and electrical conductivity that both appear in the

<figure><img src="image_2.png" /></figure>

### check for updates

**Citation:** Tian, J.; Ma, W.; Boulet, P.; Record, M.-C. Electronic and Transport Properties of Strained and Unstrained Ge₂Sb₂Te₅: A DFT Investigation. *Materials* **2023**, *16*, 5015. https://doi.org/10.3390/ma16145015

Academic Editor: Nektarios
N. Lathiotakis

Received: 20 June 2023

Revised: 9 July 2023

Accepted: 11 July 2023

Published: 15 July 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ MADIREL, Department of Chemistry, CNRS, Aix-Marseille University, 13013 Marseille, France; jing.tian@etu.univ-amu.fr (J.T.); weiliang.ma@etu.univ-amu.fr (W.M.)</footnote>

<footnote n="2">² IM2NP, Department of Chemistry, CNRS, Aix-Marseille University, 13013 Marseille, France; m-c.record@univ-amu.fr</footnote>

<footnote>* Correspondence: pascal.boulet@univ-amu.fr</footnote>
Article

# Electronic and Transport Properties of Strained and Unstrained Ge₂Sb₂Te₅: A DFT Investigation

Jing Tian <fnref n="1" />,<fnref n="2" />, Weiliang Ma <fnref n="1" />,<fnref n="2" />, Pascal Boulet <fnref n="1" />,* and Marie-Christine Record <fnref n="2" /> ID

**Abstract:** In recent years, layered chalcogenides have attracted interest for their appealing thermoelectric properties. We investigated the Ge₂Sb₂Te₅ compound in two different stacking sequences, named stacking 1 (S1) and stacking 2 (S2), wherein the Ge and Sb atomic positions can be interchanged in the structure. The compound unit cell, comprising nine atoms, is made of two layers separated by a gap. We show, using the quantum theory of atoms in molecules, that the bonding across the layers has characteristics of transit region bonding, though with a close resemblance to closed-shell bonding. Both S1 and S2 are shown to bear a similar small gap. The full determination of their thermoelectric properties, including the Seebeck coefficient, electrical conductivity and electronic and lattice thermal conductivities, was carried out by solving the Boltzmann transport equation. We show that stacking 1 exhibits a larger Seebeck coefficient and smaller electrical conductivity than stacking 2, which is related to their small electronic gap difference, and that S1 is more suitable for thermoelectric application than S2. Moreover, under certain conditions of temperature and doping level, it could be possible to use S1-Ge₂Sb₂Te₅ as both a *p* and *n* leg in a thermoelectric converter. Under biaxial, tensile and compressive strains, we observe that the thermoelectric properties are improved for both S1 and S2. Furthermore, the increase in the power factor of S1 in the cross-plane direction, namely perpendicular to the gap between the layers, shows that strains can counteract the electronic transport hindrance due to the gap.

**Keywords:** layered chalcogenides; thermoelectrics; strains; density functional theory; Boltzmann transport

## 1. Introduction

The rise in electrical energy needs over the decades calls for conducting deep investigations to improve storage and conversion systems. The discovery of efficient materials is one of the cornerstones that should allow for meeting the challenge of energy demand. Among the means to produce electrical power, both photovoltaic and thermoelectric conversions have regained strong interest. But both suffer from low conversion efficiency so far. Heat-to-electricity conversion has been known for more than two centuries, as this phenomenon was discovered, but remained unexplained at that time, in 1794 by Volta [1] and then rediscovered in 1821 by Seebeck [2], though erroneously interpreted by him as a manifestation of a magnetic field. Ioffe's work in the early 1960s [3] has done much to popularize the use of semiconductors as potential thermoelectric materials, and thermoelectric modules have found a niche usage in spacecrafts [4]. More recently, thermoelectric materials have been implemented for heat-to-electricity conversion in cars' exhaust-gas systems [5]. However, thermoelectricity has hardly found its way to large-scale deployment due to the lack of semiconductor materials with a figure of merit ZT greater than 2. One reason for this is that the figure of merit combines properties that counteract one another, i.e., the Seebeck coefficient and electrical conductivity that both appear in the

<figure><img src="image_4.png" /></figure>

### check for updates

**Citation:** Tian, J.; Ma, W.; Boulet, P.; Record, M.-C. Electronic and Transport Properties of Strained and Unstrained Ge₂Sb₂Te₅: A DFT Investigation. *Materials* **2023**, *16*, 5015. https://doi.org/10.3390/ma16145015

Academic Editor: Nektarios
N. Lathiotakis

Received: 20 June 2023

Revised: 9 July 2023

Accepted: 11 July 2023

Published: 15 July 2023

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ MADIREL, Department of Chemistry, CNRS, Aix-Marseille University, 13013 Marseille, France; jing.tian@etu.univ-amu.fr (J.T.); weiliang.ma@etu.univ-amu.fr (W.M.)</footnote>

<footnote n="2">² IM2NP, Department of Chemistry, CNRS, Aix-Marseille University, 13013 Marseille, France; m-c.record@univ-amu.fr</footnote>

<footnote>* Correspondence: pascal.boulet@univ-amu.fr</footnote>
numerator of ZT, and the electrical conductivity and electronic thermal conductivity that contribute to the numerator and denominator, respectively. The ZT expression then reads $\frac{S^2\sigma}{\kappa_e+\kappa_l}T$, where S is the Seebeck coefficient, σ the electrical conductivity, κe the electronic thermal conductivity, κl the lattice thermal conductivity and T the temperature. Hence, a compromise must be found to optimize materials ZT. Various families of material compounds have been shown to bear high thermoelectric (TE) properties, e.g., skutterudites, silicides, Zintl compounds and clathrates [6,7]. In the present work, the Ge₂Sb₂Te₅ layered chalcogenide compound was investigated. Although this compound is well known as a phase change material [8], it shares with thermoelectrics chemical elements that are known to bring interesting TE properties, which legitimates the investigation of Ge₂Sb₂Te₅ for its potential TE applications. In addition, the layered nature of Ge₂Sb₂Te₅ is a beneficial feature for TE performance as it allows for decreasing lattice thermal conductivity. Very recently, Miao et al. [9] investigated Ge₂Sb₂Te₅ using density functional theory approaches and found that, under certain conditions of temperature and carrier doping level, the figure of merit can reach 1.5. In their work, the ZT values were predicted by combining theoretical and experimental data. Ge₂Sb₂Te₅ has also been investigated experimentally for its thermoelectric properties. It is recognized that this compound bears a decent ZT value, though it has a high electrical conductivity (due to high hole concentration) and hence a low Seebeck coefficient [10]. The substitution of In for Ge [10], Se for Te [11] and S/Se for Te [12] has shown the possibility of both increasing the Seebeck coefficient through the modulation of the density of states and decreasing the lattice thermal conductivity through mass disorder improvement, hence increasing ZT by 30% to almost 50%, reaching the values of 0.78 in Ge₁.₈₅In₀.₁₅Sb₂Te₅ at 700 K [10], 0.41 in Ge₂Sb₂Te₃.₅Se₁.₅ at 703 K [11] and 0.74 in Ge₂Sb₂Te₄.₉S₀.₁ at 800 K [12]. According to these ZT values, it appears that Ge₂Sb₂Te₅ bears similar thermoelectric performances to other well-known TE materials (skutterudites, Zintl phases, etc.).

In the present work, a fully ab initio approach was used to determine the TE properties. To date, the improvement in the TE properties of Ge₂Sb₂Te₅ has been investigated using chemical substitution approaches, which have proven their efficiency, but not via induced strains on the material. Hence, we propose in our investigation to study the effects of strains on Ge₂Sb₂Te₅ and demonstrate the possibility of improving the TE properties through the band engineering approach.

## 2. Computational Approaches

Density functional theory [13,14] (DFT) was used throughout to calculate the electronic properties of Ge₂Sb₂Te₅. Local and generalized-gradient approximations were implemented with various exchange-correlation functionals (PZ [15], WC [16], PBE [17], PBEsol [18]) as well as the hybrid HSE06 [19] one to investigate its structural properties. For the subsequent TE property calculations, the WC functional was selected. The spin-orbit coupling was accounted for in these calculations for all the atoms. The Brillouin zone (BZ) was sampled with the *k*-meshes 12 × 12 × 2, 18 × 18 × 4 and 64 × 64 × 14 for structural optimizations, SCF calculation and transport property calculation, respectively. The total energy and atomic force convergence thresholds were defined as 0.136 meV and 0.257 meV.Å⁻¹, respectively. The cutoff energy for the core-valence separation was fixed at −6.0 Ry. The *R*<sub>MT</sub>*K*<sub>max</sub> value was set to 9.0 where *R*<sub>MT</sub> is the minimum LAPW radius and *K*<sub>max</sub> is the largest *k* vector in the Brillouin zone for the plane-wave cutoff. The radius of muffin tin (RMT) used for Ge, Sb and Te atoms in this work was set to 2.5 Å. The Wien2k package was used for these calculations [20,21].

The transport properties (Seebeck coefficient, electrical conductivity and electronic thermal conductivity) were determined using the Boltzmann transport equation (BTE) approach, as implemented in the BoltzTrap2 program [22]. The relaxation time of the electrons was calculated by applying the deformation potential theory [23,24].

The dynamic properties were calculated by means of the DFPT method by combining the Quantum-ESPRESSO package [25] and the Phono3py program [26,27]. A plane-wave
numerator of ZT, and the electrical conductivity and electronic thermal conductivity that contribute to the numerator and denominator, respectively. The ZT expression then reads $\frac{S^2\sigma}{\kappa_e+\kappa_l}T$, where S is the Seebeck coefficient, σ the electrical conductivity, κe the electronic thermal conductivity, κl the lattice thermal conductivity and T the temperature. Hence, a compromise must be found to optimize materials ZT. Various families of material compounds have been shown to bear high thermoelectric (TE) properties, e.g., skutterudites, silicides, Zintl compounds and clathrates [6,7]. In the present work, the Ge₂Sb₂Te₅ layered chalcogenide compound was investigated. Although this compound is well known as a phase change material [8], it shares with thermoelectrics chemical elements that are known to bring interesting TE properties, which legitimates the investigation of Ge₂Sb₂Te₅ for its potential TE applications. In addition, the layered nature of Ge₂Sb₂Te₅ is a beneficial feature for TE performance as it allows for decreasing lattice thermal conductivity. Very recently, Miao et al. [9] investigated Ge₂Sb₂Te₅ using density functional theory approaches and found that, under certain conditions of temperature and carrier doping level, the figure of merit can reach 1.5. In their work, the ZT values were predicted by combining theoretical and experimental data. Ge₂Sb₂Te₅ has also been investigated experimentally for its thermoelectric properties. It is recognized that this compound bears a decent ZT value, though it has a high electrical conductivity (due to high hole concentration) and hence a low Seebeck coefficient [10]. The substitution of In for Ge [10], Se for Te [11] and S/Se for Te [12] has shown the possibility of both increasing the Seebeck coefficient through the modulation of the density of states and decreasing the lattice thermal conductivity through mass disorder improvement, hence increasing ZT by 30% to almost 50%, reaching the values of 0.78 in Ge₁.₈₅In₀.₁₅Sb₂Te₅ at 700 K [10], 0.41 in Ge₂Sb₂Te₃.₅Se₁.₅ at 703 K [11] and 0.74 in Ge₂Sb₂Te₄.₉S₀.₁ at 800 K [12]. According to these ZT values, it appears that Ge₂Sb₂Te₅ bears similar thermoelectric performances to other well-known TE materials (skutterudites, Zintl phases, etc.).

In the present work, a fully ab initio approach was used to determine the TE properties. To date, the improvement in the TE properties of Ge₂Sb₂Te₅ has been investigated using chemical substitution approaches, which have proven their efficiency, but not via induced strains on the material. Hence, we propose in our investigation to study the effects of strains on Ge₂Sb₂Te₅ and demonstrate the possibility of improving the TE properties through the band engineering approach.

## 2. Computational Approaches

Density functional theory [13,14] (DFT) was used throughout to calculate the electronic properties of Ge₂Sb₂Te₅. Local and generalized-gradient approximations were implemented with various exchange-correlation functionals (PZ [15], WC [16], PBE [17], PBEsol [18]) as well as the hybrid HSE06 [19] one to investigate its structural properties. For the subsequent TE property calculations, the WC functional was selected. The spin-orbit coupling was accounted for in these calculations for all the atoms. The Brillouin zone (BZ) was sampled with the *k*-meshes 12 × 12 × 2, 18 × 18 × 4 and 64 × 64 × 14 for structural optimizations, SCF calculation and transport property calculation, respectively. The total energy and atomic force convergence thresholds were defined as 0.136 meV and 0.257 meV.Å⁻¹, respectively. The cutoff energy for the core-valence separation was fixed at −6.0 Ry. The *R*<sub>MT</sub>*K*<sub>max</sub> value was set to 9.0 where *R*<sub>MT</sub> is the minimum LAPW radius and *K*<sub>max</sub> is the largest *k* vector in the Brillouin zone for the plane-wave cutoff. The radius of muffin tin (RMT) used for Ge, Sb and Te atoms in this work was set to 2.5 Å. The Wien2k package was used for these calculations [20,21].

The transport properties (Seebeck coefficient, electrical conductivity and electronic thermal conductivity) were determined using the Boltzmann transport equation (BTE) approach, as implemented in the BoltzTrap2 program [22]. The relaxation time of the electrons was calculated by applying the deformation potential theory [23,24].

The dynamic properties were calculated by means of the DFPT method by combining the Quantum-ESPRESSO package [25] and the Phono3py program [26,27]. A plane-wave
energy cutoff of 70 Ry (952 eV) was employed. Total energies were minimized, with a convergence criterium of $10^{-7}$ Ry and a total force threshold of $10^{-4}$ Ry/bohr. A supercell of $2 \times 2 \times 2$ was considered, which consists of a total of 72 atoms for Ge$_2$Sb$_2$Te$_5$, with $4 \times 4 \times 1$ q-mesh sampling. In subsequent postprocessing, phonon lifetimes were sampled using a finer $19 \times 19 \times 5$ mesh.

### 3. Results and Discussion

In this section, we present and interpret the results of this investigation concerning the structural, electronic and transport features of the stacking 1 and stacking 2 Ge$_2$Sb$_2$Te$_5$ compounds, both with and without applied strains.

#### 3.1. Unstrained Ge$_2$Sb$_2$Te$_5$ Compounds

##### (a) Electronic Transport

The layered Ge$_2$Sb$_2$Te$_5$ compounds crystallize in the rhombohedral Bravais system and belong to the space group number 164 (P$\bar{3}$m1). A hexagonal representation (conventional cell) of the stacking 1 and stacking 2 sequences of Ge$_2$Sb$_2$Te$_5$ is depicted in Figure 1. From experimental crystallographic data, the c/a ratio amounts to 4.080 [28] and 4.299 [29] for stacking 1 and 4.080 [28] and 4.038 [30] for stacking 2 (Table 1). Averaged over all the functionals used, our results for the c/a ratio deviate by 0.6% and 5.6% for stacking 1 compared with Ref. [28] and Ref. [29], respectively, and by 2.3% and 4.3% for stacking 2, respectively. The lattice parameters are reasonably well reproduced by the functionals (Table 1), with average deviations for stacking 1 of 0.9% and 0.8% on *a* with respect to Ref. [28] and Ref. [29] and of 1.2% and 5.6% on *c*. For stacking 2, these deviations amount to 0.7% and 0.8% on *a* and 2.6% and 3.4% on *c*. As it turns out, the van der Waals-corrected PBE functional overestimates both the *a* and *c* parameters; therefore, considering the good performance of WC to reproduce the lattice parameters, this functional was chosen as the GGA one for further investigations. Nonetheless, our estimated band gap energies, including spin-orbit coupling or its absence, are seemingly underestimated with respect to the experimental findings, a feature that holds true for the hybrid HSE06 functional (Table 2). In any case, stacking 1 (S1) is found to be a small gap semiconductor (SC) while stacking 2 (S2) is a metal with the WC functional, whereas HSE06 describes S1 as small gap SC either with or without SOC, and S2 is either a metal without SOC or a small gap SC with SOC. The SOC interaction noticeably decreases the band gap energy of S1, but one does not observe specific features of SOC such as band splitting (Figures 2 and 3). It is obvious that both the SOC and the exact exchange bring significant changes in the features of the electronic structures of S1 and S2 (Table 3). Assuming that HSE06-SOC is the most reliable approximation, one can summarize the results as follows: (1) both S1 and S2 are small gap SCs, (2) the gap of S1-Ge$_2$Sb$_2$Te$_5$ is a direct one at $\Gamma$ and that of S2-Ge$_2$Sb$_2$Te$_5$ is an indirect one and (3) the energy gap of S1 is larger than that of S2. We note in passing that the WC-SOC approach also yields a direct band gap, though at a different *k*-point in the Brillouin zone. As transport property calculations are quite demanding in regard to Brillouin zone *k*-point sampling, the investigation of thermoelectric properties is intractable with the HSE06 functional, and we resorted to the WC+SOC one.

<table><caption>**Table 1.** Lattice parameters (in Å) and atomic positions (in crystallographic coordinates) of the stacking 1 and stacking 2 Ge$_2$Sb$_2$Te$_5$. Only the z crystallographic coordinates are provided, as the *x* and *y* ones are fixed by the crystal symmetry.</caption><thead><tr><th></th><th></th><th>LDA</th><th>WC</th><th>PBE</th><th>PBEsol</th><th>PBE-vdw</th><th>Exp.</th></tr></thead><tbody><tr><td rowspan="5">Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub><br>S1</td><td><i>a</i></td><td>4.20</td><td>4.24</td><td>4.30</td><td>4.23</td><td>4.28</td><td>4.225 [28] 4.25 [29]</td></tr><tr><td><i>c</i></td><td>17.09</td><td>16.92</td><td>17.51</td><td>17.23</td><td>17.47</td><td>17.239 [28] 18.27 [29]</td></tr><tr><td>Ge</td><td>0.0998</td><td>0.1006</td><td>0.0988</td><td>0.0995</td><td>0.0997</td><td></td></tr><tr><td>Te2</td><td>0.1974</td><td>0.1990</td><td>0.1953</td><td>0.1968</td><td>0.1971</td><td></td></tr><tr><td>Sb</td><td>0.3153</td><td>0.3178</td><td>0.3113</td><td>0.3141</td><td>0.3142</td><td></td></tr><tr><td></td><td>Te3</td><td>0.4161</td><td>0.4198</td><td>0.4091</td><td>0.4143</td><td>0.4132</td><td></td></tr></tbody></table>
energy cutoff of 70 Ry (952 eV) was employed. Total energies were minimized, with a convergence criterium of $10^{-7}$ Ry and a total force threshold of $10^{-4}$ Ry/bohr. A supercell of $2 \times 2 \times 2$ was considered, which consists of a total of 72 atoms for Ge$_2$Sb$_2$Te$_5$, with $4 \times 4 \times 1$ q-mesh sampling. In subsequent postprocessing, phonon lifetimes were sampled using a finer $19 \times 19 \times 5$ mesh.

### 3. Results and Discussion

In this section, we present and interpret the results of this investigation concerning the structural, electronic and transport features of the stacking 1 and stacking 2 Ge$_2$Sb$_2$Te$_5$ compounds, both with and without applied strains.

#### 3.1. *Unstrained Ge$_2$Sb$_2$Te$_5$ Compounds*

##### (a) Electronic Transport

The layered Ge$_2$Sb$_2$Te$_5$ compounds crystallize in the rhombohedral Bravais system and belong to the space group number 164 (P$\bar{3}$m1). A hexagonal representation (conventional cell) of the stacking 1 and stacking 2 sequences of Ge$_2$Sb$_2$Te$_5$ is depicted in Figure 1. From experimental crystallographic data, the c/a ratio amounts to 4.080 [28] and 4.299 [29] for stacking 1 and 4.080 [28] and 4.038 [30] for stacking 2 (Table 1). Averaged over all the functionals used, our results for the c/a ratio deviate by 0.6% and 5.6% for stacking 1 compared with Ref. [28] and Ref. [29], respectively, and by 2.3% and 4.3% for stacking 2, respectively. The lattice parameters are reasonably well reproduced by the functionals (Table 1), with average deviations for stacking 1 of 0.9% and 0.8% on *a* with respect to Ref. [28] and Ref. [29] and of 1.2% and 5.6% on *c*. For stacking 2, these deviations amount to 0.7% and 0.8% on *a* and 2.6% and 3.4% on *c*. As it turns out, the van der Waals-corrected PBE functional overestimates both the *a* and *c* parameters; therefore, considering the good performance of WC to reproduce the lattice parameters, this functional was chosen as the GGA one for further investigations. Nonetheless, our estimated band gap energies, including spin-orbit coupling or its absence, are seemingly underestimated with respect to the experimental findings, a feature that holds true for the hybrid HSE06 functional (Table 2). In any case, stacking 1 (S1) is found to be a small gap semiconductor (SC) while stacking 2 (S2) is a metal with the WC functional, whereas HSE06 describes S1 as small gap SC either with or without SOC, and S2 is either a metal without SOC or a small gap SC with SOC. The SOC interaction noticeably decreases the band gap energy of S1, but one does not observe specific features of SOC such as band splitting (Figures 2 and 3). It is obvious that both the SOC and the exact exchange bring significant changes in the features of the electronic structures of S1 and S2 (Table 3). Assuming that HSE06-SOC is the most reliable approximation, one can summarize the results as follows: (1) both S1 and S2 are small gap SCs, (2) the gap of S1-Ge$_2$Sb$_2$Te$_5$ is a direct one at $\Gamma$ and that of S2-Ge$_2$Sb$_2$Te$_5$ is an indirect one and (3) the energy gap of S1 is larger than that of S2. We note in passing that the WC-SOC approach also yields a direct band gap, though at a different *k*-point in the Brillouin zone. As transport property calculations are quite demanding in regard to Brillouin zone *k*-point sampling, the investigation of thermoelectric properties is intractable with the HSE06 functional, and we resorted to the WC+SOC one.

<table><caption>**Table 1.** Lattice parameters (in Å) and atomic positions (in crystallographic coordinates) of the stacking 1 and stacking 2 Ge$_2$Sb$_2$Te$_5$. Only the z crystallographic coordinates are provided, as the *x* and *y* ones are fixed by the crystal symmetry.</caption><thead><tr><th></th><th></th><th>LDA</th><th>WC</th><th>PBE</th><th>PBEsol</th><th>PBE-vdw</th><th>Exp.</th></tr></thead><tbody><tr><td rowspan="5">Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub><br>S1</td><td><i>a</i></td><td>4.20</td><td>4.24</td><td>4.30</td><td>4.23</td><td>4.28</td><td>4.225 [28] 4.25 [29]</td></tr><tr><td><i>c</i></td><td>17.09</td><td>16.92</td><td>17.51</td><td>17.23</td><td>17.47</td><td>17.239 [28] 18.27 [29]</td></tr><tr><td>Ge</td><td>0.0998</td><td>0.1006</td><td>0.0988</td><td>0.0995</td><td>0.0997</td><td></td></tr><tr><td>Te2</td><td>0.1974</td><td>0.1990</td><td>0.1953</td><td>0.1968</td><td>0.1971</td><td></td></tr><tr><td>Sb</td><td>0.3153</td><td>0.3178</td><td>0.3113</td><td>0.3141</td><td>0.3142</td><td></td></tr><tr><td></td><td>Te3</td><td>0.4161</td><td>0.4198</td><td>0.4091</td><td>0.4143</td><td>0.4132</td><td></td></tr></tbody></table>
<table><caption>**Table 1.** *Cont.*</caption><thead><tr><th></th><th></th><th>LDA</th><th>WC</th><th>PBE</th><th>PBEsol</th><th>PBE-vdw</th><th>Exp.</th></tr></thead><tbody><tr><td rowspan="5">Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub><br>S2</td><td><i>a</i></td><td>4.18</td><td>4.22</td><td>4.27</td><td>4.21</td><td>4.25</td><td>4.225 [28] 4.20 [30]</td></tr><tr><td><i>c</i></td><td>17.49</td><td>17.13</td><td>18.01</td><td>17.65</td><td>17.90</td><td>17.239 [28] 16.96 [30]</td></tr><tr><td>Ge</td><td>0.3353</td><td>0.3387</td><td>0.3329</td><td>0.3337</td><td>0.3361</td><td></td></tr><tr><td>Te2</td><td>0.2151</td><td>0.2188</td><td>0.2094</td><td>0.2139</td><td>0.2114</td><td></td></tr><tr><td>Sb</td><td>0.1150</td><td>0.1167</td><td>0.1124</td><td>0.1143</td><td>0.1135</td><td></td></tr><tr><td></td><td>Te3</td><td>0.4161</td><td>0.4213</td><td>0.4116</td><td>0.4143</td><td>0.4158</td><td></td></tr></tbody></table>

<table><caption>**Table 2.** Band gap energies (eV) of stacking 1 and stacking 2 Ge₂Sb₂Te₅ compounds with and without spin-orbit coupling (SOC and non-SOC).</caption><thead><tr><th></th><th>WC</th><th>HSE06</th><th>Exp.</th><th>Other Calculations</th></tr></thead><tbody><tr><td>S1 (Non-SOC)</td><td>0.147</td><td>0.384</td><td>0.5 [31–33]</td><td>0.26 (PBE) [34], 0.3 (PBE-D2) [35], 0.49 (HSE06-D2) [35]</td></tr><tr><td>S1 (SOC)</td><td>0.039</td><td>0.090</td><td></td><td></td></tr><tr><td>S2 (Non-SOC)</td><td>0</td><td>0.122</td><td></td><td>0 (PBE) [34]</td></tr><tr><td>S2 (SOC)</td><td>0</td><td>0.086</td><td></td><td></td></tr></tbody></table>

<table><caption>**Table 3.** Effects of the spin-orbit coupling (SOC) and exact exchange interaction, through the HSE06 functional, on the band structures of stacking 1 and stacking 2 of Ge₂Sb₂Te₅ as observed from Figures 2 and 3.</caption><thead><tr><th rowspan="2"></th><th colspan="2">WC Band Gap Location or Closure</th><th colspan="2">HSE06 Band Gap Location or Closure</th></tr><tr><th>Without SOC</th><th>With SOC</th><th>Without SOC</th><th>With SOC</th></tr></thead><tbody><tr><td>S1</td><td>Opened, Γ → Γ-K path</td><td>Opened, Γ-A → Γ-A paths</td><td>Opened, Γ → Γ-K path</td><td>Opened, Γ → Γ</td></tr><tr><td>S2</td><td>Closed, A-K and Γ-K paths</td><td>Closed, Γ-M and Γ-K paths</td><td>Closed, A-K, Γ-K and Γ-M paths</td><td>Opened, Γ-M → Γ-K paths</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Hexagonal conventional cell representation of stacking 1 (a) and stacking 2 (b) polymorphs of Ge₂Sb₂Te₅.</figcaption></figure>
<table><caption>**Table 1.** *Cont.*</caption><thead><tr><th></th><th></th><th>LDA</th><th>WC</th><th>PBE</th><th>PBEsol</th><th>PBE-vdw</th><th>Exp.</th></tr></thead><tbody><tr><td rowspan="5">Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub><br>S2</td><td><i>a</i></td><td>4.18</td><td>4.22</td><td>4.27</td><td>4.21</td><td>4.25</td><td>4.225 [28] 4.20 [30]</td></tr><tr><td><i>c</i></td><td>17.49</td><td>17.13</td><td>18.01</td><td>17.65</td><td>17.90</td><td>17.239 [28] 16.96 [30]</td></tr><tr><td>Ge</td><td>0.3353</td><td>0.3387</td><td>0.3329</td><td>0.3337</td><td>0.3361</td><td></td></tr><tr><td>Te2</td><td>0.2151</td><td>0.2188</td><td>0.2094</td><td>0.2139</td><td>0.2114</td><td></td></tr><tr><td>Sb</td><td>0.1150</td><td>0.1167</td><td>0.1124</td><td>0.1143</td><td>0.1135</td><td></td></tr><tr><td></td><td>Te3</td><td>0.4161</td><td>0.4213</td><td>0.4116</td><td>0.4143</td><td>0.4158</td><td></td></tr></tbody></table>

<table><caption>**Table 2.** Band gap energies (eV) of stacking 1 and stacking 2 Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> compounds with and without spin-orbit coupling (SOC and non-SOC).</caption><thead><tr><th></th><th>WC</th><th>HSE06</th><th>Exp.</th><th>Other Calculations</th></tr></thead><tbody><tr><td>S1 (Non-SOC)</td><td>0.147</td><td>0.384</td><td>0.5 [31–33]</td><td>0.26 (PBE) [34], 0.3 (PBE-D2) [35], 0.49 (HSE06-D2) [35]</td></tr><tr><td>S1 (SOC)</td><td>0.039</td><td>0.090</td><td></td><td></td></tr><tr><td>S2 (Non-SOC)</td><td>0</td><td>0.122</td><td></td><td>0 (PBE) [34]</td></tr><tr><td>S2 (SOC)</td><td>0</td><td>0.086</td><td></td><td></td></tr></tbody></table>

<table><caption>**Table 3.** Effects of the spin-orbit coupling (SOC) and exact exchange interaction, through the HSE06 functional, on the band structures of stacking 1 and stacking 2 of Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> as observed from Figures 2 and 3.</caption><thead><tr><th rowspan="2"></th><th colspan="2">WC Band Gap Location or Closure</th><th colspan="2">HSE06 Band Gap Location or Closure</th></tr><tr><th>Without SOC</th><th>With SOC</th><th>Without SOC</th><th>With SOC</th></tr></thead><tbody><tr><td>S1</td><td>Opened, Γ → Γ-K path</td><td>Opened, Γ-A → Γ-A paths</td><td>Opened, Γ → Γ-K path</td><td>Opened, Γ → Γ</td></tr><tr><td>S2</td><td>Closed, A-K and Γ-K paths</td><td>Closed, Γ-M and Γ-K paths</td><td>Closed, A-K, Γ-K and Γ-M paths</td><td>Opened, Γ-M → Γ-K paths</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 1.** Hexagonal conventional cell representation of stacking 1 (a) and stacking 2 (b) polymorphs of Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub>.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Electronic band structures in the Brillouin zone of the Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> compound with stacking 1 (a) and stacking 2 (b) calculated with the WC functional with spin-orbit coupling (red) and without (blue). The energies are scaled using the Fermi energy (E<sub>f</sub>).</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Electronic band structures in the Brillouin zone of the Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> compound with stacking 1 (a) and stacking 2 (b) calculated with the HSE06 hybrid functional with spin-orbit coupling (red) and without (blue). The energies are scaled using the Fermi energy (E<sub>f</sub>).</figcaption></figure>

The holes' and electrons' effective masses are, among other properties, important characteristics for electronic transport, and relate to the inverse curvature of the highest valence band and the lowest conduction band, respectively. The curvature of the valence and conduction bands of S1 in the *ab* plane is smaller than that in the *c* direction (Table 4), leading to larger mobility of the charge carriers. In the simple picture of the Drude model, the electrical conductivity ($\sigma = ne^2\tau/m^*$) should then be higher in the basal plane of the structure. The same feature is observed for stacking S2 of Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub>.

<table><caption>**Table 4.** Relative effective masses ($m^*$ in unit of the electron rest mass) of holes and electrons and potential deformation energies (E<sub>d</sub>), calculated in the *c*-axis direction and *ab* plane for stacking 1 and stacking 2 of Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> with the WC+SOC approach.</caption><thead><tr><th rowspan="2">Structure</th><th rowspan="2">Charge</th><th colspan="2">m*</th><th colspan="2">E<sub>d</sub> (eV)</th></tr><tr><th>c</th><th>a,b</th><th>c</th><th>a,b</th></tr></thead><tbody><tr><td>Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> S1</td><td>holes</td><td>-0.0974</td><td>-0.0187</td><td>7.01</td><td>13.39</td></tr><tr><td>Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> S1</td><td>electrons</td><td>0.0797</td><td>0.0175</td><td>10.09</td><td>14.88</td></tr><tr><td>Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> S2</td><td>holes</td><td>-0.0215</td><td>-0.0065</td><td>6.92</td><td>14.16</td></tr><tr><td>Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> S2</td><td>electrons</td><td>0.0580</td><td>0.0448</td><td>9.52</td><td>12.66</td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Figure 2.** Electronic band structures in the Brillouin zone of the Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> compound with stacking 1 (a) and stacking 2 (b) calculated with the WC functional with spin-orbit coupling (red) and without (blue). The energies are scaled using the Fermi energy (E<sub>f</sub>).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Electronic band structures in the Brillouin zone of the Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> compound with stacking 1 (a) and stacking 2 (b) calculated with the HSE06 hybrid functional with spin-orbit coupling (red) and without (blue). The energies are scaled using the Fermi energy (E<sub>f</sub>).</figcaption></figure>

The holes' and electrons' effective masses are, among other properties, important characteristics for electronic transport, and relate to the inverse curvature of the highest valence band and the lowest conduction band, respectively. The curvature of the valence and conduction bands of S1 in the *ab* plane is smaller than that in the *c* direction (Table 4), leading to larger mobility of the charge carriers. In the simple picture of the Drude model, the electrical conductivity ($\sigma = ne^2\tau/m^*$) should then be higher in the basal plane of the structure. The same feature is observed for stacking S2 of Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub>.

<table><caption>**Table 4.** Relative effective masses ($m^*$ in unit of the electron rest mass) of holes and electrons and potential deformation energies (E<sub>d</sub>), calculated in the *c*-axis direction and *ab* plane for stacking 1 and stacking 2 of Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> with the WC+SOC approach.</caption><thead><tr><th rowspan="2">Structure</th><th rowspan="2">Charge</th><th colspan="2"><i>m</i><sup>*</sup></th><th colspan="2"><i>E</i><sub>d</sub> (eV)</th></tr><tr><th><i>c</i></th><th><i>a,b</i></th><th><i>c</i></th><th><i>a,b</i></th></tr></thead><tbody><tr><td>Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> S1</td><td>holes</td><td>-0.0974</td><td>-0.0187</td><td>7.01</td><td>13.39</td></tr><tr><td>Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> S1</td><td>electrons</td><td>0.0797</td><td>0.0175</td><td>10.09</td><td>14.88</td></tr><tr><td>Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> S2</td><td>holes</td><td>-0.0215</td><td>-0.0065</td><td>6.92</td><td>14.16</td></tr><tr><td>Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> S2</td><td>electrons</td><td>0.0580</td><td>0.0448</td><td>9.52</td><td>12.66</td></tr></tbody></table>
The bonding structures of S1- and S2-Ge₂Sb₂Te₅ can be characterized by the electron density gradient (Figure 4). The zero-gradient lines delineating the so-called atomic basins allow for visualizing the atomic volumes, which are noticeably larger for Te and Ge than for Sb and are far from spherical in shape. The line of maximum gradient joining two atoms and passing through a bond critical point characterizes a bond path. In this respect, the Te atoms next to and on either side of the interlayer space are connected by a bond path (BCP, labeled b1 in Figure 4). The interlayer interaction is characterized by a low density, a slightly negative total energy density and a negative, though small, bond degree (H/ρ) (Table 5). The positive value of the electron density Laplacian (∇²ρ) underscores the tendency of the electron density to escape from the BCP region. Strictly speaking, the interlayer bonding interaction corresponds to transit region bonding, as defined by Espinosa [36], with the characteristics 1<|V|/G<2, H/ρ<0 and ∇²ρ>0, although this interaction is very close to that corresponding to the closed-shell region, which includes the van der Waals interactions. All the other bond interactions, labeled b2 to b5, are characteristic of the transit region bonding interaction (Table 5). The variety of bonds in S1- and S2-Ge₂Sb₂Te₅ structures is depicted in Figure 5, which represents the bond degree with respect to the |V|/G ratio [37]. The nearly-closed-shell S1-b1 and S2-b1 interactions are clearly well separated from the other bonding interactions, bearing both small bond degrees and |V|/G ratios, which is a sign of rather weak interactions. Interestingly, all the bonding interactions in these compounds share the same rigidity, as defined by Yang et al. [37,38], with G/ρ ≈ 0.41 Ha/e, which corresponds to the slope of the line (Figure 5). The weak interlayer interaction, characterized by a low bond degree and low electron density having a tendency to escape from this region, suggests a strong resistance to electron conduction along the c-axis of the compounds; hence, a low electrical conductivity σ<sub>zz</sub> and a high Seebeck coefficient are expected in this direction.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Electron density gradient of S1-Ge₂Sb₂Te₅ (a) and S2-Ge₂Sb₂Te₅ (b) plotted in the (110) plane. b: bond critical point; r: ring critical point; c: cage critical point.</figcaption></figure>

**Table 5.** Bond critical point (BCP) properties of S1- and S2-Sb₂Te₂Ge₅ calculated from the topological analysis of the electron density according to the quantum theory of atoms in molecules. b1 to b5: bond critical point (see Figure 4); r₁ and r₂: distances (in pm) between the BCP and the closest atomic nuclei; θ: angle (in degrees) between the BCP and the two nuclei it connects (the BCP being the apex); ρ: electron density at the BCP (in electron/Å³); ∇²ρ: electron density Laplacian at the BCP (in electron/Å⁵); G, V and H: kinetic, potential and total energy density at the BCP (in Ha/Å³); H/ρ: bond degree (in Ha/electron).

<table><thead><tr><th>BCP</th><th>r<sub>1</sub></th><th>r<sub>2</sub></th><th>r<sub>1</sub>/r<sub>2</sub></th><th>θ</th><th>ρ</th><th>∇<sup>2</sup>ρ</th><th>G</th><th>V</th><th>H</th><th>|V|/G</th><th>H/ρ</th></tr></thead><tbody><tr><td>S1-b1</td><td>190.8</td><td>190.8</td><td>1.000</td><td>180.00</td><td>1.145</td><td>1.779</td><td>0.46</td><td>-0.48</td><td>-0.02</td><td>1.04</td><td>-0.02</td></tr><tr><td>S1-b2</td><td>163.7</td><td>153.9</td><td>1.063</td><td>179.76</td><td>3.301</td><td>2.343</td><td>1.37</td><td>-2.15</td><td>-0.78</td><td>1.57</td><td>-0.24</td></tr></tbody></table>
The bonding structures of S1- and S2-Ge₂Sb₂Te₅ can be characterized by the electron density gradient (Figure 4). The zero-gradient lines delineating the so-called atomic basins allow for visualizing the atomic volumes, which are noticeably larger for Te and Ge than for Sb and are far from spherical in shape. The line of maximum gradient joining two atoms and passing through a bond critical point characterizes a bond path. In this respect, the Te atoms next to and on either side of the interlayer space are connected by a bond path (BCP, labeled b1 in Figure 4). The interlayer interaction is characterized by a low density, a slightly negative total energy density and a negative, though small, bond degree (H/ρ) (Table 5). The positive value of the electron density Laplacian (∇²ρ) underscores the tendency of the electron density to escape from the BCP region. Strictly speaking, the interlayer bonding interaction corresponds to transit region bonding, as defined by Espinosa [36], with the characteristics 1<|V|/G<2, H/ρ<0 and ∇²ρ>0, although this interaction is very close to that corresponding to the closed-shell region, which includes the van der Waals interactions. All the other bond interactions, labeled b2 to b5, are characteristic of the transit region bonding interaction (Table 5). The variety of bonds in S1- and S2-Ge₂Sb₂Te₅ structures is depicted in Figure 5, which represents the bond degree with respect to the |V|/G ratio [37]. The nearly-closed-shell S1-b1 and S2-b1 interactions are clearly well separated from the other bonding interactions, bearing both small bond degrees and |V|/G ratios, which is a sign of rather weak interactions. Interestingly, all the bonding interactions in these compounds share the same rigidity, as defined by Yang et al. [37,38], with G/ρ ≈ 0.41 Ha/e, which corresponds to the slope of the line (Figure 5). The weak interlayer interaction, characterized by a low bond degree and low electron density having a tendency to escape from this region, suggests a strong resistance to electron conduction along the c-axis of the compounds; hence, a low electrical conductivity σ<sub>zz</sub> and a high Seebeck coefficient are expected in this direction.

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Electron density gradient of S1-Ge₂Sb₂Te₅ (a) and S2-Ge₂Sb₂Te₅ (b) plotted in the (110) plane. b: bond critical point; r: ring critical point; c: cage critical point.</figcaption></figure>

**Table 5.** Bond critical point (BCP) properties of S1- and S2-Sb₂Te₂Ge₅ calculated from the topological analysis of the electron density according to the quantum theory of atoms in molecules. b1 to b5: bond critical point (see Figure 4); r₁ and r₂: distances (in pm) between the BCP and the closest atomic nuclei; θ: angle (in degrees) between the BCP and the two nuclei it connects (the BCP being the apex); ρ: electron density at the BCP (in electron/Å³); ∇²ρ: electron density Laplacian at the BCP (in electron/Å⁵); G, V and H: kinetic, potential and total energy density at the BCP (in Ha/Å³); H/ρ: bond degree (in Ha/electron).

<table><thead><tr><th>BCP</th><th>r<sub>1</sub></th><th>r<sub>2</sub></th><th>r<sub>1</sub>/r<sub>2</sub></th><th>θ</th><th>ρ</th><th>∇<sup>2</sup>ρ</th><th>G</th><th>V</th><th>H</th><th>|V|/G</th><th>H/ρ</th></tr></thead><tbody><tr><td>S1-b1</td><td>190.8</td><td>190.8</td><td>1.000</td><td>180.00</td><td>1.145</td><td>1.779</td><td>0.46</td><td>-0.48</td><td>-0.02</td><td>1.04</td><td>-0.02</td></tr><tr><td>S1-b2</td><td>163.7</td><td>153.9</td><td>1.063</td><td>179.76</td><td>3.301</td><td>2.343</td><td>1.37</td><td>-2.15</td><td>-0.78</td><td>1.57</td><td>-0.24</td></tr></tbody></table>
<table><caption>**Table 5. Cont.**</caption><thead><tr><th>BCP</th><th>r<sub>1</sub></th><th>r<sub>2</sub></th><th>r<sub>1</sub>/r<sub>2</sub></th><th>θ</th><th>ρ</th><th>∇<sup>2</sup>ρ</th><th>G</th><th>V</th><th>H</th><th>|V|/G</th><th>H/ρ</th></tr></thead><tbody><tr><td>S1-b3</td><td>136.8</td><td>159.8</td><td>0.856</td><td>179.56</td><td>3.791</td><td>1.952</td><td>1.55</td><td>-2.62</td><td>-1.07</td><td>1.69</td><td>-0.28</td></tr><tr><td>S1-b4</td><td>155.5</td><td>144.0</td><td>1.080</td><td>179.41</td><td>4.608</td><td>1.399</td><td>1.93</td><td>-3.52</td><td>-1.58</td><td>1.82</td><td>-0.34</td></tr><tr><td>S1-b5</td><td>138.1</td><td>160.8</td><td>0.859</td><td>179.97</td><td>3.645</td><td>2.034</td><td>1.49</td><td>-2.47</td><td>-0.98</td><td>1.66</td><td>-0.27</td></tr><tr><td>S2-b1</td><td>192.6</td><td>192.6</td><td>1.000</td><td>180.00</td><td>1.089</td><td>1.717</td><td>0.44</td><td>-0.45</td><td>-0.01</td><td>1.02</td><td>-0.01</td></tr><tr><td>S2-b2</td><td>144.5</td><td>155.6</td><td>0.929</td><td>179.55</td><td>4.530</td><td>1.558</td><td>1.91</td><td>-3.44</td><td>-1.52</td><td>1.80</td><td>-0.34</td></tr><tr><td>S2-b3</td><td>170.7</td><td>152.9</td><td>1.116</td><td>177.53</td><td>2.368</td><td>2.370</td><td>0.96</td><td>-1.32</td><td>-0.36</td><td>1.38</td><td>-0.15</td></tr><tr><td>S2-b4</td><td>128.7</td><td>153.2</td><td>0.840</td><td>178.99</td><td>5.012</td><td>0.487</td><td>2.04</td><td>-3.95</td><td>-1.92</td><td>1.94</td><td>-0.38</td></tr><tr><td>S2-b5</td><td>153.3</td><td>163.0</td><td>0.941</td><td>179.96</td><td>3.384</td><td>2.299</td><td>1.40</td><td>-2.22</td><td>-0.82</td><td>1.59</td><td>-0.24</td></tr></tbody></table>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Bond degree vs. |V|/G ratio of the bonding interactions in S1- and S2-Sb₂Te₂Ge₅.</figcaption></figure>

The average Seebeck coefficient (Figure 6), as calculated with the WC-SOC and BTE approach, shows a broad peak for both *n*- and *p*-doped S1- and S2-Ge₂Sb₂Te₅, the maximum of which is being shifted towards a higher doping level as the temperature increases. The magnitude of the maxima for the S1 stacking is larger (in absolute value) than that of the S2 stacking, which is a consequence of the larger gap of S1 compared with that of S2. Overall, a high maximum Seebeck coefficient is found that ranges between 120 µV K⁻¹ and 300 µV K⁻¹, in absolute values. As expected from the bonding analysis, the cross-plane Seebeck coefficient is noticeably higher than that in the *ab* plane (see Figures S1 and S2 in the supplemental data), especially for stacking 1 at temperatures above 300 K. This confirms that the gap between the layers hinders electron conduction, hence improving S in this direction.

As the linearized Boltzmann transport equation leads to τ-scaled electrical conductivity σ/τ, where τ is the carrier relaxation time, the values of τ for both electrons and holes were determined using the deformation potential theory [23,24]. According to this theory, the relaxation time can be determined from the following relation:

$$ \tau_i = \frac{2\sqrt{2\pi}\hbar^4}{3(k_B T)^{3/2}} \times \frac{C_{ii}}{m_i^{*3/2} E_d^{i2}} \quad (1) $$

where $C_{ii}$ is the elastic constant in the *i* direction, $m_i^*$ is the carrier effective mass and $E_d^i$ is the valence (for holes) or conduction (for electrons) band orbital energy variation when the lattice is subjected to a distortion along the direction *i*:

$$ E_d^i = \frac{dE}{d(l^i/l_0^i)} \quad (2) $$

$l_0$ being the undistorted lattice parameter. The values of τ calculated at 300 K, 500 K, 700 K and 900 K are gathered in Table 6. Overall, the relaxation times of the electrons lie in the
<table><caption>**Table 5.** *Cont.*</caption><thead><tr><th>BCP</th><th>r<sub>1</sub></th><th>r<sub>2</sub></th><th>r<sub>1</sub>/r<sub>2</sub></th><th>θ</th><th>ρ</th><th>∇<sup>2</sup>ρ</th><th>G</th><th>V</th><th>H</th><th>|V|/G</th><th>H/ρ</th></tr></thead><tbody><tr><td>S1-b3</td><td>136.8</td><td>159.8</td><td>0.856</td><td>179.56</td><td>3.791</td><td>1.952</td><td>1.55</td><td>-2.62</td><td>-1.07</td><td>1.69</td><td>-0.28</td></tr><tr><td>S1-b4</td><td>155.5</td><td>144.0</td><td>1.080</td><td>179.41</td><td>4.608</td><td>1.399</td><td>1.93</td><td>-3.52</td><td>-1.58</td><td>1.82</td><td>-0.34</td></tr><tr><td>S1-b5</td><td>138.1</td><td>160.8</td><td>0.859</td><td>179.97</td><td>3.645</td><td>2.034</td><td>1.49</td><td>-2.47</td><td>-0.98</td><td>1.66</td><td>-0.27</td></tr><tr><td>S2-b1</td><td>192.6</td><td>192.6</td><td>1.000</td><td>180.00</td><td>1.089</td><td>1.717</td><td>0.44</td><td>-0.45</td><td>-0.01</td><td>1.02</td><td>-0.01</td></tr><tr><td>S2-b2</td><td>144.5</td><td>155.6</td><td>0.929</td><td>179.55</td><td>4.530</td><td>1.558</td><td>1.91</td><td>-3.44</td><td>-1.52</td><td>1.80</td><td>-0.34</td></tr><tr><td>S2-b3</td><td>170.7</td><td>152.9</td><td>1.116</td><td>177.53</td><td>2.368</td><td>2.370</td><td>0.96</td><td>-1.32</td><td>-0.36</td><td>1.38</td><td>-0.15</td></tr><tr><td>S2-b4</td><td>128.7</td><td>153.2</td><td>0.840</td><td>178.99</td><td>5.012</td><td>0.487</td><td>2.04</td><td>-3.95</td><td>-1.92</td><td>1.94</td><td>-0.38</td></tr><tr><td>S2-b5</td><td>153.3</td><td>163.0</td><td>0.941</td><td>179.96</td><td>3.384</td><td>2.299</td><td>1.40</td><td>-2.22</td><td>-0.82</td><td>1.59</td><td>-0.24</td></tr></tbody></table>

<figure><img src="image_10.png" /><figcaption>**Figure 5.** Bond degree vs. |V|/G ratio of the bonding interactions in S1- and S2-Sb<sub>2</sub>Te<sub>2</sub>Ge<sub>5</sub>.</figcaption></figure>

The average Seebeck coefficient (Figure 6), as calculated with the WC-SOC and BTE approach, shows a broad peak for both *n*- and *p*-doped S1- and S2-Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub>, the maximum of which is being shifted towards a higher doping level as the temperature increases. The magnitude of the maxima for the S1 stacking is larger (in absolute value) than that of the S2 stacking, which is a consequence of the larger gap of S1 compared with that of S2. Overall, a high maximum Seebeck coefficient is found that ranges between 120 µV K<sup>-1</sup> and 300 µV K<sup>-1</sup>, in absolute values. As expected from the bonding analysis, the cross-plane Seebeck coefficient is noticeably higher than that in the *ab* plane (see Figures S1 and S2 in the supplemental data), especially for stacking 1 at temperatures above 300 K. This confirms that the gap between the layers hinders electron conduction, hence improving S in this direction.

As the linearized Boltzmann transport equation leads to τ-scaled electrical conductivity σ/τ, where τ is the carrier relaxation time, the values of τ for both electrons and holes were determined using the deformation potential theory [23,24]. According to this theory, the relaxation time can be determined from the following relation:

$$ \tau_i = \frac{2\sqrt{2\pi}\hbar^4}{3(k_B T)^{3/2}} \times \frac{C_{ii}}{m_i^{*3/2} E_d^{i2}} \quad (1) $$

where $C_{ii}$ is the elastic constant in the *i* direction, $m_i^*$ is the carrier effective mass and $E_d^i$ is the valence (for holes) or conduction (for electrons) band orbital energy variation when the lattice is subjected to a distortion along the direction *i*:

$$ E_d^i = \frac{dE}{d(l^i/l_0^i)} \quad (2) $$

$l_0$ being the undistorted lattice parameter. The values of τ calculated at 300 K, 500 K, 700 K and 900 K are gathered in Table 6. Overall, the relaxation times of the electrons lie in the
range $10^{-13}$–$10^{-11}$ s. The scattering time in the z direction is shorter than that in the *ab* plane, irrespective of the stacking type and temperature, highlighting the increased scattering by the interlayer interface. It is hence expected that there is a larger in-plane electrical conductivity and a smaller cross-plane one.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Seebeck coefficient at 300 K, 500 K, 700 K and 900 K of Ge₂Sb₂Te₅ with (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level.</figcaption></figure>

<table><caption>**Table 6.** Relaxation time (in s) of the electron (e) and hole (h) carriers in the *ab* plane and *z* direction for stackings S1 and S2 at 300 K, 500 K, 700 K and 900 K.</caption><thead><tr><th>T (K)</th><th>Stacking</th><th>τ<sub>e</sub><sup>ab</sup></th><th>τ<sub>e</sub><sup>z</sup></th><th>τ<sub>h</sub><sup>ab</sup></th><th>τ<sub>h</sub><sup>z</sup></th></tr></thead><tbody><tr><td rowspan="2">300</td><td>S1</td><td>6.70 × 10<sup>-12</sup></td><td>1.25 × 10<sup>-12</sup></td><td>7.49 × 10<sup>-12</sup></td><td>1.92 × 10<sup>-12</sup></td></tr><tr><td>S2</td><td>2.08 × 10<sup>-12</sup></td><td>1.98 × 10<sup>-12</sup></td><td>3.00 × 10<sup>-11</sup></td><td>1.66 × 10<sup>-11</sup></td></tr><tr><td rowspan="2">500</td><td>S1</td><td>3.11 × 10<sup>-12</sup></td><td>5.81 × 10<sup>-13</sup></td><td>3.48 × 10<sup>-12</sup></td><td>8.92 × 10<sup>-13</sup></td></tr><tr><td>S2</td><td>9.67 × 10<sup>-13</sup></td><td>9.20 × 10<sup>-13</sup></td><td>1.39 × 10<sup>-11</sup></td><td>7.71 × 10<sup>-12</sup></td></tr><tr><td rowspan="2">700</td><td>S1</td><td>1.88 × 10<sup>-12</sup></td><td>3.51 × 10<sup>-13</sup></td><td>2.10 × 10<sup>-12</sup></td><td>5.39 × 10<sup>-13</sup></td></tr><tr><td>S2</td><td>5.84 × 10<sup>-13</sup></td><td>5.56 × 10<sup>-13</sup></td><td>8.42 × 10<sup>-12</sup></td><td>4.66 × 10<sup>-12</sup></td></tr><tr><td rowspan="2">900</td><td>S1</td><td>1.29 × 10<sup>-12</sup></td><td>2.41 × 10<sup>-13</sup></td><td>1.44 × 10<sup>-12</sup></td><td>3.70 × 10<sup>-13</sup></td></tr><tr><td>S2</td><td>4.00 × 10<sup>-13</sup></td><td>3.81 × 10<sup>-13</sup></td><td>5.77 × 10<sup>-12</sup></td><td>3.19 × 10<sup>-12</sup></td></tr></tbody></table>

The average electrical conductivity (Figure 7) shows similar trends for both S1 and S2, irrespective of the type of doping; a plateau is observed up to $10^{18}$–$10^{19}$ carriers/cm⁻³ and then a steady increase when the compound becomes a degenerate semiconductor. More interesting are the in-plane and cross-plane features of the electrical conductivity, in line with what we expected from the bonding analysis (Figures S3 and S4). Indeed, the cross-plane electrical conductivity is substantially lower by about one order of magnitude for stacking 1, irrespective of the *n*- or *p*-doping type, confirming the barrier-like behavior of the gap between the layers. It is obvious that the gap behavior has much more impact on the electrical conductivity than on the Seebeck coefficient. Assuming the Mott formula for S holds ($S \propto d \ln(\sigma)$), this behavior might be explained by the fact that the large effect on σ is damped by the log function in S.
range $10^{-13}$–$10^{-11}$ s. The scattering time in the z direction is shorter than that in the *ab* plane, irrespective of the stacking type and temperature, highlighting the increased scattering by the interlayer interface. It is hence expected that there is a larger in-plane electrical conductivity and a smaller cross-plane one.

<figure><img src="image_11.png" /><figcaption>**Figure 6.** Seebeck coefficient at 300 K, 500 K, 700 K and 900 K of Ge₂Sb₂Te₅ with (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level.</figcaption></figure>

<table><caption>**Table 6.** Relaxation time (in s) of the electron (e) and hole (h) carriers in the *ab* plane and *z* direction for stackings S1 and S2 at 300 K, 500 K, 700 K and 900 K.</caption><thead><tr><th>T (K)</th><th>Stacking</th><th>τ<sub>e</sub><sup>ab</sup></th><th>τ<sub>e</sub><sup>z</sup></th><th>τ<sub>h</sub><sup>ab</sup></th><th>τ<sub>h</sub><sup>z</sup></th></tr></thead><tbody><tr><td rowspan="2">300</td><td>S1</td><td>6.70 × 10<sup>-12</sup></td><td>1.25 × 10<sup>-12</sup></td><td>7.49 × 10<sup>-12</sup></td><td>1.92 × 10<sup>-12</sup></td></tr><tr><td>S2</td><td>2.08 × 10<sup>-12</sup></td><td>1.98 × 10<sup>-12</sup></td><td>3.00 × 10<sup>-11</sup></td><td>1.66 × 10<sup>-11</sup></td></tr><tr><td rowspan="2">500</td><td>S1</td><td>3.11 × 10<sup>-12</sup></td><td>5.81 × 10<sup>-13</sup></td><td>3.48 × 10<sup>-12</sup></td><td>8.92 × 10<sup>-13</sup></td></tr><tr><td>S2</td><td>9.67 × 10<sup>-13</sup></td><td>9.20 × 10<sup>-13</sup></td><td>1.39 × 10<sup>-11</sup></td><td>7.71 × 10<sup>-12</sup></td></tr><tr><td rowspan="2">700</td><td>S1</td><td>1.88 × 10<sup>-12</sup></td><td>3.51 × 10<sup>-13</sup></td><td>2.10 × 10<sup>-12</sup></td><td>5.39 × 10<sup>-13</sup></td></tr><tr><td>S2</td><td>5.84 × 10<sup>-13</sup></td><td>5.56 × 10<sup>-13</sup></td><td>8.42 × 10<sup>-12</sup></td><td>4.66 × 10<sup>-12</sup></td></tr><tr><td rowspan="2">900</td><td>S1</td><td>1.29 × 10<sup>-12</sup></td><td>2.41 × 10<sup>-13</sup></td><td>1.44 × 10<sup>-12</sup></td><td>3.70 × 10<sup>-13</sup></td></tr><tr><td>S2</td><td>4.00 × 10<sup>-13</sup></td><td>3.81 × 10<sup>-13</sup></td><td>5.77 × 10<sup>-12</sup></td><td>3.19 × 10<sup>-12</sup></td></tr></tbody></table>

The average electrical conductivity (Figure 7) shows similar trends for both S1 and S2, irrespective of the type of doping; a plateau is observed up to $10^{18}$–$10^{19}$ carriers/cm⁻³ and then a steady increase when the compound becomes a degenerate semiconductor. More interesting are the in-plane and cross-plane features of the electrical conductivity, in line with what we expected from the bonding analysis (Figures S3 and S4). Indeed, the cross-plane electrical conductivity is substantially lower by about one order of magnitude for stacking 1, irrespective of the *n*- or *p*-doping type, confirming the barrier-like behavior of the gap between the layers. It is obvious that the gap behavior has much more impact on the electrical conductivity than on the Seebeck coefficient. Assuming the Mott formula for S holds ($S \propto d \ln(\sigma)$), this behavior might be explained by the fact that the large effect on σ is damped by the log function in S.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Electrical conductivity at 300 K, 500 K, 700 K and 900 K of Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> with (**a**,**b**) stacking 1 and (**c**,**d**) stacking 2 with respect to (**a**,**c**) *n*-type doping level and (**b**,**d**) *p*-type doping level.</figcaption></figure>

Ultimately, the power factor characterizes the electrical transport properties of a thermoelectric material. The power factor (PF) reads $S^2\sigma$. For S1, the PF exhibits a wide peak at large carriers' concentrations of about $10^{21} \text{ cm}^{-3}$ for *n*-type doping (Figure 8a) and $10^{20} \text{ cm}^{-3}$ for *p*-type doping (Figure 8b). As it turns out, the overall power factor of S1 Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> is mostly contributed by the in-plane PF, both for *n*-type (Figure S5a) and *p*-type doping (Figure S5b). The cross-plane PF (Figures S6a and S6b) is much weaker than the in-plane one, which can be related to the much smaller electrical conductivity that overkills the otherwise higher Seebeck coefficient.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Power factor ($S^2\sigma$) at 300 K, 500 K, 700 K and 900 K of Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> with (**a**,**b**) stacking 1 and (**c**,**d**) stacking 2 with respect to (**a**,**c**) *n*-type doping level and (**b**,**d**) *p*-type doping level.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 7.** Electrical conductivity at 300 K, 500 K, 700 K and 900 K of Ge₂Sb₂Te₅ with (**a**,**b**) stacking 1 and (**c**,**d**) stacking 2 with respect to (**a**,**c**) *n*-type doping level and (**b**,**d**) *p*-type doping level.</figcaption></figure>

Ultimately, the power factor characterizes the electrical transport properties of a thermoelectric material. The power factor (PF) reads $S^2\sigma$. For S1, the PF exhibits a wide peak at large carriers' concentrations of about $10^{21} \text{ cm}^{-3}$ for *n*-type doping (Figure 8a) and $10^{20} \text{ cm}^{-3}$ for *p*-type doping (Figure 8b). As it turns out, the overall power factor of S1 Ge₂Sb₂Te₅ is mostly contributed by the in-plane PF, both for *n*-type (Figure S5a) and *p*-type doping (Figure S5b). The cross-plane PF (Figures S6a and S6b) is much weaker than the in-plane one, which can be related to the much smaller electrical conductivity that overkills the otherwise higher Seebeck coefficient.

<figure><img src="image_13.png" /><figcaption>**Figure 8.** Power factor ($S^2\sigma$) at 300 K, 500 K, 700 K and 900 K of Ge₂Sb₂Te₅ with (**a**,**b**) stacking 1 and (**c**,**d**) stacking 2 with respect to (**a**,**c**) *n*-type doping level and (**b**,**d**) *p*-type doping level.</figcaption></figure>
(b) Thermal Transport

Using the DFPT approach, the anharmonic force constants were calculated for the two stacking sequences in the *ab* plane and along the *c*-axis, and the corresponding phonon spectrum and DOS are presented in Figure 9a,b. We observed a small frequency gap between 3.9 and 4.0 THz for stacking 2, which is comparable to but smaller in extent than the gap observed in the Pb₂Bi₂Te₅ compound [39]. The DOS shows the prominent contribution of the Te atoms over the whole spectrum range, especially in the acoustic and low-energy (1.5–3 THz) optical phonons. The Sb contribution is appreciable except between 2.2 and 3 THz, and the Ge contribution is appreciable above 2.2 THz only. The lattice thermal conductivity $κ_l$ was evaluated using the single-mode relaxation time approximation (RTA):

$$κ_l = \frac{1}{N V_0} \sum_{\lambda} C_{\lambda} v_{\lambda} \otimes v_{\lambda} \tau_{\lambda} \quad (3)$$

where *N* is the number of *q*-points, $V_0$ is the volume of the unit cell and $C_λ$, $v_λ$ and $τ_λ$ are heat capacity, group velocity and phonon lifetime, respectively.

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Phonon spectrum and total and projected phonon density of states of (a) stacking 1 and (b) stacking 2 at equilibrium. (c) Calculated lattice thermal conductivity of stackings 1 and 2. The evolution is plotted for the in-layer *a*-axis, cross-layer *c*-axis and average values.</figcaption></figure>

The calculated lattice thermal conductivity for the two stacking sequences in the *ab* plane ($κ_{xx}$) and along the *c*-axis $κ_{zz}$ and the average thermal conductivity $κ_{l,ave} = (2 × κ_{xx} + κ_{zz})/3$ are shown in Figure 9c. Overall, the lattice thermal conductivity is lower for S1 than for S2, which is caused by the smaller in-plane contribution of S1 compared with that of S2; the cross-plane contribution is the same for both sequences. Therefore, in S2, one can state that the frequency gap shows only marginal effects on the total lattice thermal conductivity; a very small plateau is observed in $κ_{xx}$ between 3.9 and 4.0 THz, which is slightly reflected in $κ_{l,ave}$ (Figure 10a). Interestingly, in stacking 2 the cross-plane lattice thermal conductivity exhibits a wide plateau at 0.7 W/(m·K) from 2.5 THz onward, but this behavior is not strong enough to counterbalance the steep increase in the in-plane contribution at low energy (acoustic phonons) that, by far, surpasses the S1 $κ_{xx}$ and $κ_{zz}$ contributions. The wide plateau of $κ_{zz}$ for S2 could be related to the small velocities of the optical phonons, especially above 3 THz (Figure 10d). Clearly, the overall smaller lattice thermal conductivity of S1 is caused by the smaller acoustic phonon scattering times (Figure 10b). Indeed, although we observe that the velocities of the acoustic phonons of S1 are about the same as those of S2 and that the velocities of the optical phonons of S1
(b) Thermal Transport

Using the DFPT approach, the anharmonic force constants were calculated for the two stacking sequences in the *ab* plane and along the *c*-axis, and the corresponding phonon spectrum and DOS are presented in Figure 9a,b. We observed a small frequency gap between 3.9 and 4.0 THz for stacking 2, which is comparable to but smaller in extent than the gap observed in the Pb₂Bi₂Te₅ compound [39]. The DOS shows the prominent contribution of the Te atoms over the whole spectrum range, especially in the acoustic and low-energy (1.5–3 THz) optical phonons. The Sb contribution is appreciable except between 2.2 and 3 THz, and the Ge contribution is appreciable above 2.2 THz only. The lattice thermal conductivity κ₁ was evaluated using the single-mode relaxation time approximation (RTA):

$$ \kappa_l = \frac{1}{N V_0} \sum_{\lambda} C_{\lambda} \nu_{\lambda} \otimes \nu_{\lambda} \tau_{\lambda} \quad (3) $$

where N is the number of *q*-points, V₀ is the volume of the unit cell and C_λ, ν_λ and τ_λ are heat capacity, group velocity and phonon lifetime, respectively.

<figure><img src="image_14.png" /><figcaption>**Figure 9.** Phonon spectrum and total and projected phonon density of states of (a) stacking 1 and (b) stacking 2 at equilibrium. (c) Calculated lattice thermal conductivity of stackings 1 and 2. The evolution is plotted for the in-layer *a*-axis, cross-layer *c*-axis and average values.</figcaption></figure>

The calculated lattice thermal conductivity for the two stacking sequences in the *ab* plane (κₓₓ) and along the *c*-axis κzz and the average thermal conductivity κl,ave = (2 × κₓₓ + κzz)/3 are shown in Figure 9c. Overall, the lattice thermal conductivity is lower for S1 than for S2, which is caused by the smaller in-plane contribution of S1 compared with that of S2; the cross-plane contribution is the same for both sequences. Therefore, in S2, one can state that the frequency gap shows only marginal effects on the total lattice thermal conductivity; a very small plateau is observed in κₓₓ between 3.9 and 4.0 THz, which is slightly reflected in κl,ave (Figure 10a). Interestingly, in stacking 2 the cross-plane lattice thermal conductivity exhibits a wide plateau at 0.7 W/(m·K) from 2.5 THz onward, but this behavior is not strong enough to counterbalance the steep increase in the in-plane contribution at low energy (acoustic phonons) that, by far, surpasses the S1 κₓₓ and κzz contributions. The wide plateau of κzz for S2 could be related to the small velocities of the optical phonons, especially above 3 THz (Figure 10d). Clearly, the overall smaller lattice thermal conductivity of S1 is caused by the smaller acoustic phonon scattering times (Figure 10b). Indeed, although we observe that the velocities of the acoustic phonons of S1 are about the same as those of S2 and that the velocities of the optical phonons of S1
are much higher than those of stacking 2, the phonon lifetimes of S2 are much longer than those of S1 in the low-frequency range of 0–3 THz (acoustic and optical modes), which contributes to 70% and 91% of the total $\kappa_l$ of S1 and S2, respectively. In a sense, the phonons of S2 transport heat more slowly but more surely, hence more efficiently.

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Calculated cumulative lattice thermal conductivity for Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> stacking 1 (blue) and stacking 2 (orange) with RTA (a). Phonon lifetimes for S1 and S2 (b). In-layer group velocity (c) and cross-layer group velocity (d) for both S1 and S2. Results obtained at 500 K.</figcaption></figure>

It has been stated in the literature [40,41] that low lattice thermal conductivity could be related to a strong coupling between low-energy optical phonons and acoustic phonons, leading to avoided crossing in the phonon band structure. As a tentative explanation of the lower thermal conductivity of S1 compared with that of S2, one may invoke such a larger coupling in S1 than in S2. Several avoided band crossings between acoustic and optical phonons may possibly appear (Figure 9a) in S1, e.g., at the A k-point around 1 THz, and on the K-Γ line (at $\xi \approx 1/3$ along this line around 2 THz and $\xi \approx 3/4$ around 1.2 THz), the former coupling being seemingly much weaker and the latter one being absent in S2 (Figure 9b). The strong phonon coupling has been attributed to the distortion of the structure around the central atom, namely the deviation, for instance, from an octahedral or a tetrahedral symmetry: the larger the deviation, the larger the mode anharmonicity. However, this effect cannot be evoked for Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> as the larger deviation of the central atom (Sb in our case) from the center of the octahedral cage (see Table S1) is observed for the S2 structure for which the lattice thermal conductivity is higher than that of S1. Hence, the metallic state of S2 (in the limit of the WC functional used here) is likely the reason for the higher lattice thermal conductivity of S2.

### (c) Figure of Merit

The figure of merit ZT can be evaluated from the power factor and the lattice thermal conductivities. The results of ZT as a function of carrier concentrations and temperatures for stackings 1 and 2 are shown in Figure 11. The ranges of conditions to obtain the maximum ZT for stacking 1 are 200–700 K/10<sup>18</sup>–10<sup>20</sup> holes/cm<sup>3</sup> for *p*-type doping (Figure 11a), and 150–300 K/10<sup>17</sup>–10<sup>19</sup> electrons/cm<sup>3</sup> and 900–1000 K/10<sup>20</sup>–10<sup>21</sup> electrons/cm<sup>3</sup> for *n*-type doping (Figure 11b), hence the suitable region for *n*-type doping splits into two parts, namely one located at a high temperature and high doping level and the other one at a low temperature and low doping level. The favorable working temperature is rather wide for the *p*-type, but for high temperatures *n*-type doping is more suitable. Interestingly,
are much higher than those of stacking 2, the phonon lifetimes of S2 are much longer than those of S1 in the low-frequency range of 0–3 THz (acoustic and optical modes), which contributes to 70% and 91% of the total $\kappa_l$ of S1 and S2, respectively. In a sense, the phonons of S2 transport heat more slowly but more surely, hence more efficiently.

<figure><img src="image_15.png" /><figcaption>**Figure 10.** Calculated cumulative lattice thermal conductivity for Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> stacking 1 (blue) and stacking 2 (orange) with RTA (a). Phonon lifetimes for S1 and S2 (b). In-layer group velocity (c) and cross-layer group velocity (d) for both S1 and S2. Results obtained at 500 K.</figcaption></figure>

It has been stated in the literature [40,41] that low lattice thermal conductivity could be related to a strong coupling between low-energy optical phonons and acoustic phonons, leading to avoided crossing in the phonon band structure. As a tentative explanation of the lower thermal conductivity of S1 compared with that of S2, one may invoke such a larger coupling in S1 than in S2. Several avoided band crossings between acoustic and optical phonons may possibly appear (Figure 9a) in S1, e.g., at the A k-point around 1 THz, and on the K-Γ line (at $\xi \approx 1/3$ along this line around 2 THz and $\xi \approx 3/4$ around 1.2 THz), the former coupling being seemingly much weaker and the latter one being absent in S2 (Figure 9b). The strong phonon coupling has been attributed to the distortion of the structure around the central atom, namely the deviation, for instance, from an octahedral or a tetrahedral symmetry: the larger the deviation, the larger the mode anharmonicity. However, this effect cannot be evoked for Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> as the larger deviation of the central atom (Sb in our case) from the center of the octahedral cage (see Table S1) is observed for the S2 structure for which the lattice thermal conductivity is higher than that of S1. Hence, the metallic state of S2 (in the limit of the WC functional used here) is likely the reason for the higher lattice thermal conductivity of S2.

### (c) Figure of Merit

The figure of merit ZT can be evaluated from the power factor and the lattice thermal conductivities. The results of ZT as a function of carrier concentrations and temperatures for stackings 1 and 2 are shown in Figure 11. The ranges of conditions to obtain the maximum ZT for stacking 1 are 200–700 K/10<sup>18</sup>–10<sup>20</sup> holes/cm<sup>3</sup> for *p*-type doping (Figure 11a), and 150–300 K/10<sup>17</sup>–10<sup>19</sup> electrons/cm<sup>3</sup> and 900–1000 K/10<sup>20</sup>–10<sup>21</sup> electrons/cm<sup>3</sup> for *n*-type doping (Figure 11b), hence the suitable region for *n*-type doping splits into two parts, namely one located at a high temperature and high doping level and the other one at a low temperature and low doping level. The favorable working temperature is rather wide for the *p*-type, but for high temperatures *n*-type doping is more suitable. Interestingly,
S1-Ge₂Sb₂Te₅ could be used both as a *p* and *n* leg in a thermoelectric converter at high temperatures (about 800 K) where their ZT values are predicted to be similar (about 1.7). The highest value is 2.3 for stacking 1 with *p*-type doping. By contrast, the suitable conditions for stacking 2 to obtain high ZT are unfavorable for applications as they correspond to temperatures well below the room temperature, irrespective of the doping type.

<figure><img src="image_14.png" /><figcaption>**Figure 11.** Calculated figures of merit ZT of Ge₂Sb₂Te₅ stacking 1 (**a**,**b**) and stacking 2 (**c**,**d**). Both *p*-type (**a**,**c**) and *n*-type (**b**,**d**) doping are considered.</figcaption></figure>

## 3.2. Electronic Transport Properties of the Ge₂Sb₂Te₅ Compounds under Tensile and Compressive Strains

As the calculation of the lattice thermal conductivity is still a computationally heavy task, it was not possible to calculate $\kappa_l$ for S1 and S2 for all the applied strains. Therefore, we focus in this part on the electronic transport properties and more specifically on the Seebeck coefficient and the power factor.

Biaxial, compressive ($\eta < 0$) and tensile ($\eta > 0$) strains were applied to the S1- and S2-Ge₂Sb₂Te₅ compounds along the basal *ab* plane. As expected, the electronic structure of the compounds is affected by the strains, which is manifested by the change in the energy gap that either slightly widens or strongly shrinks up to the point where the compounds become metallic (Figure 12). A concomitant evolution of the Seebeck coefficient is observed as it reaches a maximum value when the gap is largest; this is particularly evident for the S2 stacking (Figure 12b) for which the maximum gap comes up at +1% tensile strain that correlates exactly with the maximum of S for both *n*- and *p*-doped compounds. For S1, the behavior is somewhat different. We note first that the gap is almost constant in the range of -1.75% to +2% strains. At these extremes, the Seebeck coefficient reaches a maximum of about 380 µV K⁻¹ under tensile strain for *p*-type doping and a maximum of about -350 µV K⁻¹ under compressive strain for *n*-type doping. Clearly, the Seebeck coefficient can be noticeably improved by applying biaxial strains, which could be achieved experimentally by depositing Ge₂Sb₂Te₅ on a suitably selected support.

Under tensile strains with electron carriers and compressive strains with hole carriers, the power factor (PF) of the S1-Ge₂Sb₂Te₅ compound is improved (Figures 13–16), irrespective of the direction (*ab* plane or *z* direction). We note here that the electron scattering times are those of the unstrained structure. Still, noticeable differences can be underlined. The improvement is small when S1 is doped with electrons, whereas it is impressively increased when it is doped with holes. This behavior holds true whether transport occurs in-plane or cross-plane; thus, strains can counteract the fact that the cross-plane gap between the
S1-Ge₂Sb₂Te₅ could be used both as a *p* and *n* leg in a thermoelectric converter at high temperatures (about 800 K) where their ZT values are predicted to be similar (about 1.7). The highest value is 2.3 for stacking 1 with *p*-type doping. By contrast, the suitable conditions for stacking 2 to obtain high ZT are unfavorable for applications as they correspond to temperatures well below the room temperature, irrespective of the doping type.

<figure><img src="image_16.png" /><figcaption>**Figure 11.** Calculated figures of merit ZT of Ge₂Sb₂Te₅ stacking 1 (**a**,**b**) and stacking 2 (**c**,**d**). Both *p*-type (**a**,**c**) and *n*-type (**b**,**d**) doping are considered.</figcaption></figure>

## 3.2. Electronic Transport Properties of the Ge₂Sb₂Te₅ Compounds under Tensile and Compressive Strains

As the calculation of the lattice thermal conductivity is still a computationally heavy task, it was not possible to calculate $\kappa_l$ for S1 and S2 for all the applied strains. Therefore, we focus in this part on the electronic transport properties and more specifically on the Seebeck coefficient and the power factor.

Biaxial, compressive ($\eta < 0$) and tensile ($\eta > 0$) strains were applied to the S1- and S2-Ge₂Sb₂Te₅ compounds along the basal *ab* plane. As expected, the electronic structure of the compounds is affected by the strains, which is manifested by the change in the energy gap that either slightly widens or strongly shrinks up to the point where the compounds become metallic (Figure 12). A concomitant evolution of the Seebeck coefficient is observed as it reaches a maximum value when the gap is largest; this is particularly evident for the S2 stacking (Figure 12b) for which the maximum gap comes up at +1% tensile strain that correlates exactly with the maximum of S for both *n*- and *p*-doped compounds. For S1, the behavior is somewhat different. We note first that the gap is almost constant in the range of -1.75% to +2% strains. At these extremes, the Seebeck coefficient reaches a maximum of about 380 µV K⁻¹ under tensile strain for *p*-type doping and a maximum of about -350 µV K⁻¹ under compressive strain for *n*-type doping. Clearly, the Seebeck coefficient can be noticeably improved by applying biaxial strains, which could be achieved experimentally by depositing Ge₂Sb₂Te₅ on a suitably selected support.

Under tensile strains with electron carriers and compressive strains with hole carriers, the power factor (PF) of the S1-Ge₂Sb₂Te₅ compound is improved (Figures 13–16), irrespective of the direction (*ab* plane or *z* direction). We note here that the electron scattering times are those of the unstrained structure. Still, noticeable differences can be underlined. The improvement is small when S1 is doped with electrons, whereas it is impressively increased when it is doped with holes. This behavior holds true whether transport occurs in-plane or cross-plane; thus, strains can counteract the fact that the cross-plane gap between the
layers hinders electronic transport. Hence, on average, under compressive strains and with hole doping, S1 exhibits an excellent power factor ranging from 2 to 4 W m⁻¹ K⁻² between 700 K and 300 K (Figure S8). By contrast, under tensile strains and with electron doping, S1 shows a similar PF to that when unstrained (Figure S9). Overall, the same behavior is observed for the S2-Ge₂Sb₂Te₅ compound when strained (Figures S10–S15).

<figure><img src="image_15.png" /><figcaption>**Figure 12.** Energy gap and Seebeck coefficient at 300 K vs. compressive (η < 0) and tensile (η > 0) biaxial strains applied to (a) S1 and (b) S2 Ge₂Sb₂Te₅ compounds.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 13.** *ab* plane power factor of S1-Ge₂Sb₂Te₅ calculated under compressive strains (η < 0) for electron carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K.</figcaption></figure>
layers hinders electronic transport. Hence, on average, under compressive strains and with hole doping, S1 exhibits an excellent power factor ranging from 2 to 4 W m⁻¹ K⁻² between 700 K and 300 K (Figure S8). By contrast, under tensile strains and with electron doping, S1 shows a similar PF to that when unstrained (Figure S9). Overall, the same behavior is observed for the S2-Ge₂Sb₂Te₅ compound when strained (Figures S10–S15).

<figure><img src="image_17.png" /><figcaption>**Figure 12.** Energy gap and Seebeck coefficient at 300 K vs. compressive (η < 0) and tensile (η > 0) biaxial strains applied to (a) S1 and (b) S2 Ge₂Sb₂Te₅ compounds.</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Figure 13.** *ab* plane power factor of S1-Ge₂Sb₂Te₅ calculated under compressive strains (η < 0) for electron carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Figure 14.** ab plane power factor of S1-Ge₂Sb₂Te₅ calculated under tensile strains (η > 0) for hole carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K.</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Figure 15.** z direction power factor of S1-Ge₂Sb₂Te₅ calculated under compressive strains (η < 0) for electron carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Figure 14.** ab plane power factor of S1-Ge₂Sb₂Te₅ calculated under tensile strains (η > 0) for hole carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K.</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>**Figure 15.** z direction power factor of S1-Ge₂Sb₂Te₅ calculated under compressive strains (η < 0) for electron carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Figure 16.** z direction power factor of S1-Ge₂Sb₂Te₅ calculated under tensile strains (η > 0) for hole carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K.</figcaption></figure>

4. Conclusions

In this paper, a thorough investigation of the electronic and thermoelectric properties of the Ge₂Sb₂Te₅ compound is presented. The two different stacking sequences, namely stacking 1 (S1) and stacking 2 (S2), for which Ge and Sb can be swapped, were investigated by combining density functional theory, the quantum theory of atoms in molecules (QTAIM) and the Boltzmann transport theory. At the hybrid functional level including spin-orbit interaction, both S1 and S2 are shown to be semiconducting materials with close, small gap energies (0.090 eV and 0.086 eV, respectively). S1 bears a direct gap whereas S2 has an indirect one. From the QTAIM theory, the analysis of the electron density Laplacian and energy densities at the bond critical points allows us to conclude that the interlayer bonding interaction is of transit region nature, though sharing features of closed-shell interaction. The thermoelectric properties of both unstrained and strained S1 and S2 stackings were calculated. It was found that the S1 stacking is more suitable for use in a thermogenerator, both as a *p* and an *n* leg, as it exhibits high ZT values under certain ranges of temperatures and doping levels. Under strains, we observe that the thermoelectric properties of both S1 and S2 are improved. The improvement is more obvious when the material is *p*-doped under compressive strains than when it is *n*-doped under tensile strains. Remarkably, the strains are amenable to preventing the cross-plane electronic transport hindrance caused by the interlayer gap.

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ma16145015/s1: Figure S1: Seebeck coefficient in the ab plane of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S2: Seebeck coefficient in the cross-plane (c-axis direction) of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S3: Electrical conductivity in the ab plane of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S4:
<figure><img src="image_21.png" /><figcaption>**Figure 16.** z direction power factor of S1-Ge₂Sb₂Te₅ calculated under tensile strains (η > 0) for hole carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K.</figcaption></figure>

**4. Conclusions**

In this paper, a thorough investigation of the electronic and thermoelectric properties of the Ge₂Sb₂Te₅ compound is presented. The two different stacking sequences, namely stacking 1 (S1) and stacking 2 (S2), for which Ge and Sb can be swapped, were investigated by combining density functional theory, the quantum theory of atoms in molecules (QTAIM) and the Boltzmann transport theory. At the hybrid functional level including spin-orbit interaction, both S1 and S2 are shown to be semiconducting materials with close, small gap energies (0.090 eV and 0.086 eV, respectively). S1 bears a direct gap whereas S2 has an indirect one. From the QTAIM theory, the analysis of the electron density Laplacian and energy densities at the bond critical points allows us to conclude that the interlayer bonding interaction is of transit region nature, though sharing features of closed-shell interaction. The thermoelectric properties of both unstrained and strained S1 and S2 stackings were calculated. It was found that the S1 stacking is more suitable for use in a thermogenerator, both as a *p* and an *n* leg, as it exhibits high ZT values under certain ranges of temperatures and doping levels. Under strains, we observe that the thermoelectric properties of both S1 and S2 are improved. The improvement is more obvious when the material is *p*-doped under compressive strains than when it is *n*-doped under tensile strains. Remarkably, the strains are amenable to preventing the cross-plane electronic transport hindrance caused by the interlayer gap.

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ma16145015/s1: Figure S1: Seebeck coefficient in the ab plane of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S2: Seebeck coefficient in the cross-plane (c-axis direction) of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S3: Electrical conductivity in the ab plane of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S4:
Electrical conductivity in the cross-plane (c-axis direction) of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S5: Power factor (S²σ) in the ab plane of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S6: Power factor (S²σ) in the cross-plane (c-axis direction) of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S7: Deformed octahedral environment of (a) Ge in S1 stacking; (b) Sb in S2 stacking; Figure S8: Average power factor of *p*-doped S1-Ge₂Sb₂Te₅ with respect to tensile strains at the temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S9: Average power factor of *n*-doped S1-Ge₂Sb₂Te₅ with respect to compressive strains at the temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S10: *ab* plane power factor of S2-Ge₂Sb₂Te₅ calculated under compressive strains (η < 0) for electron carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S11: *ab* plane power factor of S2-Ge₂Sb₂Te₅ calculated under tensile strains (η > 0) for hole carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S12: *z* direction power factor of S2-Ge₂Sb₂Te₅ calculated under compressive strains (η < 0) for electron carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S13: *z* direction power factor of S2-Ge₂Sb₂Te₅ calculated under tensile strains (η > 0) for hole carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S14: Average power factor of *p*-doped S2-Ge₂Sb₂Te₅ with respect to tensile strains at the temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S15: Average power factor of *n*-doped S2-Ge₂Sb₂Te₅ with respect to compressive strains at the temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Table S1: Geometrical characteristics of the octahedral cage around the Ge atom in the S1 stacking and Sb atom in the S2 stacking.

**Author Contributions:** Conceptualization, P.B. and M.-C.R.; methodology, P.B. and M.-C.R.; software, J.T. and W.M.; formal analysis, J.T. and W.M.; data curation, J.T. and W.M.; writing—original draft preparation, J.T.; writing—review and editing, P.B. and M.-C.R.; supervision, P.B. and M.-C.R.; project administration, M.-C.R. All authors have read and agreed to the published version of the manuscript.

**Funding:** The PhD theses of J. Tian and W. Ma are financially supported by the China Scholarship Council (CSC).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data are available upon request to the authors.

**Acknowledgments:** This work was granted access to the HPC/AI resources of TGCC under the allocation A0130806881 made by GENCI. It was also granted access to the HPC resources of Aix-Marseille University.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Volta, A. Nueva memoria sull'elettricità animale. *Ann. Chem. Nat. Hist.* **1794**, *5*, 132–144.

2. Seebeck, T.J. Magnetic polarization of metals and ores by temperature differences. Abh. Königlichen Akad. Wiss. Berl. **1822**, 265–373.

3. Ioffe, I.A. *Semiconductor Thermoelements, and Thermoelectric Cooling*; Inforsearch: London, UK, 1957.

4. Furlong, R.R.; Wahlquist, E.J. US space missions using radioisotope power systems. *Nucl. News* **1999**, *42*, 26–35.

5. Orr, B.; Akbarzadeh, A.; Mochizuki, M.; Singh, R. A review of car waste heat recovery systems utilising thermoelectric generators and heat pipes. *Appl. Therm. Eng.* **2016**, *101*, 490–495. [CrossRef]

6. Rowe, D.M. (Ed.) *Thermoelectrics Handbook: Macro to Nano*; CRC Press Taylor & Francis: Boca Raton, FL, USA, 2006.

7. Rowe, D.M. (Ed.) *Thermoelectrics and Its Energy Harvesting: Materials, Preparation, and Characterization in Thermoelectrics*; CRC Press Taylor & Francis: Boca Raton, FL, USA, 2012.

8. Kolobov, A.V.; Fons, P.; Frenkel, A.I.; Ankudinov, A.L.; Tominaga, J.; Uruga, T. Understanding the phase-change mechanism of rewritable optical media. *Nat. Mater.* **2004**, *3*, 703–708. [CrossRef] [PubMed]

9. Miao, J.; Wang, P.; Zhou, P.; Huang, S.; Qian, D.; Yuan, Y.; Lan, R. Remarkable anisotropy in rhombohedral Ge₂Sb₂Te₅ compound: A promising thermoelectric material with multiple conduction bands and acoustic-optical branches coupling. *J. Alloys Compd.* **2022**, *900*, 163471. [CrossRef]

10. Hu, P.; Wei, T.-R.; Qi, P.; Cao, Y.; Yang, J.; Shi, X.; Chen, L. Largely Enhanced Seebeck Coefficient and Thermoelectric Performance by the Distortion of Electronic Density of States in Ge₂Sb₂Te₅. *ACS Appl. Mater. Interfaces* **2019**, *11*, 34046–34052. [CrossRef]
Electrical conductivity in the cross-plane (c-axis direction) of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S5: Power factor (S²σ) in the ab plane of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S6: Power factor (S²σ) in the cross-plane (c-axis direction) of Ge₂Sb₂Te₅ at 300 K, 500 K, 700 K and 900 K for (a,b) stacking 1 and (c,d) stacking 2 with respect to (a,c) *n*-type doping level and (b,d) *p*-type doping level; Figure S7: Deformed octahedral environment of (a) Ge in S1 stacking; (b) Sb in S2 stacking; Figure S8: Average power factor of *p*-doped S1-Ge₂Sb₂Te₅ with respect to tensile strains at the temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S9: Average power factor of *n*-doped S1-Ge₂Sb₂Te₅ with respect to compressive strains at the temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S10: *ab* plane power factor of S2-Ge₂Sb₂Te₅ calculated under compressive strains (η < 0) for electron carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S11: *ab* plane power factor of S2-Ge₂Sb₂Te₅ calculated under tensile strains (η > 0) for hole carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S12: *z* direction power factor of S2-Ge₂Sb₂Te₅ calculated under compressive strains (η < 0) for electron carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S13: *z* direction power factor of S2-Ge₂Sb₂Te₅ calculated under tensile strains (η > 0) for hole carriers and temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S14: Average power factor of *p*-doped S2-Ge₂Sb₂Te₅ with respect to tensile strains at the temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Figure S15: Average power factor of *n*-doped S2-Ge₂Sb₂Te₅ with respect to compressive strains at the temperatures of (a) 300 K, (b) 500 K, (c) 700 K and (d) 900 K; Table S1: Geometrical characteristics of the octahedral cage around the Ge atom in the S1 stacking and Sb atom in the S2 stacking.

**Author Contributions:** Conceptualization, P.B. and M.-C.R.; methodology, P.B. and M.-C.R.; software, J.T. and W.M.; formal analysis, J.T. and W.M.; data curation, J.T. and W.M.; writing—original draft preparation, J.T.; writing—review and editing, P.B. and M.-C.R.; supervision, P.B. and M.-C.R.; project administration, M.-C.R. All authors have read and agreed to the published version of the manuscript.

**Funding:** The PhD theses of J. Tian and W. Ma are financially supported by the China Scholarship Council (CSC).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data are available upon request to the authors.

**Acknowledgments:** This work was granted access to the HPC/AI resources of TGCC under the allocation A0130806881 made by GENCI. It was also granted access to the HPC resources of Aix-Marseille University.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Volta, A. Nueva memoria sull'elettricità animale. *Ann. Chem. Nat. Hist.* **1794**, *5*, 132-144.

2. Seebeck, T.J. Magnetic polarization of metals and ores by temperature differences. Abh. Königlichen Akad. Wiss. Berl. **1822**, 265-373.

3. Ioffe, I.A. *Semiconductor Thermoelements, and Thermoelectric Cooling*; Inforsearch: London, UK, 1957.

4. Furlong, R.R.; Wahlquist, E.J. US space missions using radioisotope power systems. *Nucl. News* **1999**, *42*, 26-35.

5. Orr, B.; Akbarzadeh, A.; Mochizuki, M.; Singh, R. A review of car waste heat recovery systems utilising thermoelectric generators and heat pipes. *Appl. Therm. Eng.* **2016**, *101*, 490-495. [CrossRef]

6. Rowe, D.M. (Ed.) *Thermoelectrics Handbook: Macro to Nano*; CRC Press Taylor & Francis: Boca Raton, FL, USA, 2006.

7. Rowe, D.M. (Ed.) *Thermoelectrics and Its Energy Harvesting: Materials, Preparation, and Characterization in Thermoelectrics*; CRC Press Taylor & Francis: Boca Raton, FL, USA, 2012.

8. Kolobov, A.V.; Fons, P.; Frenkel, A.I.; Ankudinov, A.L.; Tominaga, J.; Uruga, T. Understanding the phase-change mechanism of rewritable optical media. *Nat. Mater.* **2004**, *3*, 703-708. [CrossRef] [PubMed]

9. Miao, J.; Wang, P.; Zhou, P.; Huang, S.; Qian, D.; Yuan, Y.; Lan, R. Remarkable anisotropy in rhombohedral Ge₂Sb₂Te₅ compound: A promising thermoelectric material with multiple conduction bands and acoustic-optical branches coupling. *J. Alloys Compd.* **2022**, *900*, 163471. [CrossRef]

10. Hu, P.; Wei, T.-R.; Qi, P.; Cao, Y.; Yang, J.; Shi, X.; Chen, L. Largely Enhanced Seebeck Coefficient and Thermoelectric Performance by the Distortion of Electronic Density of States in Ge₂Sb₂Te₅. *ACS Appl. Mater. Interfaces* **2019**, *11*, 34046-34052. [CrossRef]
11. Du, W.; Gu, Y.; Wang, K.; Yang, X.; Xing, J.; Guo, K.; Luo, J.; Zhao, J.-T. Effective Mass Enhancement and Thermal Conductivity Reduction for Improving the Thermoelectric Properties of Pseudo-Binary Ge₂Sb₂Te₅. *Ann. Phys.* **2020**, *532*, 1900390. [CrossRef]

12. Hu, P.; Wei, T.-R.; Huang, S.-J.; Xia, X.-G.; Qiu, P.-F.; Yang, J.; Chen, L.-D.; Shi, X. Anion-site-modulated thermoelectric properties in Ge₂Sb₂Te₅-based compounds. *Rare Met.* **2020**, *39*, 1127–1133. [CrossRef]

13. Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. *Phys. Rev.* **1964**, *136*, B86–B871. [CrossRef]

14. Kohn, W.; Sham, L.J. Self-Consistent Equations Including Exchange and Correlation Effects. *Phys. Rev.* **1965**, *140*, A1133–A1138. [CrossRef]

15. Perdew, J.P.; Zunger, A. Self-interaction correction to density-functional approximations for many-electron systems. *Phys. Rev. B* **1981**, *23*, 5048–5079. [CrossRef]

16. Wu, Z.; Cohen, R.E. More accurate generalized gradient approximation for solids. *Phys. Rev. B* **2006**, *73*, 235116. [CrossRef]

17. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. *Phys. Rev. Lett.* **1996**, *77*, 3865–3868. [CrossRef]

18. Perdew, P.; Ruzsinszky, A.; Csonka, G.I.; Vydrov, O.A.; Scuseria, G.E.; Constantin, L.A.; Zhou, X.; Burke, K. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces. *Phys. Rev. Lett.* **2008**, *100*, 136406. [CrossRef]

19. Heyd, J.; Scuseria, G.E.; Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. *J. Chem. Phys.* **2003**, *118*, 8207–8215. [CrossRef]

20. Schwarz, K.; Blaha, P.; Madsen, G.K.H. Electronic structure calculations of solids using the WIEN2k package for material sciences. *Comput. Phys. Commun.* **2002**, *147*, 71–76.

21. Blaha, P.; Schwarz, K.; Tran, F.; Laskowski, R.; Madsen, G.K.H.; Marks, L.D. WIEN2k, an APW + LO Program for Calculating The Properties of Solids. *J. Chem. Phys.* **2020**, *152*, 074101. [CrossRef] [PubMed]

22. Madsen, G.K.H.; Carrete, J.; Verstraete, M.J. BoltzTraP2, a program for interpolating band structures and calculating semi-classical transport coefficients. *Computer. Phys. Commun.* **2018**, *231*, 140–145. [CrossRef]

23. Bardeen, J.; Shockley, W. Deformation Potentials and Mobilities in Non-Polar Crystals. *Phys. Rev.* **1950**, *80*, 72–80. [CrossRef]

24. Herring, C.; Vogt, E. Transport and Deformation-Potential Theory for Many-Valley Semiconductors with Anisotropic Scattering. *Phys. Rev.* **1956**, *101*, 944–961. [CrossRef]

25. Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G.L.; Cococcioni, M.; Dabo, I.; et al. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. *J. Phys.* **2009**, *21*, 395502. [CrossRef] [PubMed]

26. Togo, A.; Chaput, L.; Tanaka, I. Distributions of phonon lifetimes in Brillouin zones. *Phys. Rev. B* **2015**, *91*, 094306. [CrossRef]

27. Togo, A. First-principles Phonon Calculations with Phonopy and Phono3py. *J. Phys. Soc. Jpn.* **2023**, *92*, 012001. [CrossRef]

28. Matsunaga, T.; Yamada, N.; Kubota, Y. Structures of stable and metastable Ge₂Sb₂Te₅, an intermetallic compound in GeTe–Sb₂Te₃ pseudobinary systems. *Acta Cryst.* **2004**, *B60*, 685–691. [CrossRef]

29. Kooi, B.J.; De Hosson, J.T.M. Electron diffraction and high-resolution transmission electron microscopy of the high temperature crystal structures of GeₓSb₂Te₃ₓ (x=1,2,3) phase change material. *J. Appl. Phys.* **2002**, *92*, 3584–3590. [CrossRef]

30. Petrov, I.I.; Imamov, R.M.; Pinsker, Z.G. Electron-diffraction determination of the structures of Ge₂Sb₂Te₅ and GeSb₄Te₇. *Sov. Phys. Crystallogr.* **1968**, *13*, 339–342.

31. Lee, B.-S.; Abelson, J.R.; Bishop, S.G.; Kang, D.-H.; Cheong, B.-K.; Kim, K.-B. Investigation of the optical and electronic properties of Ge₂Sb₂Te₅ phase change material in its amorphous, cubic, and hexagonal phases. *J. Appl. Phys.* **2005**, *97*, 093509. [CrossRef]

32. Kato, T.; Tanaka, K. Electronic Properties of Amorphous and Crystalline Ge₂Sb₂Te₅ Films. *Jpn. J. Appl. Phys.* **2005**, *44*, 7430. [CrossRef]

33. Park, J.-W.; Baek, S.H.; Kang, T.D.; Lee, H.; Kang, Y.-S.; Lee, T.-Y.; Suh, D.-S.; Kim, K.J.; Kim, C.K.; Khang, Y.H.; et al. Optical properties of (GeTe, Sb₂Te₃) pseudobinary thin films studied with spectroscopic ellipsometry. *Appl. Phys. Lett.* **2008**, *93*, 021914. [CrossRef]

34. Lee, G.; Jhi, S.-H. Ab initio studies of structural and electronic properties of the crystalline Ge₂Sb₂Te₅. *Phys. Rev. B* **2008**, *77*, 153201. [CrossRef]

35. Sa, B.; Zhou, J.; Ahuja, R.; Sun, Z. First-principles investigations of electronic and mechanical properties for stable Ge₂Sb₂Te₅ with van der Waals corrections. *Comput. Mater. Sci.* **2014**, *82*, 66–69. [CrossRef]

36. Espinosa, E.; Alkorta, I.; Elguero, J.; Molins, E. From weak to strong interactions: A comprehensive analysis of the topological and energetic properties of the electron density distribution involving X–H…F–Y systems. *J. Chem. Phys.* **2002**, *117*, 5529–5542. [CrossRef]

37. Yang, H.; Boulet, P.; Record, M.-C. A rapid method for analyzing the chemical bond from energy densities calculations at the bond critical point. *Comput. Theor. Chem.* **2020**, *1178*, 112784. [CrossRef]

38. Yang, H.; Boulet, P.; Record, M.-C. New insight into the structure-property relationships from chemical bonding analysis: Application to thermoelectric materials. *J. Solid State Chem.* **2020**, *286*, 121266. [CrossRef]

39. Ma, W.; Record, M.-C.; Tian, J.; Boulet, P. Influence of the stacking sequence on layered-chalcogenide properties: First principles investigation of Pb₂Bi₂Te₅. *Phys. Chem. Chem. Phys.* **2021**, *23*, 11300–11313. [CrossRef] [PubMed]
11. Du, W.; Gu, Y.; Wang, K.; Yang, X.; Xing, J.; Guo, K.; Luo, J.; Zhao, J.-T. Effective Mass Enhancement and Thermal Conductivity Reduction for Improving the Thermoelectric Properties of Pseudo-Binary Ge₂Sb₂Te₅. *Ann. Phys.* **2020**, *532*, 1900390. [CrossRef]

12. Hu, P.; Wei, T.-R.; Huang, S.-J.; Xia, X.-G.; Qiu, P.-F.; Yang, J.; Chen, L.-D.; Shi, X. Anion-site-modulated thermoelectric properties in Ge₂Sb₂Te₅-based compounds. *Rare Met.* **2020**, *39*, 1127–1133. [CrossRef]

13. Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. *Phys. Rev.* **1964**, *136*, B86–B871. [CrossRef]

14. Kohn, W.; Sham, L.J. Self-Consistent Equations Including Exchange and Correlation Effects. *Phys. Rev.* **1965**, *140*, A1133–A1138. [CrossRef]

15. Perdew, J.P.; Zunger, A. Self-interaction correction to density-functional approximations for many-electron systems. *Phys. Rev. B* **1981**, *23*, 5048–5079. [CrossRef]

16. Wu, Z.; Cohen, R.E. More accurate generalized gradient approximation for solids. *Phys. Rev. B* **2006**, *73*, 235116. [CrossRef]

17. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. *Phys. Rev. Lett.* **1996**, *77*, 3865–3868. [CrossRef]

18. Perdew, P.; Ruzsinszky, A.; Csonka, G.I.; Vydrov, O.A.; Scuseria, G.E.; Constantin, L.A.; Zhou, X.; Burke, K. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces. *Phys. Rev. Lett.* **2008**, *100*, 136406. [CrossRef]

19. Heyd, J.; Scuseria, G.E.; Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. *J. Chem. Phys.* **2003**, *118*, 8207–8215. [CrossRef]

20. Schwarz, K.; Blaha, P.; Madsen, G.K.H. Electronic structure calculations of solids using the WIEN2k package for material sciences. *Comput. Phys. Commun.* **2002**, *147*, 71–76.

21. Blaha, P.; Schwarz, K.; Tran, F.; Laskowski, R.; Madsen, G.K.H.; Marks, L.D. WIEN2k, an APW + LO Program for Calculating The Properties of Solids. *J. Chem. Phys.* **2020**, *152*, 074101. [CrossRef] [PubMed]

22. Madsen, G.K.H.; Carrete, J.; Verstraete, M.J. BoltzTraP2, a program for interpolating band structures and calculating semi-classical transport coefficients. *Computer. Phys. Commun.* **2018**, *231*, 140–145. [CrossRef]

23. Bardeen, J.; Shockley, W. Deformation Potentials and Mobilities in Non-Polar Crystals. *Phys. Rev.* **1950**, *80*, 72–80. [CrossRef]

24. Herring, C.; Vogt, E. Transport and Deformation-Potential Theory for Many-Valley Semiconductors with Anisotropic Scattering. *Phys. Rev.* **1956**, *101*, 944–961. [CrossRef]

25. Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G.L.; Cococcioni, M.; Dabo, I.; et al. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. *J. Phys.* **2009**, *21*, 395502. [CrossRef] [PubMed]

26. Togo, A.; Chaput, L.; Tanaka, I. Distributions of phonon lifetimes in Brillouin zones. *Phys. Rev. B* **2015**, *91*, 094306. [CrossRef]

27. Togo, A. First-principles Phonon Calculations with Phonopy and Phono3py. *J. Phys. Soc. Jpn.* **2023**, *92*, 012001. [CrossRef]

28. Matsunaga, T.; Yamada, N.; Kubota, Y. Structures of stable and metastable Ge₂Sb₂Te₅, an intermetallic compound in GeTe–Sb₂Te₃ pseudobinary systems. *Acta Cryst.* **2004**, *B60*, 685–691. [CrossRef]

29. Kooi, B.J.; De Hosson, J.T.M. Electron diffraction and high-resolution transmission electron microscopy of the high temperature crystal structures of GeₓSb₂Te₃ₓ (x=1,2,3) phase change material. *J. Appl. Phys.* **2002**, *92*, 3584–3590. [CrossRef]

30. Petrov, I.I.; Imamov, R.M.; Pinsker, Z.G. Electron-diffraction determination of the structures of Ge₂Sb₂Te₅ and GeSb₄Te₇. *Sov. Phys. Crystallogr.* **1968**, *13*, 339–342.

31. Lee, B.-S.; Abelson, J.R.; Bishop, S.G.; Kang, D.-H.; Cheong, B.-K.; Kim, K.-B. Investigation of the optical and electronic properties of Ge₂Sb₂Te₅ phase change material in its amorphous, cubic, and hexagonal phases. *J. Appl. Phys.* **2005**, *97*, 093509. [CrossRef]

32. Kato, T.; Tanaka, K. Electronic Properties of Amorphous and Crystalline Ge₂Sb₂Te₅ Films. *Jpn. J. Appl. Phys.* **2005**, *44*, 7430. [CrossRef]

33. Park, J.-W.; Baek, S.H.; Kang, T.D.; Lee, H.; Kang, Y.-S.; Lee, T.-Y.; Suh, D.-S.; Kim, K.J.; Kim, C.K.; Khang, Y.H.; et al. Optical properties of (GeTe, Sb₂Te₃) pseudobinary thin films studied with spectroscopic ellipsometry. *Appl. Phys. Lett.* **2008**, *93*, 021914. [CrossRef]

34. Lee, G.; Jhi, S.-H. Ab initio studies of structural and electronic properties of the crystalline Ge₂Sb₂Te₅. *Phys. Rev. B* **2008**, *77*, 153201. [CrossRef]

35. Sa, B.; Zhou, J.; Ahuja, R.; Sun, Z. First-principles investigations of electronic and mechanical properties for stable Ge₂Sb₂Te₅ with van der Waals corrections. *Comput. Mater. Sci.* **2014**, *82*, 66–69. [CrossRef]

36. Espinosa, E.; Alkorta, I.; Elguero, J.; Molins, E. From weak to strong interactions: A comprehensive analysis of the topological and energetic properties of the electron density distribution involving X–H…F–Y systems. *J. Chem. Phys.* **2002**, *117*, 5529–5542. [CrossRef]

37. Yang, H.; Boulet, P.; Record, M.-C. A rapid method for analyzing the chemical bond from energy densities calculations at the bond critical point. *Comput. Theor. Chem.* **2020**, *1178*, 112784. [CrossRef]

38. Yang, H.; Boulet, P.; Record, M.-C. New insight into the structure-property relationships from chemical bonding analysis: Application to thermoelectric materials. *J. Solid State Chem.* **2020**, *286*, 121266. [CrossRef]

39. Ma, W.; Record, M.-C.; Tian, J.; Boulet, P. Influence of the stacking sequence on layered-chalcogenide properties: First principles investigation of Pb₂Bi₂Te₅. *Phys. Chem. Chem. Phys.* **2021**, *23*, 11300–11313. [CrossRef] [PubMed]
40. Delaire, O.; Ma, J.; Marty, K.; May, A.F.; McGuire, M.A.; Du, M.-H.; Singh, D.J.; Podlesnyak, A.; Ehlers, G.; Lumsden, M.D.; et al. Giant anharmonic phonon scattering in PbTe. *Nature Mater.* **2011**, *10*, 614–619. [CrossRef]

41. Xie, H. The role of off-centering behavior and acoustic-optical phonon coupling in heat transport. *Mater. Lab* **2023**, *2*, 220051. [CrossRef]

**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.
40. Delaire, O.; Ma, J.; Marty, K.; May, A.F.; McGuire, M.A.; Du, M.-H.; Singh, D.J.; Podlesnyak, A.; Ehlers, G.; Lumsden, M.D.; et al. Giant anharmonic phonon scattering in PbTe. *Nature Mater.* **2011**, *10*, 614–619. [CrossRef]

41. Xie, H. The role of off-centering behavior and acoustic-optical phonon coupling in heat transport. *Mater. Lab* **2023**, *2*, 220051. [CrossRef]

**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.