pith. sign in

arxiv: 2211.13346 · v2 · pith:DOE722FOnew · submitted 2022-11-23 · ❄️ cond-mat.mtrl-sci

Topological phonon analysis of the 2D buckled honeycomb lattice: an application to real materials

Pith reviewed 2026-05-24 10:21 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords topological phononsbuckled honeycomb2D materialsphonon topologytopological quantum chemistryforce constantsWilson loops
0
0 comments X

The pith

The buckled honeycomb lattice supports nine topological phonon phases in its analytical model, yet all real 2D materials fall into the trivial phase.

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

The paper determines how many distinct phonon band topologies are possible in the 2D buckled honeycomb lattice and which of them are topological. Group theory and topological quantum chemistry identify eleven phases with five necessarily non-trivial. An analytical model with force constants up to third neighbors finds four more topological phases via Wilson loops. First-principles calculations for Si, Ge, P, As and Sb place them all in the trivial phase. Monte Carlo sampling of the force constant space shows why topological phases are physically difficult to realize.

Core claim

Eleven distinct phases are possible, five of which necessarily have non-trivial topology according to topological quantum chemistry. Another four of them are also identified as topological using Wilson loops in an analytical model that includes all the symmetry allowed force constants up to third nearest neighbors, making a total of nine topological phases. All computed real crystals lie in a trivial phase. Monte Carlo calculations elucidate why topological phonon phases are physically difficult to realize in real materials with this crystal structure.

What carries the argument

Analytical model of all symmetry-allowed force constants up to third nearest neighbors, combined with Monte Carlo sampling of the resulting parameter space to determine accessible phases.

If this is right

  • Topological phonon phases are difficult to realize in the buckled honeycomb structure.
  • The phase diagram is dominated by the trivial phase for physically plausible force constants.
  • Real 2D crystals of Si, Ge, P, As, and Sb exhibit trivial phonon topology.
  • Five phases are guaranteed topological by symmetry alone, independent of specific force constants.

Where Pith is reading between the lines

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

  • Materials with this structure may require external tuning like strain to reach topological phases.
  • Analyses of other 2D lattices could reveal structures where topological phonons are easier to realize.
  • The truncation at third neighbors may miss effects that open topological regions in real systems.

Load-bearing premise

The analytical model truncated at third-nearest-neighbor force constants and the Monte Carlo sampling of that space are sufficient to determine which phases are physically accessible.

What would settle it

Finding a buckled honeycomb material whose phonon spectrum shows non-trivial topology, either through longer-range interactions or in a new compound.

Figures

Figures reproduced from arXiv: 2211.13346 by Ion Errea, J. L. Ma\~nes, Maia G. Vergniory, Martin Gutierrez-Amigo.

Figure 1
Figure 1. Figure 1: FIG. 1. Buckled honeycomb lattice with real an reciprocal [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Matrices of force constants consistent with the BHL symmetries up to third nearest neighbors. As described in [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Wilson loop analysis for phase 8 in Tables I and [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. DFPT phonon dispersion bands of P, Sb, As, Si and Ge in the BHL. The irrep content of each subset of bands [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Phonon dispersion bands for Germanium BHL (or [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Finally, a relation can be established between the ex￾istence of topological phases and the need for sizeable couplings between far apart neighbors. In any real mate￾rial the couplings between different neighbors are bound to decrease as the distance is increased, due to the lo￾calized nature of the atomic wave functions. In the case of the planar honeycomb lattice, where the boundaries between different p… view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. DFPT dispersion phonons for the “blue phosphorus” [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Symmetry operations of the [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Phase 1 with Winding 1 in the green subset. [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Phase 3 with Winding 1 in the red subset. [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Phase 4 with Winding 2 in red and blue subsets. [PITH_FULL_IMAGE:figures/full_fig_p013_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Phase 5 with Winding 2 in the red subset. [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Phase 6 without any winding [PITH_FULL_IMAGE:figures/full_fig_p013_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Phase 7 without any winding [PITH_FULL_IMAGE:figures/full_fig_p013_14.png] view at source ↗
Figure 19
Figure 19. Figure 19: FIG. 19. Non-zero force constants up to third nearest neigh [PITH_FULL_IMAGE:figures/full_fig_p014_19.png] view at source ↗
read the original abstract

By means of group theory, topological quantum chemistry, first-principles and Monte Carlo calculations, we analyze the topology of the 2D buckled honeycomb lattice phonon spectra. Taking the pure crystal structure as an input, we show that eleven distinct phases are possible, five of which necessarily have non-trivial topology according to topological quantum chemistry. Another four of them are also identified as topological using Wilson loops in an analytical model that includes all the symmetry allowed force constants up to third nearest neighbors, making a total of nine topological phases. We then compute the ab initio phonon spectra for the two-dimensional crystals of Si, Ge, P, As and Sb in this structure and construct its phase diagram. Despite the large proportion of topological phases found in the analytical model, all of the crystals lie in a trivial phase. By analyzing the force constants space using Monte Carlo calculations, we elucidate why topological phonon phases are physically difficult to realize in real materials with this crystal structure.

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

1 major / 1 minor

Summary. The manuscript analyzes the topology of phonon spectra on the 2D buckled honeycomb lattice via group theory, topological quantum chemistry (TQC), first-principles calculations, and Monte Carlo sampling. It identifies eleven symmetry-allowed phases, five of which are necessarily non-trivial by TQC. An analytical model with all symmetry-allowed force constants up to third nearest neighbors identifies four additional phases as topological via Wilson-loop calculations, for a total of nine topological phases. Ab initio phonon dispersions for 2D Si, Ge, P, As and Sb all fall in the trivial phase. Monte Carlo exploration of the 3NN force-constant space is used to argue that topological phases are physically difficult to realize in this structure.

Significance. If the central claims hold, the work supplies a model-independent TQC classification of five phases together with an explicit enumeration of four more phases inside a 3NN model, plus concrete ab initio results for five real materials. The Monte Carlo sampling quantifies the relative volume of topological versus trivial regions in the allowed parameter space, which is a useful diagnostic for material design. These elements are strengths when the truncation assumptions are justified.

major comments (1)
  1. [Abstract and Monte Carlo section] The headline conclusion that topological phases are physically difficult to realize (abstract; Monte Carlo analysis section) is load-bearing on the claim that the Monte Carlo sampling within the symmetry-allowed 3NN force-constant space is representative. No explicit test is reported that extending the model to 4NN or longer-range interactions (or adding anharmonic terms) leaves the accessible topological volume negligible or that the real-material points remain inside the trivial region.
minor comments (1)
  1. [Analytical model section] Clarify the precise definition and normalization of the force-constant parameters used in the Monte Carlo sampling (e.g., how the symmetry-allowed basis is constructed and scaled).

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comment. We respond point-by-point below.

read point-by-point responses
  1. Referee: [Abstract and Monte Carlo section] The headline conclusion that topological phases are physically difficult to realize (abstract; Monte Carlo analysis section) is load-bearing on the claim that the Monte Carlo sampling within the symmetry-allowed 3NN force-constant space is representative. No explicit test is reported that extending the model to 4NN or longer-range interactions (or adding anharmonic terms) leaves the accessible topological volume negligible or that the real-material points remain inside the trivial region.

    Authors: We agree that an explicit extension to 4NN (or beyond) would provide a stronger test of representativeness. The five phases identified as necessarily non-trivial by topological quantum chemistry are independent of interaction range. The four additional phases were found only within the 3NN model via Wilson-loop calculations; the Monte Carlo sampling was performed exclusively in that 3NN parameter space. Real-material ab initio calculations already incorporate all interaction ranges present in the DFT supercell and place Si, Ge, P, As and Sb firmly in the trivial phase. We did not extend the Monte Carlo to 4NN because the number of independent symmetry-allowed force constants increases substantially, rendering exhaustive sampling impractical. Anharmonic terms lie outside the harmonic phonon framework of the work. We will revise the Monte Carlo section and abstract to explicitly acknowledge the 3NN truncation and to qualify the claim that topological phases are physically difficult to realize as applying within the 3NN model and consistent with the ab initio results. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation uses external TQC and independent MC sampling over model parameters

full rationale

The paper applies topological quantum chemistry (an external framework) to identify five necessarily non-trivial phases from the crystal structure alone, computes ab initio force constants for real materials (Si, Ge, P, As, Sb), builds an analytical dynamical matrix with all symmetry-allowed terms up to third-nearest neighbors, evaluates Wilson loops inside that finite-dimensional model to flag four additional topological phases, and then Monte Carlo samples the symmetry-allowed coefficient space to map volumes. None of these steps reduce by construction to the target claim via self-definition, fitted-input renaming, or a self-citation chain whose load-bearing premise is unverified; the MC exploration is performed directly on the model's parameters rather than being tuned to reproduce the final conclusion. The 3NN truncation is an explicit modeling choice whose sufficiency is debatable but does not create definitional circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the applicability of topological quantum chemistry to phonon spectra and on the completeness of the third-nearest-neighbor force-constant model; no new particles or forces are postulated.

axioms (2)
  • domain assumption Topological quantum chemistry correctly classifies phonon band topology from the crystal structure and symmetry
    Invoked to assert that five phases are necessarily non-trivial
  • domain assumption The space of symmetry-allowed force constants up to third nearest neighbors is sufficient to capture all relevant topology
    Used to construct the analytical model whose Wilson loops identify four additional topological phases

pith-pipeline@v0.9.0 · 5715 in / 1592 out tokens · 22619 ms · 2026-05-24T10:21:40.164272+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

80 extracted references · 80 canonical work pages · 2 internal anchors

  1. [1]

    Topological phonon analysis of the 2D buckled honeycomb lattice: an application to real materials

    applying TQC to 3D materials. In that work they point out the almost absence of fragile cumulative topol- ogy for phonons. Inspired by these results, we deepen in the phonon study by performing Monte Carlo calcula- tions and analysing the force constants in the 2D buckled honeycomb lattice. The TQC analysis relies on detecting an obstruction to a localize...

  2. [2]

    calculations were done using the Quantum ESPRESSO package [46, 47]. We parametrize the exchange-correlation functional assuming the Perdew- Burke-Ernzerhof [66] parametrization and model the electron-ion interaction with projector augmented wave pseudopotentials [67, 68] including four electrons in the valence for Si and Ge, and three electrons for P, As,...

  3. [3]

    (B2) the matrix Φ 11(0) must be in- variant under any operation belonging to the site sym- metry group C3v of atom 1, Φ11(0) =V (S)Φ11(0)V (S)t ∀S∈C3v

    On-site couplings According to Eq. (B2) the matrix Φ 11(0) must be in- variant under any operation belonging to the site sym- metry group C3v of atom 1, Φ11(0) =V (S)Φ11(0)V (S)t ∀S∈C3v. (B3) Applying this equation with S equal to C+ 3 and σd1 (see Fig. (7) for notation), which together generate the group C3v, shows that the on-site matrix must take the f...

  4. [4]

    First nearest neighbors Here we consider the couplings between atoms 1 and 2 in Fig. 1. We can obtain a first constraint on Φ 11(0) by using Eq. (B2) with S = σd1, which leaves atoms 1 and 2 invariant Φ12(0) =V (σd1)Φ12(0)V (σd1)t. (B6) A second constraint is obtained by combining Eq. (B2) for S =I with Eq. (6) Φ12(0) =V (I)Φ21(0)V (I)t = Φ21(0) = Φ12(0)t,...

  5. [5]

    (B2) with S =σd1, that exchanges two sec- ond neighbors, followed by Eq

    Second nearest neighbors Using Eq. (B2) with S =σd1, that exchanges two sec- ond neighbors, followed by Eq. (6), gives σd1Φ11(a2− a1)σ−1 d1 = Φ11(a1− a2) = Φt 11(a2− a1), (B10) which directly leads to: Φ11(a2− a1) =   a2 +b2 −d2 −g2 d2 a2−b2 −h2 g2 −h2 e2  . (B11) Then the remaing second nearest neighbor matrices can be obtained using Φ 22(a1− a2) = Φ...

  6. [6]

    2, the geometry of third nearest neigh- bors is closely related to the one for first neighbors, with σd1 and the inversion I playing analogous roles here

    Third nearest neighbors As seen in Fig. 2, the geometry of third nearest neigh- bors is closely related to the one for first neighbors, with σd1 and the inversion I playing analogous roles here. In- stead of Eqs. (B6) and (B7) we have now Φ12(a1 + a2) =V (σd1)Φ12(a1 + a2)V (σd1)t (B13) and Φ12(a1 + a2) = Φ21(−a1− a2) = Φt 12(a1 + a2), (B14) and these two c...

  7. [7]

    The vector representation is re- ducible, and according to the BCS V =A1(z) +E(x,y )

    Symmetry adapted modes The mechanical representation is induced from the vec- tor representation of C3v, which is the site-symmetry group for the WP 2 b. The vector representation is re- ducible, and according to the BCS V =A1(z) +E(x,y ). (C1) As a consequence, we may compute separately the off- plane modes, induced from A1 and involving atomic dis- place...

  8. [8]

    Eigenvalues of the model dynamical matrix D(k) Changing to the symmetry-adapted basis turns the dy- namical matrix into a block-diagonal form, where the di- mension of each block equals the multiplicity of the corre- sponding irrep. Thus the change of basis yields the eigen- values for all the multiplicity one irreps, while to compute the frequencies asso...

  9. [9]

    4 show strong similarities within the families of pentavalent (P,Sb,As) and tetravalent (Si,Ge) materials

    Mapping DFPT data to the model Even at first glance the DFPT phonon bands in Figs. 4 show strong similarities within the families of pentavalent (P,Sb,As) and tetravalent (Si,Ge) materials. Moreover, one can double check this assumption by comparing the force constants extracted by our DFPT calculations. As shown in Fig 19, the properly rescaled DFPT force...

  10. [10]

    Decay ratios of force constants The decay ratios of force constants used in section IV and Table IV are defined by taking the mean absolute values of the nonzero elements of the matrices of force constants in Fig. 2. Concretely, if we define m0 = 1 3 (2|a0| +|e0|) m1 = 1 5 (|a1 +b1| +|a1−b1| +|e1| + 2|h1|) m2 = 1 9 (|a2 +b2| +|a2−b2| +|e2| + 2|d2| + 2|g2| +...

  11. [11]

    Zhang, P

    T. Zhang, P. Cheng, X. Chen, J.-F. Jia, X. Ma, K. He, L. Wang, H. Zhang, X. Dai, Z. Fang, X. Xie, and Q.-K. Xue, Phys. Rev. Lett. 103, 266803 (2009)

  12. [12]

    B. I. Halperin, Phys. Rev. B 25, 2185 (1982)

  13. [13]

    K¨ onig, S

    M. K¨ onig, S. Wiedmann, C. Br¨ une, A. Roth, H. Buh- mann, L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 318 (2007)

  14. [14]

    Gu and X.-G

    Z.-C. Gu and X.-G. Wen, 80, 155131 (2009)

  15. [15]

    Pollmann, A

    F. Pollmann, A. M. Turner, E. Berg, and M. Oshikawa, Phys. Rev. B 81, 064439 (2010)

  16. [16]

    Chen, Z.-C

    X. Chen, Z.-C. Gu, and X.-G. Wen, Phys. Rev. B 83, 035107 (2011)

  17. [17]

    Fu, Phys

    L. Fu, Phys. Rev. Lett. 106, 106802 (2011)

  18. [18]

    Slager, A

    R.-J. Slager, A. Mesaros, V. Juriˇ ci´ c, and J. Zaanen, Nature Physics 9, 98 (2013)

  19. [19]

    Ando and L

    Y. Ando and L. Fu, Annual Review of Condensed Matter Physics 6, 361 (2015), https://doi.org/10.1146/annurev- conmatphys-031214-014501

  20. [20]

    Z. Wang, A. Alexandradinata, R. J. Cava, and B. A. Bernevig, Nature 532, 189 (2016)

  21. [21]

    W. A. Benalcazar, B. A. Bernevig, and T. L. Hughes, Science 357, 61 (2017)

  22. [22]

    W. A. Benalcazar, B. A. Bernevig, and T. L. Hughes, Phys. Rev. B 96, 245115 (2017)

  23. [23]

    Z. Song, Z. Fang, and C. Fang, Phys. Rev. Lett. 119, 246402 (2017)

  24. [24]

    Langbehn, Y

    J. Langbehn, Y. Peng, L. Trifunovic, F. von Oppen, and P. W. Brouwer, Phys. Rev. Lett. 119, 246401 (2017)

  25. [25]

    B. J. Wieder, B. Bradlyn, Z. Wang, J. Cano, Y. Kim, H.-S. D. Kim, A. M. Rappe, C. L. Kane, and B. A. Bernevig, Science 361, 246 (2018)

  26. [26]

    Schindler, A

    F. Schindler, A. M. Cook, M. G. Vergniory, Z. Wang, S. S. P. Parkin, B. A. Bernevig, and T. Neupert, Science Advances 4 (2018), 10.1126/sciadv.aat0346

  27. [27]

    Khalaf, Phys

    E. Khalaf, Phys. Rev. B 97, 205136 (2018)

  28. [28]

    Fu and C

    L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007)

  29. [29]

    T. L. Hughes, E. Prodan, and B. A. Bernevig, Phys. Rev. B 83, 245132 (2011)

  30. [30]

    A. M. Turner, Y. Zhang, R. S. K. Mong, and A. Vish- wanath, Phys. Rev. B 85, 165120 (2012)

  31. [31]

    C. Fang, M. J. Gilbert, and B. A. Bernevig, Phys. Rev. B 86, 115112 (2012)

  32. [32]

    H. C. Po, A. Vishwanath, and H. Watanabe, Nature Communications 8, 50 (2017)

  33. [33]

    Khalaf, H

    E. Khalaf, H. C. Po, A. Vishwanath, and H. Watanabe, Phys. Rev. X 8, 031070 (2018)

  34. [34]

    Kruthoff, J

    J. Kruthoff, J. de Boer, J. van Wezel, C. L. Kane, and R.-J. Slager, Phys. Rev. X 7, 041069 (2017)

  35. [35]

    Bradlyn, L

    B. Bradlyn, L. Elcoro, J. Cano, M. G. Vergniory, Z. Wang, C. Felser, M. I. Aroyo, and B. A. Bernevig, Nature 547, 298 (2017)

  36. [36]

    M. G. Vergniory, L. Elcoro, C. Felser, N. Regnault, B. A. Bernevig, and Z. Wang, Nature 566, 480 (2019)

  37. [37]

    M. G. Vergniory, B. J. Wieder, L. Elcoro, S. S. P. Parkin, C. Felser, B. A. Bernevig, and N. Regnault, Science 376 (2022)

  38. [38]

    J. Cano, B. Bradlyn, Z. Wang, L. Elcoro, M. G. Vergniory, C. Felser, M. I. Aroyo, and B. A. Bernevig, Phys. Rev. Lett. 120, 266401 (2018)

  39. [39]

    Bradlyn, Z

    B. Bradlyn, Z. Wang, J. Cano, and B. A. Bernevig, Phys. Rev. B 99, 045140 (2019)

  40. [40]

    H. C. Po, H. Watanabe, and A. Vishwanath, Phys. Rev. Lett. 121, 126402 (2018)

  41. [41]

    S¨ usstrunk and S

    R. S¨ usstrunk and S. D. Huber, Science349, 47 (2015)

  42. [42]

    S¨ usstrunk and S

    R. S¨ usstrunk and S. D. Huber, Proceedings of the Na- tional Academy of Sciences 113, E4767 (2016)

  43. [43]

    S. D. Huber, Nature Physics 12, 621 (2016)

  44. [44]

    C. L. Kane and T. C. Lubensky, Nature Physics 10, 39 (2014)

  45. [45]

    Stenull, C

    O. Stenull, C. L. Kane, and T. C. Lubensky, Phys. Rev. Lett. 117, 068001 (2016)

  46. [46]

    J. Mei, Y. Wu, C. T. Chan, and Z.-Q. Zhang, Phys. Rev. B 86, 035141 (2012)

  47. [47]

    Zhang, Z

    T. Zhang, Z. Song, A. Alexandradinata, H. Weng, C. Fang, L. Lu, and Z. Fang, Phys. Rev. Lett. 120, 016401 (2018)

  48. [48]

    J. Li, Q. Xie, S. Ullah, R. Li, H. Ma, D. Li, Y. Li, and X.-Q. Chen, Phys. Rev. B 97, 054305 (2018)

  49. [49]

    Esmann, F

    M. Esmann, F. R. Lamberti, P. Senellart, I. Favero, O. Krebs, L. Lanco, C. Gomez Carbonell, A. Lemaˆ ıtre, 15 and N. D. Lanzillotti-Kimura, Phys. Rev. B 97, 155422 (2018)

  50. [50]

    Singh, Q

    S. Singh, Q. Wu, C. Yue, A. H. Romero, and A. A. Soluyanov, Phys. Rev. Materials 2, 114204 (2018)

  51. [51]

    B. Peng, Y. Hu, S. Murakami, T. Zhang, and B. Monser- rat, Science Advances 6 (2020), 10.1126/sciadv.abd1618

  52. [52]

    J. Li, J. Liu, S. A. Baronett, M. Liu, L. Wang, R. Li, Y. Chen, D. Li, Q. Zhu, and X.-Q. Chen, Nature Com- munications 12, 1204 (2021), arXiv:2006.00705 [cond- mat.mtrl-sci]

  53. [53]

    Y. Xu, M. G. Vergniory, D.-S. Ma, J. L. Ma˜ nes, Z.-D. Song, B. A. Bernevig, N. Regnault, and L. Elcoro, arXiv e-prints , arXiv:2211.11776 (2022)

  54. [54]

    J. L. Ma˜ nes, Phys. Rev. B102, 024307 (2020)

  55. [55]

    Baroni, S

    S. Baroni, S. de Gironcoli, A. Dal Corso, and P. Gian- nozzi, Rev. Mod. Phys. 73, 515 (2001)

  56. [56]

    Giannozzi, S

    P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- cioni, I. Dabo, et al., Journal of physics: Condensed mat- ter 21, 395502 (2009)

  57. [57]

    Giannozzi, O

    P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. B. Nardelli, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, M. Cococcioni, N. Colonna, I. Carnimeo, A. D. Corso, S. de Gironcoli, P. Delugas, R. A. D. Jr, A. Ferretti, A. Floris, G. Fratesi, G. Fugallo, R. Gebauer, U. Gerstmann, F. Giustino, T. Gorni, J. Jia, M. Kawa- mura, H.-Y. Ko, A. Kokalj, E. K¨ u¸...

  58. [58]

    Zak, Phys

    J. Zak, Phys. Rev. Lett. 45, 1025 (1980)

  59. [59]

    Zak, Phys

    J. Zak, Phys. Rev. B 23, 2824 (1981)

  60. [60]

    J. Cano, B. Bradlyn, Z. Wang, L. Elcoro, M. G. Vergniory, C. Felser, M. I. Aroyo, and B. A. Bernevig, Phys. Rev. B 97, 035139 (2018)

  61. [61]

    J. Cano, L. Elcoro, M. I. Aroyo, B. A. Bernevig, and B. Bradlyn, Phys. Rev. B 105, 125115 (2022)

  62. [62]

    Michel and J

    L. Michel and J. Zak, Physics Reports 341, 377 (2001), symmetry, invariants, topology

  63. [63]

    M. B. Walker and J. Zak, Phys. Rev. Lett. 74, 3824 (1995)

  64. [64]

    Kohn, Phys

    W. Kohn, Phys. Rev. B 7, 2285 (1973)

  65. [65]

    Giustino and A

    F. Giustino and A. Pasquarello, Phys. Rev. Lett. 96, 216403 (2006)

  66. [66]

    M. I. Aroyo, J. M. Perez-Mato, C. Capillas, E. Kroumova, S. Ivantchev, G. Madariaga, A. Kirov, and H. Won- dratschek, 221, 15 ()

  67. [67]

    M. I. Aroyo, A. Kirov, C. Capillas, J. M. Perez-Mato, and H. Wondratschek, 62, 115 ()

  68. [68]

    Bilbao crystallographic server,

    “Bilbao crystallographic server,” http://www.cryst.ehu.es

  69. [69]

    Wilczek and A

    F. Wilczek and A. Zee, Phys. Rev. Lett. 52, 2111 (1984)

  70. [70]

    A. A. Soluyanov and D. Vanderbilt, Phys. Rev. B 83, 235401 (2011)

  71. [71]

    Alexandradinata, Z

    A. Alexandradinata, Z. Wang, and B. A. Bernevig, Phys. Rev. X 6, 021008 (2016)

  72. [72]

    C. J. Bradley and A. P. Cracknell, The mathematical theory of symmetry in solids (Clarendon Press, 1972)

  73. [73]

    Python tight binding (pythtb),

    “Python tight binding (pythtb),”

  74. [74]

    Taherinejad, K

    M. Taherinejad, K. F. Garrity, and D. Vanderbilt, Phys. Rev. B 89, 115102 (2014)

  75. [75]

    Z. Peng, X. Chen, Y. Fan, D. J. Srolovitz, and D. Lei, Light: Science & Applications 9, 190 (2020)

  76. [76]

    J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Phys. Rev. Lett. 100, 136406 (2008)

  77. [77]

    Kresse and J

    G. Kresse and J. Furthm¨ uller, Phys. Rev. B 54, 11169 (1996)

  78. [78]

    Kresse and J

    G. Kresse and J. Furthm¨ uller, Computational Materials Science 6, 15 (1996)

  79. [79]

    Methfessel and A

    M. Methfessel and A. T. Paxton, Phys. Rev. B 40, 3616 (1989)

  80. [80]

    A. A. Maradudin and S. H. Vosko, Rev. Mod. Phys. 40, 1 (1968)