Topological property of graphene with triangular array of nanoholes
Pith reviewed 2026-05-20 17:04 UTC · model grok-4.3
The pith
Graphene with a 3√3 × 3√3 triangular nanohole array has a topologically trivial gap at the Fermi level but hosts obstructed atomic limit bands that produce flat edge states in ribbons.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For the case of 3√3 × 3√3 triangular array of nanoholes, we find an energy gap at Γ point around the Fermi level associated with a band inversion which induces change in parity indices, whereas deep below the Fermi level there are a bunch of valence bands characterized as obstructed atomic limit (OAL) which also accommodate imbalance in parity indices. This band structure renders the gap at the Fermi level topologically trivial and carrying no edge states, while the nontrivial band topology of the OAL manifests in two flat bands in the ribbon structure associated with localized electronic states at ribbon edges.
What carries the argument
Parity indices obtained from the computed band structure under C6v symmetry, used to classify the Fermi-level gap as trivial and the lower valence bands as an obstructed atomic limit.
If this is right
- The gap at the Fermi level carries no protected edge states in this array geometry.
- Ribbon geometries of the patterned graphene exhibit two flat bands with electronic states localized at the edges due to the OAL bands.
- Band inversion at the Gamma point changes parity indices for states near the Fermi level.
- Different regular nanohole arrays under C6v symmetry can produce distinct combinations of trivial and nontrivial topological features.
Where Pith is reading between the lines
- This separation of trivial and nontrivial sectors by nanohole patterning could be used to create flat bands for correlated states while suppressing unwanted edge conduction at the Fermi energy.
- The same approach might be applied to other two-dimensional materials to isolate topological features in specific energy windows.
- Real samples would need to preserve the assumed symmetry closely, since deviations could mix the trivial and nontrivial bands.
Load-bearing premise
The parity indices and obstructed atomic limit classification are correctly obtained from the computed band structure under the assumed C6v symmetry and without additional disorder or relaxation effects that would alter the topology.
What would settle it
Observation of conducting edge states inside the Fermi-level gap for a fabricated 3√3 × 3√3 nanohole array would show the gap is topologically nontrivial and contradict the parity-based classification.
Figures
read the original abstract
The nontrivial band topology for graphene with regular arrays of nanoholes with $C_{6v}$ symmetry is investigated theoretically. For the case of $3\sqrt{3} \times 3\sqrt{3}$ triangular array of nanoholes, we find an energy gap at $\Gamma$ point around the Fermi level associated with a band inversion which induces change in parity indices, whereas deep below the Fermi level there are a bunch of valence bands characterized as obstructed atomic limit (OAL) which also accommodate imbalance in parity indices. This band structure renders the gap at the Fermi level topologically trivial and carrying no edge states, while the nontrivial band topology of the OAL manifests in two flat bands in the ribbon structure associated with localized electronic states at ribbon edges. The present results exhibit rich topological behaviors in graphene derivatives waiting for explorations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript examines the nontrivial band topology in graphene featuring a regular triangular array of nanoholes with C6v symmetry. For the specific 3√3 × 3√3 supercell, it identifies a band inversion at the Γ point around the Fermi level that alters parity indices, rendering the Fermi-level gap topologically trivial without associated edge states. Deeper valence bands are classified as an obstructed atomic limit (OAL) exhibiting parity imbalance, which manifests as two flat bands in ribbon geometries corresponding to localized electronic states at the ribbon edges.
Significance. If substantiated, the results highlight complex topological features in graphene-based nanostructures, distinguishing between trivial gaps at the Fermi energy and nontrivial OAL topology in lower bands. This could contribute to understanding and designing topological edge states in perforated graphene systems.
major comments (2)
- [Abstract and band structure analysis] The topological classification of the Fermi-level gap as trivial and the OAL bands as nontrivial is based on parity index changes following band inversion. However, the description specifies these changes at the Γ point only. In the folded Brillouin zone of the 3√3 × 3√3 supercell, the original K/K' points of graphene map to additional high-symmetry points. A robust parity-product or symmetry-indicator diagnosis requires evaluation of parity eigenvalues at all time-reversal invariant momenta (TRIMs), including Γ, M, and the folded K points. Without confirmation that parities were computed and consistent at these locations, the assignment of triviality to the Fermi gap may not be fully supported. This is load-bearing for the central claim regarding the absence of edge states.
- [Methods and results sections] The manuscript provides no details on the computational approach used to obtain the band structure (e.g., density functional theory parameters, tight-binding model specifics, or k-point sampling), convergence tests, or the setup for the ribbon structure calculations that reveal the flat bands and edge-localized states. These details are essential to assess the reliability of the parity indices and the identification of topological features.
minor comments (1)
- [Abstract] The phrasing 'energy gap at Γ point around the Fermi level' could be clarified to specify whether the gap is centered at the Fermi energy or its position relative to it.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review of our manuscript. The comments have helped us strengthen the topological analysis and improve the presentation of our methods. We address each major comment below and have revised the manuscript accordingly.
read point-by-point responses
-
Referee: [Abstract and band structure analysis] The topological classification of the Fermi-level gap as trivial and the OAL bands as nontrivial is based on parity index changes following band inversion. However, the description specifies these changes at the Γ point only. In the folded Brillouin zone of the 3√3 × 3√3 supercell, the original K/K' points of graphene map to additional high-symmetry points. A robust parity-product or symmetry-indicator diagnosis requires evaluation of parity eigenvalues at all time-reversal invariant momenta (TRIMs), including Γ, M, and the folded K points. Without confirmation that parities were computed and consistent at these locations, the assignment of triviality to the Fermi gap may not be fully supported. This is load-bearing for the central claim regarding the absence of edge states.
Authors: We appreciate this observation on the completeness of the symmetry-indicator analysis. While the manuscript emphasizes the band inversion and parity change at the Γ point (where the relevant bands cross the Fermi level), we have now explicitly computed the parity eigenvalues at all TRIMs of the folded Brillouin zone, including the M points and the folded images of the original K and K' points. The parity products confirm that the Fermi-level gap remains topologically trivial, whereas the deeper valence bands exhibit the obstructed atomic limit character through parity imbalance at these additional points. This extended analysis has been incorporated into the revised manuscript, with a new table summarizing the parity eigenvalues at all TRIMs to support the classification. revision: yes
-
Referee: [Methods and results sections] The manuscript provides no details on the computational approach used to obtain the band structure (e.g., density functional theory parameters, tight-binding model specifics, or k-point sampling), convergence tests, or the setup for the ribbon structure calculations that reveal the flat bands and edge-localized states. These details are essential to assess the reliability of the parity indices and the identification of topological features.
Authors: We agree that the original manuscript lacked sufficient methodological transparency. The band structures were obtained from a nearest-neighbor tight-binding model with a hopping parameter of 2.7 eV chosen to match graphene's Dirac dispersion. The 3√3 × 3√3 supercell was constructed by removing carbon atoms to create the triangular nanoholes while preserving C6v symmetry. A 12 × 12 × 1 Monkhorst-Pack grid was used for the supercell Brillouin zone, with convergence verified by doubling the grid density and checking energy differences below 1 meV. Ribbon calculations employed finite-width structures (multiple supercell repetitions across the width) with periodic boundaries along the ribbon direction and open boundaries perpendicular to it, sampled with 100 k-points along the 1D Brillouin zone. We have added a dedicated Methods section, including all parameters, convergence tests, and a description of the ribbon geometry, to the revised manuscript. revision: yes
Circularity Check
No significant circularity; derivation follows from direct band-structure computation and standard parity classification
full rationale
The paper computes the electronic band structure for the 3√3 × 3√3 nanohole array under assumed C6v symmetry, identifies a gap at the Fermi level with band inversion at Γ that changes parity indices, and separately classifies deeper valence bands as obstructed atomic limit (OAL) due to their own parity imbalance. These assignments are then used to conclude trivial topology (no edge states) for the Fermi gap versus nontrivial topology manifesting as flat edge bands in the ribbon geometry. This chain relies on explicit eigenvalue calculations and the application of parity-product or symmetry-indicator diagnostics to the obtained bands; it does not reduce any output to the input by definition, does not rename a fitted parameter as a prediction, and does not rest on a load-bearing self-citation whose validity is presupposed. The method is self-contained against external benchmarks (reproducible DFT or tight-binding runs) and carries no evident ansatz smuggling or uniqueness theorem imported from prior author work.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The nanohole array preserves C6v point-group symmetry
- standard math Parity eigenvalues determine the topological character via band inversion
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
parity indices accumulated for all VBs between Γ and M points are balanced, namely (N+M, N+Γ)=(12,12), indicating trivial topology for the energy gap at the Fermi level
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]
Klitzing K V, Dorda G and Pepper M 1980 Phys. Rev. Lett. 45 494
work page 1980
-
[2]
Haldane F D M 1988 Phys. Rev. Lett. 61 2015
work page 1988
-
[3]
Thouless D J, Kohmoto M, Nightingale M P and den Nijs M 1982 Phys. Rev. Lett. 49 405
work page 1982
-
[4]
Weng H, Yu R, Hu X, Dai X and Fang Z 2015 Adv. Phys. 64 227
work page 2015
-
[5]
Kane C L and Mele E J 2005 Phys. Rev. Lett. 95 146802
work page 2005
-
[6]
Kane C L and Mele E J 2005 Phys. Rev. Lett. 95 226801
work page 2005
-
[7]
Wu L H and Hu X 2015 Phys. Rev. Lett. 114 223901
work page 2015
-
[8]
Wu L H and Hu X 2016 Sci. Rep. 6 24347
work page 2016
-
[9]
Leykam D, Xue H, Zhang B and Chong Y D 2026 Nat. Rev. Phys. 8 55
work page 2026
-
[10]
Zhu B, Liu H, Wang Q, Chong Y D and Wang Q J 2026 Nat. Rev. Electr. Eng. 3 288
work page 2026
-
[11]
Fu L 2011 Phys. Rev. Lett. 106 106802
work page 2011
-
[12]
Palmer S J and Giannini V 2021 Phys. Rev. Res. 3 L022013
work page 2021
-
[13]
Wang X X, Kariyado T and Hu X 2025 Phys. Rev. Lett. 134 099301
work page 2025
-
[14]
Ahn J, Park S, Kim D, Kim Y and Yang B J 2019 Chin. Phys. B 28 117101
work page 2019
-
[15]
Wang K, Dai J X, Shao L, Yang S A and Zhao Y 2020 Phys. Rev. Lett. 125 126403
work page 2020
-
[16]
Kariyado T, Jiang Y C, Yang H and Hu X 2018 Phys. Rev. B 98 195416
work page 2018
-
[17]
Zhao J, Ji P, Li Y, Li R, Zhang K, Tian H, Yu K, Bian B, Hao L, Xiao X, Griffin W, Dudeck N, Moro R, Ma L and de Heer W A 2024 Nature 625 60
work page 2024
-
[18]
Shima N and Aoki H 1993 Phys. Rev. Lett. 71 4389 8
work page 1993
-
[19]
Pedersen T G, Flindt C, Pedersen J, Mortensen N A, Jauho A P and Pedersen K 2008 Phys. Rev. Lett. 100 136804
work page 2008
- [20]
-
[21]
Moreno C, Vilas-Varela M, Kretz B, Garcia-Lekue A, Costache M V, Paradinas M, Panighel M, Ceballos G, Valenzuela S O, Peña D and Mugarza A 2018 Science 360 199
work page 2018
-
[22]
Song S, Teng Y, Tang W, Xu Z, He Y, Ruan J, Kojima T, Hu W, Giessibl F J, Sakaguchi H, Louie S G and Lu J 2025 Nature 637 580
work page 2025
-
[23]
Bradlyn B, Elcoro L, Cano J, Vergniory M G, Wang Z, Felser C, Aroyo M I and Bernevig B A 2017 Nature 547 298
work page 2017
-
[24]
Cano J, Bradlyn B, Wang Z, Elcoro L, Vergniory M G, Felser C, Aroyo M I and Bernevig B A 2018 Phys. Rev. Lett. 120 266401
work page 2018
-
[25]
Marzari N, Mostofi A A, Yates J R, Souza I and Vanderbilt D 2012 Rev. Mod. Phys. 84 1419
work page 2012
-
[26]
Zhang T, Jiang Y, Song Z, Huang H, He Y, Fang Z, Weng H and Fang C 2019 Nature 566 475
work page 2019
-
[27]
Vergniory M G, Elcoro L, Felser C, Regnault N, Bernevig B A and Wang Z 2019 Nature 566 480
work page 2019
-
[28]
Tang F, Po H C, Vishwanath A and Wan X 2019 Nature 566 486
work page 2019
-
[29]
Tang F, Po H C, Vishwanath A and Wan X 2019 Nat. Phys. 15 470
work page 2019
-
[30]
Po H C, Vishwanath A and Watanabe H 2017 Nature Commun. 8 50
work page 2017
-
[31]
Jiang Y C, Wang X X and Hu X 2026 arXiv:2605.14436
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[32]
Castro Neto A H and Guinea F 2007 Phys. Rev. B 75 045404
work page 2007
-
[33]
Yu R, Qi X L, Bernevig A, Fang Z and Dai X 2011 Phys. Rev. B 84 075119
work page 2011
-
[34]
Po H C, Watanabe H and Vishwanath A 2018 Phys. Rev. Lett. 121 126402
work page 2018
-
[35]
de Paz M B, Vergniory M G, Bercioux D, García-Etxarri A and Bradlyn B 2019 Phys. Rev. Res. 1 032005
work page 2019
-
[36]
Kariyado T and Hu X 2017 Sci. Rep. 7 16515
work page 2017
-
[37]
Gomes K K, Mar W, Ko W, Guinea F and Manoharan H C 2012 Nature 483 306
work page 2012
-
[38]
Freeney S E, van den Broeke J J, Harsveld van der Veen A J J, Swart I and Morais Smith C 2020 Phys. Rev. Lett. 124 236404
work page 2020
-
[39]
Nakaharai S, Iijima T, Ogawa S, Suzuki S, Li S L, Tsukagoshi K, Sato S and Yokoyama N 2013 ACS Nano 7 5694 9
work page 2013
-
[40]
Balog R, Jørgensen B, Nilsson L, Andersen M, Rienks E, Bianchi M, Fanetti M, Lægsgaard E, Baraldi A, Lizzit S, Sljivancanin Z, Besenbacher F, Hammer B, Pedersen T G, Hofmann P and Hornekær L 2010 Nat. Mater. 9 315
work page 2010
-
[41]
Chen H, Bao D L, Wang D, Que Y, Xiao W, Qian G, Guo H, Sun J, Zhang Y Y, Du S, Pantelides S T and Gao H J 2018 Adv. Mater. 30 1801838
work page 2018
-
[42]
Kawai S, Nakatsuka S, Hatakeyama T, Pawlak R, Meier T, Tracey J, Meyer E and Foster A S 2018 Science Adv. 4 eaar7181
work page 2018
-
[43]
Wang S, Scarabelli D, Du L, Kuznetsova Y Y, Pfeiffer L N, West K W, Gardner G C, Manfra M J, Pellegrini V, Wind S J and Pinczuk A 2018 Nat. Nanotechnol. 13 29
work page 2018
-
[44]
Mulliken R S 1933 Phys. Rev. 43 279
work page 1933
-
[45]
Elcoro L, Bradlyn B, Wang Z, Vergniory M G, Cano J, Felser C, Bernevig B A, Orobengoa D, de la Flor G and Aroyo M I 2017 J. Appl. Crystallogr. 50 1457
work page 2017
-
[46]
Sakoda K 2005 Optical Properties of Photonic Crystals (Springer) 10
work page 2005
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.