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
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.
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
- 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
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.
Referee Report
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)
- [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)
- [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
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
-
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
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
axioms (2)
- domain assumption Topological quantum chemistry correctly classifies phonon band topology from the crystal structure and symmetry
- domain assumption The space of symmetry-allowed force constants up to third nearest neighbors is sufficient to capture all relevant topology
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking (D=3 forcing) unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
eleven distinct phases... five of which necessarily have non-trivial topology according to topological quantum chemistry. Another four... Wilson loops in an analytical model that includes all the symmetry allowed force constants up to third nearest neighbors
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel (J-cost uniqueness) unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Monte Carlo calculations... force constants decay faster than according to a (100:35:4) ratio... percentage of space occupied by each phase
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
-
[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...
work page internal anchor Pith review Pith/arXiv arXiv 2022
-
[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]
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]
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]
(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]
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]
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]
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]
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]
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]
-
[12]
B. I. Halperin, Phys. Rev. B 25, 2185 (1982)
work page 1982
-
[13]
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)
work page 2007
- [14]
-
[15]
F. Pollmann, A. M. Turner, E. Berg, and M. Oshikawa, Phys. Rev. B 81, 064439 (2010)
work page 2010
- [16]
- [17]
- [18]
-
[19]
Y. Ando and L. Fu, Annual Review of Condensed Matter Physics 6, 361 (2015), https://doi.org/10.1146/annurev- conmatphys-031214-014501
-
[20]
Z. Wang, A. Alexandradinata, R. J. Cava, and B. A. Bernevig, Nature 532, 189 (2016)
work page 2016
-
[21]
W. A. Benalcazar, B. A. Bernevig, and T. L. Hughes, Science 357, 61 (2017)
work page 2017
-
[22]
W. A. Benalcazar, B. A. Bernevig, and T. L. Hughes, Phys. Rev. B 96, 245115 (2017)
work page 2017
-
[23]
Z. Song, Z. Fang, and C. Fang, Phys. Rev. Lett. 119, 246402 (2017)
work page 2017
-
[24]
J. Langbehn, Y. Peng, L. Trifunovic, F. von Oppen, and P. W. Brouwer, Phys. Rev. Lett. 119, 246401 (2017)
work page 2017
-
[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)
work page 2018
-
[26]
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]
- [28]
-
[29]
T. L. Hughes, E. Prodan, and B. A. Bernevig, Phys. Rev. B 83, 245132 (2011)
work page 2011
-
[30]
A. M. Turner, Y. Zhang, R. S. K. Mong, and A. Vish- wanath, Phys. Rev. B 85, 165120 (2012)
work page 2012
-
[31]
C. Fang, M. J. Gilbert, and B. A. Bernevig, Phys. Rev. B 86, 115112 (2012)
work page 2012
-
[32]
H. C. Po, A. Vishwanath, and H. Watanabe, Nature Communications 8, 50 (2017)
work page 2017
- [33]
-
[34]
J. Kruthoff, J. de Boer, J. van Wezel, C. L. Kane, and R.-J. Slager, Phys. Rev. X 7, 041069 (2017)
work page 2017
-
[35]
B. Bradlyn, L. Elcoro, J. Cano, M. G. Vergniory, Z. Wang, C. Felser, M. I. Aroyo, and B. A. Bernevig, Nature 547, 298 (2017)
work page 2017
-
[36]
M. G. Vergniory, L. Elcoro, C. Felser, N. Regnault, B. A. Bernevig, and Z. Wang, Nature 566, 480 (2019)
work page 2019
-
[37]
M. G. Vergniory, B. J. Wieder, L. Elcoro, S. S. P. Parkin, C. Felser, B. A. Bernevig, and N. Regnault, Science 376 (2022)
work page 2022
-
[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)
work page 2018
-
[39]
B. Bradlyn, Z. Wang, J. Cano, and B. A. Bernevig, Phys. Rev. B 99, 045140 (2019)
work page 2019
-
[40]
H. C. Po, H. Watanabe, and A. Vishwanath, Phys. Rev. Lett. 121, 126402 (2018)
work page 2018
- [41]
-
[42]
R. S¨ usstrunk and S. D. Huber, Proceedings of the Na- tional Academy of Sciences 113, E4767 (2016)
work page 2016
-
[43]
S. D. Huber, Nature Physics 12, 621 (2016)
work page 2016
-
[44]
C. L. Kane and T. C. Lubensky, Nature Physics 10, 39 (2014)
work page 2014
-
[45]
O. Stenull, C. L. Kane, and T. C. Lubensky, Phys. Rev. Lett. 117, 068001 (2016)
work page 2016
-
[46]
J. Mei, Y. Wu, C. T. Chan, and Z.-Q. Zhang, Phys. Rev. B 86, 035141 (2012)
work page 2012
- [47]
-
[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)
work page 2018
- [49]
- [50]
-
[51]
B. Peng, Y. Hu, S. Murakami, T. Zhang, and B. Monser- rat, Science Advances 6 (2020), 10.1126/sciadv.abd1618
- [52]
-
[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)
work page internal anchor Pith review Pith/arXiv arXiv 2022
-
[54]
J. L. Ma˜ nes, Phys. Rev. B102, 024307 (2020)
work page 2020
- [55]
-
[56]
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)
work page 2009
-
[57]
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¸...
work page 2017
- [58]
- [59]
-
[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)
work page 2018
-
[61]
J. Cano, L. Elcoro, M. I. Aroyo, B. A. Bernevig, and B. Bradlyn, Phys. Rev. B 105, 125115 (2022)
work page 2022
-
[62]
L. Michel and J. Zak, Physics Reports 341, 377 (2001), symmetry, invariants, topology
work page 2001
-
[63]
M. B. Walker and J. Zak, Phys. Rev. Lett. 74, 3824 (1995)
work page 1995
- [64]
- [65]
-
[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]
M. I. Aroyo, A. Kirov, C. Capillas, J. M. Perez-Mato, and H. Wondratschek, 62, 115 ()
- [68]
- [69]
-
[70]
A. A. Soluyanov and D. Vanderbilt, Phys. Rev. B 83, 235401 (2011)
work page 2011
-
[71]
A. Alexandradinata, Z. Wang, and B. A. Bernevig, Phys. Rev. X 6, 021008 (2016)
work page 2016
-
[72]
C. J. Bradley and A. P. Cracknell, The mathematical theory of symmetry in solids (Clarendon Press, 1972)
work page 1972
- [73]
-
[74]
M. Taherinejad, K. F. Garrity, and D. Vanderbilt, Phys. Rev. B 89, 115102 (2014)
work page 2014
-
[75]
Z. Peng, X. Chen, Y. Fan, D. J. Srolovitz, and D. Lei, Light: Science & Applications 9, 190 (2020)
work page 2020
-
[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)
work page 2008
- [77]
-
[78]
G. Kresse and J. Furthm¨ uller, Computational Materials Science 6, 15 (1996)
work page 1996
- [79]
-
[80]
A. A. Maradudin and S. H. Vosko, Rev. Mod. Phys. 40, 1 (1968)
work page 1968
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.