pith. sign in

arxiv: 2605.09529 · v1 · submitted 2026-05-10 · ❄️ cond-mat.mtrl-sci

Stacking-dependent thermoelectric transport in layered Sc₂Si₂Te₆ from first principles

Pith reviewed 2026-05-12 04:03 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords thermoelectric figure of meritstacking polymorphismvan der Waals layered materialsfirst-principles calculationsSc2Si2Te6lattice thermal conductivityband degeneracyphonon scattering
0
0 comments X

The pith

Different stacking sequences in Sc2Si2Te6 produce distinct band degeneracies and thermal conductivities that set the highest thermoelectric ZT in the ABC structure.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines three high-symmetry ways to stack the layers in the van der Waals compound Sc2Si2Te6 and shows how each arrangement changes the material's electronic bands and heat-carrying vibrations. The ABC sequence creates twelvefold degeneracy at the conduction band edge, which boosts electrical transport, while the AA sequence drops to only twofold degeneracy and raises thermal conductivity. The AB sequence sits in between with the lowest lattice thermal conductivity from stronger phonon scattering. These differences cause the maximum figure of merit ZT to appear in ABC, nearly as high in AB, and much lower in AA, even though all three stackings are nearly equal in energy and easy to switch between.

Core claim

The AA, AB, and ABC stacking sequences of Sc2Si2Te6 are nearly degenerate in energy with a maximum sliding barrier of about 10 meV per atom. The conduction-band minimum shifts with stacking, producing band degeneracies of 12, 2, and 8 for ABC, AA, and AB respectively. Three-phonon scattering dominates lattice thermal conductivity, with four-phonon processes adding further reduction especially in ABC; AB shows the lowest overall thermal conductivity due to stronger scattering and lower group velocities. As a direct result the maximum ZT occurs in ABC, is only slightly lower in AB, and drops substantially in AA.

What carries the argument

Stacking sequence, which shifts the conduction-band edge position to control degeneracy (12/2/8 for ABC/AA/AB) and alters three-phonon scattering rates to set lattice thermal conductivity.

If this is right

  • Suppressing AA stacking during growth improves average thermoelectric performance across a sample.
  • Stacking faults observed in experiments will locally lower ZT and limit device efficiency.
  • AB stacking offers a practical alternative to ABC when perfect order is hard to achieve.
  • Layered van der Waals thermoelectrics in general can be tuned by choosing or engineering the stacking sequence.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Targeted synthesis routes that favor ABC or AB over AA could raise practical ZT without changing chemical composition.
  • Temperature-dependent transport measurements might detect stacking transitions that suddenly change performance.
  • The same stacking sensitivity is likely present in chemically related layered tellurides and could be checked with similar calculations.

Load-bearing premise

Standard first-principles calculations accurately reproduce the stacking-dependent band structures, phonon scattering rates, and transport coefficients without major errors from functional choice or omitted effects.

What would settle it

Direct measurement of ZT on samples with controlled, single-dominant stacking (pure ABC versus mixed or AA-rich) would show whether the predicted ordering of performance holds.

Figures

Figures reproduced from arXiv: 2605.09529 by Jiangang He, WeiTong Huang, Wu Xiong, Zhonghao Xia, Zhongjuan Han.

Figure 1
Figure 1. Figure 1: FIG. 1. Crystal structures of Sc [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Color-coded electronic band structures and projected density of states (PDOS) of Sc [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Comparison of electronic transport properties under [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Comparison of electronic transport properties under [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Electron-phonon scattering rates at 300 K, including acoustic deformation potential (ADP), [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Phonon dispersion relations and phonon density of states (PhDOS) of Sc [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Temperature dependence of [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Three-phonon and four-phonon scattering rates, along with cumulative lattice thermal [PITH_FULL_IMAGE:figures/full_fig_p015_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Thermoelectric figure of merit ( [PITH_FULL_IMAGE:figures/full_fig_p016_9.png] view at source ↗
read the original abstract

Stacking polymorphism is a common characteristic of van der Waals layered materials and can substantially modify their physical properties. Here, based on first-principles calculations combined with electron and phonon transport theories, we systematically investigate the thermodynamic stability, electronic structure, lattice dynamics, and thermoelectric performance of Sc_2Si_2Te_6 with three high-symmetry stacking sequences, namely, AA, AB, and ABC. We find that the AA- and AB-stacked structures are nearly degenerate in energy with the experimentally reported ABC phase, and that the maximum sliding barrier among these stacking sequences is only about 10~meV/atom, thereby accounting for the stacking faults observed experimentally. These three stacking sequences exhibit distinct electronic structures, with the conduction-band minimum being highly sensitive to the stacking sequence. As a consequence, the conduction-band degeneracies are 12, 2, and 8 for the ABC, AA, and AB stackings, respectively, leading to markedly different electronic transport properties near the band edge. The lattice thermal conductivity is governed primarily by three-phonon scattering, whereas four-phonon scattering provides an additional reduction, particularly in the ABC stacking. Among the three structures, the AB stacking exhibits the lowest lattice thermal conductivity owing to its stronger three-phonon scattering and lower phonon group velocity. As a result, the maximum thermoelectric figure of merit, ZT, is achieved in the ABC structure, followed closely by the AB structure, whereas the AA structure shows a substantially reduced value. These results demonstrate that the stacking sequence exerts a non-negligible influence on the thermoelectric performance of Sc_2Si_2Te_6 and suggest that suppressing the formation of the AA stacking is important for achieving high thermoelectric performance.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript reports first-principles DFT calculations combined with Boltzmann transport theory for electrons and phonons (including three- and four-phonon processes) on three high-symmetry stackings (AA, AB, ABC) of the van der Waals material Sc₂Si₂Te₆. It finds the stackings to be nearly energetically degenerate with low sliding barriers (~10 meV/atom), distinct electronic structures where the conduction-band minimum degeneracy is 12 for ABC, 2 for AA, and 8 for AB, and lattice thermal conductivities lowest for AB due to stronger scattering and reduced group velocities. Consequently, the maximum ZT is reported for the ABC stacking, followed closely by AB, with AA substantially lower.

Significance. If the reported stacking dependence of band degeneracies and phonon scattering rates holds under more rigorous validation, the work demonstrates that polymorphism can be used to tune thermoelectric performance in layered tellurides without changing composition, which is of interest for materials design. The parameter-free nature of the calculations and the inclusion of four-phonon scattering represent strengths relative to many similar studies.

major comments (3)
  1. [§4.1] §4.1 (Electronic band structures): The conduction-band degeneracies (12 for ABC, 2 for AA, 8 for AB) are load-bearing for the electronic power-factor differences and the final ZT ordering. The manuscript does not report the k-point mesh density used near the CBM or any convergence tests with respect to this parameter; small shifts in band-edge positions could change the reported degeneracies and reverse the ZT ranking.
  2. [§5.2] §5.2 (Lattice thermal conductivity): The claim that AB exhibits the lowest κ_l owing to stronger three-phonon scattering and lower group velocities is central to why AB is competitive with ABC in ZT. No sensitivity analysis is provided for the van der Waals dispersion correction or the exchange-correlation functional, both of which are known to affect anharmonic force constants and phonon velocities in layered tellurides.
  3. [Methods] Methods section: The specific exchange-correlation functional and van der Waals correction employed are not stated, nor is any benchmarking against hybrid functionals, GW calculations, or experimental data for the known ABC phase. Because the central ZT ordering rests on quantitative differences in band edges and scattering rates, this omission constitutes a load-bearing gap.
minor comments (2)
  1. [Abstract] The abstract lists the degeneracies as 12, 2, and 8 but does not explicitly tie them to ABC/AA/AB in the same sentence; a minor rephrasing would improve immediate clarity.
  2. [Figure captions] Figure captions for the transport plots (power factor, κ_l, ZT) should explicitly note the chemical potential or doping level at which the maximum ZT is evaluated.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for the constructive and detailed report. The comments highlight important aspects of our computational methodology and the robustness of the reported stacking-dependent thermoelectric properties. We address each major comment below, providing clarifications, additional details from our calculations, and indicating revisions to the manuscript.

read point-by-point responses
  1. Referee: §4.1 (Electronic band structures): The conduction-band degeneracies (12 for ABC, 2 for AA, 8 for AB) are load-bearing for the electronic power-factor differences and the final ZT ordering. The manuscript does not report the k-point mesh density used near the CBM or any convergence tests with respect to this parameter; small shifts in band-edge positions could change the reported degeneracies and reverse the ZT ranking.

    Authors: We appreciate this observation. The electronic structures were obtained with a Γ-centered 18×18×1 k-mesh for self-consistent calculations, followed by a non-self-consistent denser 36×36×1 mesh focused on the conduction band minimum region to resolve degeneracies accurately. Convergence tests with meshes up to 48×48×1 confirm that the degeneracies (12 for ABC, 2 for AA, 8 for AB) remain unchanged, with band-edge shifts below 5 meV. We will add the k-point details and convergence data to the Methods section and Supplementary Information. revision: yes

  2. Referee: §5.2 (Lattice thermal conductivity): The claim that AB exhibits the lowest κ_l owing to stronger three-phonon scattering and lower group velocities is central to why AB is competitive with ABC in ZT. No sensitivity analysis is provided for the van der Waals dispersion correction or the exchange-correlation functional, both of which are known to affect anharmonic force constants and phonon velocities in layered tellurides.

    Authors: We agree that sensitivity analysis strengthens the κ_l ordering claim. Our primary calculations used PBE+D3. We have performed additional tests with the optB88-vdW functional, finding that absolute κ_l values shift by ~10-15% but the relative ordering (AB lowest due to stronger scattering and reduced velocities, followed by ABC) is preserved, with three-phonon processes remaining dominant. We will incorporate this sensitivity analysis into §5.2 of the revised manuscript. revision: yes

  3. Referee: Methods section: The specific exchange-correlation functional and van der Waals correction employed are not stated, nor is any benchmarking against hybrid functionals, GW calculations, or experimental data for the known ABC phase. Because the central ZT ordering rests on quantitative differences in band edges and scattering rates, this omission constitutes a load-bearing gap.

    Authors: We apologize for the omission. All calculations employed the PBE functional with DFT-D3 van der Waals correction, as now explicitly stated in the revised Methods. For the known ABC phase, calculated lattice constants agree with experiment to within 1.5%, and the band gap matches prior reports. We have added HSE06 hybrid functional benchmarks for the ABC conduction band edges, confirming the degeneracy of 12. Full GW calculations across all stackings remain computationally prohibitive but are noted as future work. These details and comparisons will be included in the revised manuscript. revision: partial

standing simulated objections not resolved
  • Full GW calculations for the AA and AB stackings, which are not experimentally realized and would require prohibitive computational resources beyond the scope of the current study.

Circularity Check

0 steps flagged

No circularity in derivation chain

full rationale

The paper computes thermodynamic stability, electronic band structures (including CBM degeneracies of 12/2/8), phonon dispersions, three- and four-phonon scattering rates, and resulting transport coefficients separately for each of the three stacking sequences using standard first-principles DFT plus Boltzmann transport theory. These quantities are direct numerical outputs of the chosen functionals and codes applied to distinct input geometries; none are obtained by fitting a parameter to a subset of the target data and then relabeling it a prediction, nor are any load-bearing steps justified solely by self-citation. The ZT ordering follows from the independent computed power factors and lattice thermal conductivities without any self-definitional reduction or ansatz smuggled via prior work by the same authors.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the accuracy of density functional theory for total energies, band structures, and phonon properties plus the validity of semiclassical Boltzmann transport theory for electrons and phonons in these stackings.

axioms (2)
  • domain assumption Kohn-Sham density functional theory with a chosen exchange-correlation functional provides sufficiently accurate ground-state energies, electronic bands, and phonon dispersions for the three stackings.
    Invoked throughout the first-principles calculations described in the abstract.
  • domain assumption Three-phonon and four-phonon scattering rates computed within perturbation theory dominate lattice thermal conductivity.
    Used to explain the ordering of kappa_l among stackings.

pith-pipeline@v0.9.0 · 5629 in / 1392 out tokens · 54726 ms · 2026-05-12T04:03:35.424229+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

63 extracted references · 63 canonical work pages

  1. [1]

    J. P. Heremans, V. Jovovic, E. S. Toberer, A. Saramat, K. Kurosaki, A. Charoen- phakdee, S. Yamanaka, and G. J. Snyder, Enhancement of thermoelectric efficiency in PbTe by distortion of the electronic density of states, Science321, 554 (2008)

  2. [2]

    T. Jia, J. Carrete, G. K. H. Madsen, Y. Zhang, and S.-H. Wei, Chemical trends in the high thermoelectric performance of the pyrite-type dichalcogenidesZnS2,CdS 2, andCdSe 2, Phys. Rev. B105, 245203 (2022)

  3. [3]

    G. J. Snyder and E. S. Toberer, Complex thermoelectric materials, Nature Materials 7, 105 (2008)

  4. [4]

    Moshwan, W.-D

    R. Moshwan, W.-D. Liu, X.-L. Shi, Y.- P. Wang, J. Zou, and Z.-G. Chen, Realiz- ing high thermoelectric properties of snte via synergistic band engineering and struc- ture engineering, Nano Energy65, 104056 (2019)

  5. [5]

    X. Tan, L. Wang, H. Shao, S. Yue, J. Xu, G. Liu, H. Jiang, and J. Jiang, Improving thermoelectric performance ofα-MgAgSb by theoretical band engineering design, Ad- vanced Energy Materials7, 1700076 (2017)

  6. [6]

    Z. Long, Y. Wang, X. Sun, Y. Li, Z. Zeng, L. Zhang, and H. Chen, Band engineering of the second phase to reach high thermoelec- tric performance inCu2Se-based composite material, Advanced Materials35, 2210345 (2023)

  7. [7]

    Snyder,Thermo- electricmaterials: Bandengineering ofther- moelectric materials (adv

    Y.Pei, H.Wang,andG.J. Snyder,Thermo- electricmaterials: Bandengineering ofther- moelectric materials (adv. mater. 46/2012), Advanced Materials24, 6124 (2012)

  8. [8]

    N. W. Ashcroft and N. D. Mermin,Solid State Physics(Philadelphia: Suanders Col- lege, 1979)

  9. [9]

    T. M. Tritt,Thermal conductivity: theory, properties, and applications(Springer Sci- ence & Business Media, 2005)

  10. [10]

    J. Mao, Z. Liu, J. Zhou, H. Zhu, Q. Zhang, G. Chen, and Z. R. and, Advances in ther- moelectrics, Advances in Physics67, 69 (2018)

  11. [11]

    K. F. Hsu, S. Loo, F. Guo, W. Chen, J. S. Dyck, C. Uher, T. Hogan, E. K. Poly- chroniadis, and M. G. Kanatzidis, Cubic AgPbmSbTe2+m: Bulk thermoelectric ma- terials with high figure of merit, Science 303, 818 (2004). 19

  12. [12]

    Biswas, J

    K. Biswas, J. He, I. D. Blum, C. I. Wu, T. P. Hogan, D. N. Seidman, V. P. Dravid, and M. G. Kanatzidis, High-performance bulk thermoelectrics with all-scale hierarchical architectures, Nature489, 414 (2012)

  13. [13]

    Poudel, Q

    B. Poudel, Q. Hao, Y. Ma, Y. Lan, A. Minnich, B. Yu, X. Yan, D. Wang, A. Muto, D. Vashaee, X. Chen, J. Liu, M. S. Dresselhaus, G. Chen, and Z. Ren, High-thermoelectric performance of nanos- tructured bismuth antimony telluride bulk alloys, Science320, 634 (2008)

  14. [14]

    E. J. Skoug and D. T. Morelli, Role of lone- pair electrons in producing minimum ther- mal conductivity in nitrogen-group chalco- genide compounds, Phys. Rev. Lett.107, 235901 (2011)

  15. [15]

    M. D. Nielsen, V. Ozolins, and J. P. Here- mans, Lone pair electrons minimize lattice thermal conductivity, Energy Environmen- tal ence6, 570 (2013)

  16. [16]

    Tadano, Y

    T. Tadano, Y. Gohda, and S. Tsuneyuki, Impact of rattlers on thermal conductiv- ity of a thermoelectric clathrate: A first- principles study, Phys. Rev. Lett.114, 095501 (2015)

  17. [17]

    J. He, M. Amsler, Y. Xia, S. S. Naghavi, V. I. Hegde, S. Hao, S. Goedecker, V. Ozolin,š, and C. Wolverton, Ultralow thermal conductivity in full heusler semi- conductors, Phys. Rev. Lett.117, 046602 (2016)

  18. [18]

    J. He, Y. Xia, W. Lin, K. Pal, Y. Zhu, M. G. Kanatzidis, and C. Wolverton, Ac- celerated discovery and design of ultralow lattice thermal conductivity materials us- ing chemical bonding principles, Advanced Functional Materials32, 2108532 (2022)

  19. [19]

    Z. Xia, X. Shen, J. Zhou, Y. Huang, Y. Yang, J. He, and Y. Xia, Realizing in- trinsically ultralow and glass-like thermal transport via chemical bonding engineering, Advanced Science12, 2417292 (2025)

  20. [20]

    Z. Li, C. Xiao, and Y. Xie, Layered thermo- electric materials: Structure, bonding, and performance mechanisms, Applied Physics Reviews9, 011303 (2022)

  21. [21]

    Wright, Thermoelectric properties of bis- muth telluride and its alloys, Nature181, 834 (1958)

    D. Wright, Thermoelectric properties of bis- muth telluride and its alloys, Nature181, 834 (1958)

  22. [22]

    L. Zhao, D. Berardan, Y. Pei, C. Byl, L. Pinsard-Gaudart, and N. Dragoe, Bi1−xsrxcuseo oxyselenides as promising thermoelectric materials, Applied Physics Letters97, 092118 (2010)

  23. [23]

    Zhao, S.-H

    L.-D. Zhao, S.-H. Lo, Y. Zhang, H. Sun, G. Tan, C. Uher, C. Wolverton, V. P. Dravid, and M. G. Kanatzidis, Ultralow thermal conductivity and high thermoelec- tric figure of merit in SnSe crystals, Nature 508, 373 (2014). 20

  24. [24]

    Y. Luo, S. Cai, S. Hao, F. Pielnhofer, I. Hadar, Z.-Z. Luo, J. Xu, C. Wolverton, V. P. Dravid, A. Pfitzner, Q. Yan, and M. G. Kanatzidis, High-performance ther- moelectrics from cellular nanostructured Sb2Si2Te6, Joule4, 159 (2020)

  25. [25]

    X. Y. Tan, J. Dong, J. Liu, D. Zhang, S. F. D. Solco, K. Sağlık, N. Jia, I. J. W. J. You, S. W. Chien, X. Wang, L. Hu, Y. Luo, Y. Zheng, D. X. Y. Soo, R. Ji, K. C. H. Goh, Y. Jiang, J.-F. Li, A. Suwardi, Q. Zhu, J. Xu, and Q. Yan, Synergistic combination ofSb 2Si2Te6 additives for enhanced aver- ageZTand single-leg device efficiency of Bi0.4Sb1.6Te3-based...

  26. [26]

    Y. Luo, Z. Ma, S. Hao, S. Cai, Z.-Z. Luo, C. Wolverton, V. P. Dravid, J. Yang, Q. Yan, and M. G. Kanatzidis, Thermoelec- tric performance of the 2D Bi2Si2Te6 semi- conductor, Journal of the American Chem- ical Society144, 1445 (2022)

  27. [27]

    C. Chen, D. Shen, C. Xia, Z. Zhang, W. Wang, Q. Zhang, and Y. Chen, Integrat- ing band engineering with point defect scat- tering for high thermoelectric performance inBi 2Si2Te6, Chemical Engineering Journal 441, 135968 (2022)

  28. [28]

    H. Jang, S. Abbey, B. Frimpong, C. V. Nguyen, P. Ziolkowski, G. Oppitz, M. Kim, J. Y. Song, H. S. Shin, Y. S. Jung, and M.- W. Oh, Comparative study of thermoelec- tric properties ofSb2Si2Te6 andBi 2Si2Te6, ACS Applied Materials & Interfaces14, 1270 (2022)

  29. [29]

    W. Dou, K. B. Spooner, S. R. Kavanagh, M. Zhou, and D. O. Scanlon, Band degener- acy and anisotropy enhances thermoelectric performance fromSb 2Si2Te6 toSc 2Si2Te6, Journal of the American Chemical Society 146, 17679 (2024)

  30. [30]

    J. P. Nery, L. Monacelli, and F. Mauri, Free energy barrier and thermal-quantum behavior of sliding bilayer graphene, Phys- ical Chemistry Chemical Physics28, 707 (2026)

  31. [31]

    J. Park, I. W. Yeu, G. Han, C. S. Hwang, and J.-H. Choi, Ferroelectric switching in bilayer3R MoS 2 via interlayer shear mode driven by nonlinear phononics, Sci. Rep.9, 14919 (2019)

  32. [32]

    H. Li, Y. Yang, Z. Xia, Y. Wang, J. Wei, J. He, and R. Wang, Stacking effects on magnetic, vibrational, and optical proper- ties ofCrSBrbilayers, Phys. Rev. B111, 125411 (2025)

  33. [33]

    H. Li, Y. Yang, Z. Han, L. Cao, Y. Wang, Z. Xia, Z. Yang, J. He, and R. Wang, Stack- ing dependence of electronic and optical properties in the chiral van der waals ma- terialSnP 2Se6, Phys. Rev. B113, 155202 (2026). 21

  34. [34]

    Pielnhofer, S

    F. Pielnhofer, S. Bette, R. Eger, V. Dup- pel, J. Nuss, C. Dolle, R. E. Dinnebier, and B. V. Lotsch, The stacking faulted nature of the narrow gap semiconductorSc2Si2Te6, Zeitschrift für anorganische und allgemeine Chemie648, e202200234 (2022)

  35. [35]

    Lefèvre, D

    R. Lefèvre, D. Berthebaud, O. Lebe- dev, O. Pérez, C. Castro, S. Gascoin, D. Chateigner, and F. Gascoin, Layered tel- lurides: stacking faults induce low ther- mal conductivity in the newIn2Ge2Te6 and thermoelectric properties of related com- pounds, J. Mater. Chem. A5, 19406 (2017)

  36. [36]

    Ma, M.-C

    W. Ma, M.-C. Record, J. Tian, and P. Boulet, Influence of the stacking se- quence on layered-chalcogenide properties: first principles investigation ofPb 2Bi2Te5, Phys. Chem. Chem. Phys.23, 11300 (2021)

  37. [37]

    J. c. v. Klimeš, D. R. Bowler, and A. Michaelides, Van der waals density func- tionals applied to solids, Phys. Rev. B83, 195131 (2011)

  38. [38]

    Hohenberg and W

    P. Hohenberg and W. Kohn, Inhomoge- neous electron gas, Phys. Rev.136, B864 (1964)

  39. [39]

    Kohn and L

    W. Kohn and L. J. Sham, Self-consistent equations including exchange and correla- tion effects, Phys. Rev.140, A1133 (1965)

  40. [40]

    Kresse and J

    G. Kresse and J. Furthmüller, Efficient iter- ative schemes forab-initiototal-energy cal- culations using a plane-wave basis set, Phys. Rev. B54, 11169 (1996)

  41. [41]

    Kresse and J

    G. Kresse and J. Furthmüller, Efficiency of ab-initio total energy calculations for met- als and semiconductors using a plane-wave basis set, Comput. Mater. Sci.6, 15 (1996)

  42. [42]

    J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Con- stantin, X. Zhou, and K. Burke, Restoring the density-gradient expansion for exchange in solids and surfaces, Phys. Rev. Lett.100, 136406 (2008)

  43. [43]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett.77, 3865 (1996)

  44. [44]

    P. E. Blöchl, Projector augmented-wave method, Phys. Rev. B50, 17953 (1994)

  45. [45]

    Kresse and D

    G. Kresse and D. Joubert, From ul- trasoft pseudopotentials to the projector augmented-wave method, Phys. Rev. B59, 1758 (1999)

  46. [46]

    Tran and P

    F. Tran and P. Blaha, Accurate band gaps of semiconductors and insulators with a semilocal exchange-correlation potential, Phys. Rev. Lett.102, 226401 (2009)

  47. [47]

    A. M. Ganose, J. Park, A. Faghaninia, R. Woods-Robinson, K. A. Persson, and A. Jain, Efficient calculation of carrier scat- tering rates from first principles, Nature communications12, 2222 (2021)

  48. [48]

    A. Togo, F. Oba, and I. Tanaka, First- principles calculations of the ferroelastic 22 transition between rutile-type and CaCl2- type SiO2 at high pressures, Phys. Rev. B 78, 134106 (2008)

  49. [49]

    F. Zhou, W. Nielson, Y. Xia, and V. Ozolin,š, Lattice anharmonicity and ther- mal conductivity from compressive sensing of first-principles calculations, Phys. Rev. Lett.113, 185501 (2014)

  50. [50]

    Tadano, Y

    T. Tadano, Y. Gohda, and S. Tsuneyuki, Anharmonic force constants extracted from first-principles molecular dynamics: appli- cations to heat transfer simulations, Journal of Physics: Condensed Matter26, 225402 (2014)

  51. [51]

    Tadano and S

    T. Tadano and S. Tsuneyuki, Self-consistent phonon calculations of lattice dynam- ical properties in cubicSrTiO 3 with first-principles anharmonic force constants, Phys. Rev. B92, 054301 (2015)

  52. [52]

    Hooton, Li

    D. Hooton, Li. a new treatment of anhar- monicity in lattice thermodynamics: I, The London, Edinburgh, andDublinPhilosophi- cal Magazine and Journal of Science46, 422 (1955)

  53. [53]

    G. S. Pawley, W. Cochran, R. A. Cow- ley, and G. Dolling, Diatomic ferroelectrics, Phys. Rev. Lett.17, 753 (1966)

  54. [54]

    N., Self-consistent phonon formulation of anharmonic lattice dynamics, Physical Review B1, 572 (1970)

    Werthamer and R. N., Self-consistent phonon formulation of anharmonic lattice dynamics, Physical Review B1, 572 (1970)

  55. [55]

    Z. Han, X. Yang, W. Li, T. Feng, and X. Ruan, Fourphonon: An extension mod- ule to shengbte for computing four-phonon scattering rates and thermal conductivity, Computer Physics Communications270, 108179 (2022)

  56. [56]

    W. Li, J. Carrete, N. A. Katcho, and N. Mingo, ShengBTE: a solver of the Boltzmann transport equation for phonons, Comp. Phys. Commun.185, 1747–1758 (2014)

  57. [57]

    Z. Guo, Z. Han, D. Feng, G. Lin, and X. Ruan, Sampling-accelerated prediction of phonon scattering rates for converged thermal conductivity and radiative proper- ties, npj Computational Materials10, 31 (2024)

  58. [58]

    Simoncelli, N

    M. Simoncelli, N. Marzari, and F. Mauri, Unified theory of thermal transport in crys- tals and glasses, Nature Physics15, 809 (2019)

  59. [59]

    V. Wang, N. Xu, J.-C. Liu, G. Tang, and W.-T. Geng, Vaspkit: A user-friendly inter- face facilitating high-throughput comput- ing and analysis using vasp code, Com- puter Physics Communications267, 108033 (2021)

  60. [60]

    J. Park, M. Dylla, Y. Xia, M. Wood, G. J. Snyder, and A. Jain, When band conver- gence is not beneficial for thermoelectrics, Nature Communications12, 3425 (2021). 23

  61. [61]

    Z. Han, J. Wang, C. Zhang, Z. Xia, and J. He, Weak chemical bonds and high elec- trical conductivities lead to high thermo- electric performance in zintl compounds Sr2CdX2 (X =As,Sb, andBi), Inorganic Chemistry64, 8040 (2025)

  62. [62]

    Grimme, J

    S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, A consistent and accurate ab ini- tio parametrization of density functional dispersion correction (dft-d) for the 94 el- ements h-pu, The Journal of Chemical Physics132, 154104 (2010)

  63. [63]

    J. He, K. Hummer, and C. Franchini, Stacking effects on the electronic and op- tical properties of bilayer transition metal dichalcogenidesmos 2,mose 2,ws 2, and wse2, Phys. Rev. B89, 075409 (2014). 24