Defects Potentials for Two-Dimensional Topological Materials

Amit Goft; Eric Akkermans; Nadav Orion; Yuval Abulafia

REVIEW 3 major objections 6 minor 2 cited by

Defects Potentials for Two-Dimensional Topological Materials

T0 review · 3 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper shows that a single vacancy in a honeycomb lattice creates exactly one topological zero mode with analytical index -1, while an adatom creates none, and derives the defect potentials that make this difference.

desk verdict Useful new defect-potential derivation and a testable density signature, but the claimed index -1 for a vacancy is not yet justified in the continuum model because the counted zero mode is non-normalizable and no regularization is supplied. read the letter →

arxiv 2507.01530 v1 pith:3NGINU57 submitted 2025-07-02 cond-mat.mes-hall

Yuval Abulafia , Amit Goft , Nadav Orion , Eric Akkermans This is my paper · ORCID

classification cond-mat.mes-hall

keywords analyticalindextopologicalzeromodesgraphenevacancydefectpotentialschiralsymmetryDiracHamiltonianvalleycouplingGreen'sfunctions

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper aims to show that localized defects in two-dimensional bipartite Dirac materials can be engineered to turn an otherwise non-topological material into a topological one, with topology measured by the analytical index, the difference between the number of zero modes of an operator and of its adjoint. Starting from tight-binding models, it derives first-quantized continuum potentials for two concrete defects in graphene: an adatom, which breaks chiral symmetry and yields no zero modes, and a vacancy, which preserves chiral symmetry and produces exactly one topological zero mode localized on the B sublattice, with analytical index -1. The vacancy zero mode is $\psi_{\mathrm{zm}}(r) = (1/r)(0,0,e^{i\theta},e^{-i\theta})^{T}$, a non-normalizable power-law state; the paper also shows how the same Green's function machinery gives distinct measurable density signatures for the two defects. A sympathetic reader would care because the method turns defect engineering into a design rule: off-diagonal, valley-coupling potential structure induces topology, while diagonal mass-like structure does not.

What carries the argument

The carrying object is the first-quantized Dirac Hamiltonian in the Hilbert space $L^2(\mathbb{R}^2)\otimes S^2_{KK'}\otimes S^2_{AB}$, with all defect potentials written as $4\times4$ matrix operators obtained by expanding tight-binding operators around the two valleys. The vacancy potential (34) is the key object: its off-diagonal blocks couple valleys through $\delta(r-R_0)\hat{L}$ and $\hat{L}^\dagger\delta(r-R_0)$ terms, making the off-diagonal block $\hat{D} = \hat{D}_0 + \hat{D}_{\mathrm{vac}}$ non-Hermitian, which is the necessary condition for a nonzero analytical index. The index itself, $\mathrm{Index}\,\hat{D} = \dim\ker\hat{D} - \dim\ker\hat{D}^\dagger$, counts topological zero modes, and the bulk-edge correspondence (39) ties it to the winding number $\nu$ and the zero-mode number $N_{\mathrm{zm}}$. On the solution side, the ladder relations $\hat{L}|m\rangle = |m-1\rangle$, $\hat{L}^\dagger|m\rangle = |m+1\rangle$ select exactly one surviving angular momentum channel on the B sublattice, giving the $1/r$ mode (48).

What would settle it

A decisive check is a chiral tight-binding calculation on a finite honeycomb lattice with one A-sublattice vacancy: count exact zero-energy states and their sublattice localization. If the number of B-sublattice zero modes is not exactly one, or if the localized tail is not $1/r$, the claim $\mathrm{Index}\,\hat{D} = -1$ fails. Experimentally, an STM map of the intervalley local density around a single vacancy should show one wavefront dislocation and the angular dependence $\cos(\Delta K\cdot r) - \cos(\Delta K\cdot r + 2\theta)$ with no $2k_F r$ Friedel oscillations; two wavefront dislocations would contradict the vacancy potential (34).

Watch

Extended reading notes

Core claim

The central claim is that a vacancy in graphene is topological while an adatom is trivial. In the continuum description built on a four-dimensional pseudospin Hilbert space (valley $\{K,K'\}$ times sublattice $\{A,B\}$), pristine graphene has an off-diagonal chiral Hamiltonian $\hat{H}_0$ with block $\hat{D}_0 = v_F\,\mathrm{diag}(\hat{L}^\dagger,-\hat{L})$, where $\hat{L} = -i\partial_x - \partial_y$. Adding the vacancy potential $\hat{V}_{\mathrm{vac}}$ replaces $\hat{D}_0$ by $\hat{D} = \hat{D}_0 + \hat{D}_{\mathrm{vac}}$, and because $\hat{D}_{\mathrm{vac}}$ contains derivatives acting on $\delta(r-R_0)$, this block is non-Hermitian, $\hat{D} \neq \hat{D}^\dagger$. The paper solves the zero-mode equations using angular momentum ladder relations and finds a single non-normalizable solution on the B sublattice, $\psi_{\mathrm{zm}}(r) = (1/r)(0,0,e^{i\theta},e^{-i\theta})^{T}$, so $\mathrm{Index}\,\hat{D} = -1$ and $|\nu| = N_{\mathrm{zm}} = 1$. The adatom potential is diagonal in the sublattice basis, leaves $\hat{D}$ Hermitian, and gives index zero. The accompanying Green's function calculation distinguishes the two cases in the intervalley local density: the adatom gives $2k_F r$ Friedel oscillations, while the vacancy gives a $1/r^2$ decay with angular pattern $\cos(\Delta K\cdot r) - \cos(\Delta K\cdot r + 2\theta)$ and no Friedel oscillations.

Load-bearing premise

The load-bearing premise is that the index on the infinite plane counts the non-normalizable $1/r$ zero mode, which the paper itself flags as such, as a genuine element of $\ker \hat{D}^\dagger$; if square-integrability were required, both kernels would be empty and the vacancy would have index zero exactly like the adatom.

Editorial extensions

If this is right

A single A-sublattice vacancy in chiral graphene hosts exactly one topological zero-energy mode on the B sublattice, with wavefunction $\frac{1}{r}(0,0,e^{i\theta},e^{-i\theta})^T$ and analytical index $-1$.
An adatom creates no topological zero modes and index $0$; zero modes that appear on finite flakes are supernumerary and not counted by the index.
The vacancy-induced intervalley local density decays as $1/r^2$ with angular dependence $\cos(\Delta K\cdot r) - \cos(\Delta K\cdot r + 2\theta)$ and no Friedel oscillations, cleanly distinguishing a vacancy from an adatom.
Two vacancies on the same sublattice give index $\pm 2$, one A and one B vacancy give index $0$, and an added adatom mass term does not change the vacancy's topological zero-mode count.
Topological textures such as the Kekulé distortion share the vacancy potential's matrix structure, so the same design rule for inducing topology applies.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Beyond the paper: because the index count relies on a non-normalizable mode, a finite-size or regularized version of the calculation could shift the index unless a boundary condition is specified; testing the $1/r$ tail by scaling the system size would settle this.
Beyond the paper: the design rule that off-diagonal intervalley coupling making $\hat{D}$ non-Hermitian induces topology, while diagonal mass terms do not, should transfer to other bipartite Dirac lattices such as silicene, germanene, and brickwall lattices with fermion doubling.
Beyond the paper: the predicted single wavefront dislocation and $1/r^2$ density decay are directly testable in STM experiments on isolated vacancies; a two-dislocation pattern would indicate that the continuum potential needs additional short-distance structure.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Two things you should know. The paper delivers a genuinely useful systematic derivation of continuum defect potentials from tight-binding for graphene vacancies and adatoms, and it produces a concrete, parameter-free prediction for the local density change at a vacancy, Eq. (61), which is a clean signature distinguishing the vacancy from an adatom. The second thing is that the paper's central topological claim—Index D = -1 for a vacancy—is not actually proven in the continuum model as written. The zero mode they count, Eq. (48), is non-normalizable, and the paper says so itself. The index formula (38) is a Fredholm index for elliptic operators, but on the non-compact plane with a delta-potential the operator is not Fredholm in the standard L2 domain. The formal angular-momentum calculation in Eqs. (46)-(47) writes delta(r) r^{-m}, which is not a distribution without extra regularization. So the value -1 is a property of a formal kernel, not of the operator as defined.

The new content is real. The vacancy potential matrix (34) with delta(r) times derivative operators, and the first-order density change (61) with its cos(ΔK·r) - cos(ΔK·r+2θ) angular structure, are not in the literature in this exact form. The Green's function machinery in Appendices A-C is detailed and, at the formal level, internally consistent. The comparison to the adatom (57) is instructive and the qualitative distinction—single wavefront dislocation for the vacancy versus two for the adatom—is an externally testable prediction. The paper also cites the relevant literature on graphene vacancies fairly; the prior zero mode [18,20] is acknowledged.

The soft spot is exactly the index. The authors flag the non-normalizability but then use the mode to fix the index. They suggest replacing the delta with a finite-range potential, but that changes the normalizability question and doesn't automatically give an L2 zero mode with a 1/r tail. The physical vacancy zero mode is real and well-known; the finite-lattice index N_B - N_A = 1 gives the same count. What's missing is an explicit statement that the index is defined through that finite-lattice regularization, or a weighted-space setup. That's a repair, not a fatal flaw, but it is the load-bearing piece of the topological classification story, so it needs to be written down.

Who is this for: condensed matter theorists working on defects in Dirac materials, topological classification beyond translation-invariant systems, and STM signatures of point defects. It deserves a serious referee. I would send it to review and ask for the index regularization to be addressed, and for a clearer separation between the new potential/density results and the known zero-mode fact. I'd probably cite it for the density formula.

Referee Report

3 major / 6 minor

Summary. The paper develops a systematic procedure to derive continuum defect potentials from tight-binding models for bipartite Dirac materials, and applies it to an adatom, a vacancy, and the combination of both in graphene. The central claim is that the vacancy potential produces a non-Hermitian chiral Dirac operator with analytical index Index D = -1, corresponding to exactly one topological zero mode on the B sublattice with wavefunction psi_zm(r) = (1/r)(0,0,e^{iθ},e^{-iθ})^T, while the adatom potential preserves index zero and no topological zero modes. The paper further uses these potentials in a Green's function formalism to compute first- and all-order valley-coupled local density changes, predicting distinct Friedel-oscillation behavior and a single wavefront dislocation for a vacancy.

Significance. If the index assignment can be made mathematically rigorous, the paper provides a parameter-free continuum method for classifying point-defect-induced topology in two-dimensional Dirac materials and generates concrete, falsifiable predictions: one topological zero mode at a vacancy, index-zero adatom, and specific density oscillations with a single wavefront dislocation. The derivations from tight-binding to continuum potentials are systematic and contain no fitted parameters. The Green's function density result (61) is an explicit calculation that can be tested against numerics or scanning tunneling microscopy. These strengths make the manuscript potentially valuable, provided the functional-analytic issues described below are resolved.

major comments (3)

[Section IV.B, Eqs. (46)-(47) and Appendix B] The assignment Index D = -1 counts the non-normalizable mode (48), psi_zm(r) = (1/r)(0,0,e^{iθ},e^{-iθ})^T, as an element of ker D†. On L2(R2) this function is not square-integrable, so with standard Sobolev or L2 domains ker D† is empty and the operator D on the non-compact plane is not Fredholm; the essential spectrum reaches zero because the massless Dirac symbol vanishes at zero momentum. The paper explicitly states that the mode is non-normalizable but does not supply a regularization, a weighted-space formulation, or a limiting procedure that makes Eq. (38) well-defined. Please specify the functional setting or the finite-lattice/large-radius limit in which the index is computed and show that its value is indeed -1 in that setting.
[Section IV.B, Eqs. (46)-(47) and Appendix B] The angular-momentum derivation evaluates expressions such as δ(r) r^{-m} and δ(r) L† r^m, which are not well-defined as distributions because r^{-m} is singular at the origin and the derivative operators act on the delta distribution in an order that requires a specified convention. Likewise, Appendix B, around Eqs. (B1)-(B3), states that one integrates by parts to handle derivatives of δ(r1-r2), but this is not a standard identity for products of derivative operators and delta functions unless the operator ordering is defined. As written, both the zero-mode counting and the Green's function expressions (59)-(61) rest on formal manipulations; please provide a regularization or a distributional definition that justifies these steps.
[Section V.B.2, Eqs. (60)-(61)] The text says that 'the result from (61) remains unaffected by the Fermi energy choice, whether set to zero or any non-zero value due to gating' and then immediately adds 'However, I(r) is influenced by both the Fermi energy and the lower boundary selection.' Since δρ(1)_vac(r) = I(r)[cos(ΔK·r) − cos(ΔK·r + 2θ)], these two statements are in direct tension. Please clarify whether the Fermi-energy independence refers only to the angular structure of the density, and specify the precise conditions under which the prediction (61) is expected to hold.

minor comments (6)

[Section I] There is a missing space in 'thebuildingofdistinguishablequbits' in the first paragraph of the introduction.
[Eq. (58)] The notation H^{(1)}_{0,1}(z) is ambiguous; please clarify whether this denotes Hankel functions of order 0 and 1, and define the notation explicitly.
[Figures 6 and 7] The vertical axes of Figures 6 and 7 are not labeled; please specify that they plot the normalized radial functions f_A/f_0 and f_B/f_0 (Figure 6) and I_A/I_0 and I_B/I_0 (Figure 7), and define the normalization in the captions.
[Table III] The caption states that 'zero denoting coefficients that vanish as r → ∞,' but many entries in the table are blank rather than zero; please clarify the meaning of empty entries and of the entries marked '0'.
[Section IV.B after Eq. (42)] The statement that the delta potential 'can be replaced by a real and spatially localized function that vanishes for r > r0' is important for the regularization issue raised in the major comments; please indicate how the index computation and the zero-mode wavefunction would change under such a replacement.
[References] Reference [62] is an arXiv preprint; please provide its publication status or journal reference if available.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular step found: potentials, zero modes, and density are derived from tight-binding inputs without fitting; self-citations supply prior bulk-invariant context and a comparison benchmark, not the derivation itself.

full rationale

We traced the derivation chain from tight-binding inputs to the advertised predictions. The adatom and vacancy potentials (Eqs. 20-27 and 28-35) are obtained by Fourier/valley decomposition of the original tight-binding operators, not by fitting; the zero-mode equations (46)-(47) are solved in the angular-momentum basis, yielding the 1/r mode (48); and the density result (61) follows from an explicit first-order Green's function calculation (55)-(60), with higher orders treated in Appendix C. No fitted parameter is relabeled as a prediction, and no quantity is defined into existence by the index formula: the analytical index (38) is cited from external mathematical references and then evaluated on the kernels obtained from the zero-mode calculation. The main self-citations, [16] and [62], supply the bulk winding-number interpretation and a wavefront-dislocation comparison, respectively; they do not enter as fit parameters or as the origin of the calculated index or density. The potentially serious mathematical weakness, namely that the counted zero mode (48) is non-normalizable so that (38) is not a standard L2 Fredholm index on the non-compact plane without a regularization, is a well-posedness/rigor concern rather than a circularity: the paper does not define the index to count non-L2 modes by construction, and it does not simply rename the mode as the index. A small score of 2 reflects the presence of minor, non-load-bearing self-citations, not a derivation that reduces to its own inputs.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central claim rests on the continuum low-energy approximation, the point-like delta modeling of the vacancy, and the unsettled counting of non-normalizable zero modes in the analytical index. No parameters are fitted to data; vF, a, t are physical inputs.

assumptions (4)

domain assumption The continuum low-energy limit a δk << 1 replaces the tight-binding Hamiltonian by the Dirac operator (18).
Used systematically in Section III to convert tight-binding defect potentials into first-quantized continuum potentials (Eqs. 26, 34). Standard for graphene, but ignores lattice-scale effects at the defect core.
domain assumption The vacancy is represented by a delta-function potential δ(r - R0) whose exact shape is unimportant (Section III B, after Eq. 42).
The zero-mode counting and index result rely on this point-like limit; a finite-size vacancy could alter the kernel.
ad hoc to paper The analytical index on non-compact R^2 counts non-normalizable power-law zero modes as kernel elements (Section IV B, Eq. 48).
The zero mode ψ_zm ~ 1/r is not square-integrable; the paper does not provide a regularization. This is the load-bearing assumption behind Index D = -1.
standard math The relation |Index D| = |ν| = N_zm (bulk-edge correspondence) is taken from [16] and [70].
The paper does not derive the bulk winding number ν; it relies on the Atiyah-Singer index theorem and prior work.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Defects Potentials for Two-Dimensional Topological Materials." pith.science (2026). https://pith.science/paper/3NGINU57

@misc{pith2026250701530,
  author       = {Pith},
  title        = {Pith review of: Defects Potentials for Two-Dimensional Topological Materials},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3NGINU57}},
  note         = {Machine review of arXiv:2507.01530}
}

read the original abstract

For non-topological quantum materials, introducing defects can significantly alter their properties by modifying symmetry and generating a nonzero analytical index, thus transforming the material into a topological one. We present a method to construct the potential matrix configuration with the purpose of obtaining a non-zero analytical index, akin to a topological invariant like a winding or Chern number. We establish systematic connections between these potentials, expressed in the continuum limit, and their initial tight-binding model description. We apply our method to graphene with an adatom, a vacancy, and both as key examples illustrating our comprehensive description. This method enables analytical differentiation between topological and non-topological zero-energy modes and allows for the construction of defects that induce topology.

Figures

Figures reproduced from arXiv: 2507.01530 by the authors.

**Figure 3.** FIG. 3. Two triangular sublattices [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

**Figure 4.** FIG. 4. Representation of an adatom potential in honeycomb [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗

Figures from the paper (4 more)

**Figure 5.** Figure 5: FIG. 5. A single carbon atom removed from a graphene hon [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]

**Figure 6.** Figure 6: FIG. 6. Numerical evaluation of [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]

**Figure 7.** Figure 7: FIG. 7. The spatial decay of [PITH_FULL_IMAGE:figures/full_fig_p011_7.png]

**Figure 8.** Figure 8: FIG. 8. Numerical evaluation of [PITH_FULL_IMAGE:figures/full_fig_p013_8.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

Topological Winding Numbers from Wavefront Dislocations in Local Electronic Density
cond-mat.mes-hall 2025-08 conditional novelty 4.0 of 10

Wavefront dislocations in the local electron density around a vacancy in graphene carry the chiral winding number, measurable from STM images through a Fourier filtering protocol.
Engineering Topological Materials
cond-mat.mes-hall 2025-08 conditional novelty 4.0 of 10

Localized defects with a designed phase winding can shift a Dirac material into a topological class, e.g. from BDI to BDI or CII, generating zero modes.

Reference graph

Works this paper leans on

90 extracted references · 52 canonical work pages · cited by 2 Pith papers

[1]

The adatom potentialˆVad in (27) is given by (37) with just one non-zero diagonal element Haa

Valley coupling electronic density: adatom Here, we apply the previously described method to compute the first-order termδρ(1) ad (r) resulting from val- ley coupling due to an adatom [59], highlighting essential computational steps and the outcome (further details in Appendix A). The adatom potentialˆVad in (27) is given by (37) with just one non-zero di...
[2]

This is compared to the local density correc- tion (57) for an adatom potential to differentiate between topological vacancies and non-topological adatoms

Valley coupling electronic density: vacancy We use the previously mentioned method to calculate the first-order electronic density changeδρ(1) vac (r) due to a vacancy. This is compared to the local density correc- tion (57) for an adatom potential to differentiate between topological vacancies and non-topological adatoms. The vacancy potential ˆVvac in (...
[3]

B. J. Wieder, B. Bradlyn, J. Cano, Z. Wang, M. G. Vergniory, L. Elcoro, A. A. Soluyanov, C. Felser, T. Ne- upert, N. Regnault, and A. Bernevig B., Topological ma- terials discovery from crystal symmetry, Nature Reviews Materials 7, 196 (2024)

2024
[4]

M. Z. Hasan and C. L. Kane, Colloquium: Topological insulators, Rev. Mod. Phys.82, 3045 (2010)

2010
[5]

A.I.Khan, A.Ramdas, E.Lindgren, H.-M.Kim, B.Won, X. Wu, K. Saraswat, C.-T.Chen, Y. Suzuki, F.H. da Jor- nada, I.-K. Oh, and E. Pop, Surface conduction and re- duced electrical resistivity in ultrathin noncrystalline nbp semimetal, Science 387, 62 (2025)

2025
[6]

Kitaev, Fault-tolerant quantum computation by anyons, Annals of Physics303, 2 (2003)

A. Kitaev, Fault-tolerant quantum computation by anyons, Annals of Physics303, 2 (2003)

2003
[7]

Fu and C

L. Fu and C. L. Kane, Superconducting proximity effect and majorana fermions at the surface of a topological insulator, Phys. Rev. Lett.100, 096407 (2008)

2008
[8]

R.TakahashiandS.Murakami,Thermoelectrictransport in perfectly conducting channels in quantum spin hall systems, Phys. Rev. B81, 161302 (2010)

2010

Show all 90 references

[9]

C. L. Kane and E. J. Mele, Quantum spin hall effect in graphene, Phys. Rev. Lett.95, 226801 (2005)

2005
[10]

Xu and J

C. Xu and J. E. Moore, Stability of the quantum spin hall effect: Effects of interactions, disorder, and𭟋2 topology, Phys. Rev. B73, 045322 (2006)

2006
[11]

D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, Quantized hall conductance in a two- dimensional periodic potential, Phys. Rev. Lett.49, 405 (1982)

1982
[12]

Kohmoto, Topological invariant and the quantiza- tion of the hall conductance, Annals of Physics160, 343 (1985)

M. Kohmoto, Topological invariant and the quantiza- tion of the hall conductance, Annals of Physics160, 343 (1985)

1985
[13]

Altland and M

A. Altland and M. R. Zirnbauer, Nonstandard symme- try classes in mesoscopic normal-superconducting hybrid structures, Physical Review B55, 1142 (1997)

1997
[14]

A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Lud- wig, Classification of topological insulators and super- conductors in three spatial dimensions, Phys. Rev. B78, 195125 (2008)

2008
[15]

Kitaev, Periodic table for topological insulators and superconductors, AIP Conference Proceedings1134, 22 (2009)

A. Kitaev, Periodic table for topological insulators and superconductors, AIP Conference Proceedings1134, 22 (2009)

2009
[16]

J. C. Y. Teo and C. L. Kane, Topological defects and gapless modes in insulators and superconductors, Phys. Rev. B 82, 115120 (2010)

2010
[17]

C.-K. Chiu, J. C. Y. Teo, A. P. Schnyder, and S. Ryu, Classification of topological quantum matter with sym- metries, Rev. Mod. Phys.88, 035005 (2016)

2016
[18]

A. Goft, Y. Abulafia, N. Orion, C. L. Schochet, and E. Akkermans, Defects in graphene: A topological de- scription, Phys. Rev. B108, 054101 (2023)

2023
[19]

Kelly and N

K. Kelly and N. Halas, Determination of α and β site defects on graphite using c60-adsorbed stm tips, Surface Science 416, L1085 (1998)

1998
[20]

V. M. Pereira, F. Guinea, J. M. B. Lopes dos Santos, N. M. R. Peres, and A. H. Castro Neto, Disorder induced localized states in graphene, Phys. Rev. Lett.96, 036801 (2006)

2006
[21]

Amara, S

H. Amara, S. Latil, V. Meunier, P. Lambin, and J.-C. Charlier, Scanning tunneling microscopy fingerprints of point defects in graphene: A theoretical prediction, Phys. Rev. B 76, 115423 (2007)

2007
[22]

V. M. Pereira, J. M. B. Lopes dos Santos, and A. H. Castro Neto, Modeling disorder in graphene, Phys. Rev. B 77, 115109 (2008)

2008
[23]

J. J. Palacios, J. Fernández-Rossier, and L. Brey, Vacancy-induced magnetism in graphene and graphene ribbons, Phys. Rev. B77, 195428 (2008)

2008
[24]

N. M. R. Peres, S.-W. Tsai, J. E. Santos, and R. M. Ribeiro, Scanning tunneling microscopy currents on lo- cally disordered graphene, Phys. Rev. B 79, 155442 (2009)

2009
[25]

N. M. R. Peres, L. Yang, and S.-W. Tsai, Local density of states and scanning tunneling currents in graphene, New Journal of Physics11, 095007 (2009)

2009
[26]

Dutreix, L

C. Dutreix, L. Bilteanu, A. Jagannathan, and C. Bena, Friedel oscillations at the dirac cone merging point in anisotropic graphene and graphenelike materials, Phys. Rev. B 87, 245413 (2013)

2013
[27]

Dutreix and M

C. Dutreix and M. I. Katsnelson, Friedel oscillations at the surfaces of rhombohedral n-layer graphene, Phys. Rev. B 93, 035413 (2016)

2016
[28]

B. R. K. Nanda, M. Sherafati, Z. S. Popović, and S. Sat- pathy, Corrigendum: Electronic structure of the sub- stitutional vacancy in graphene: density-functional and green’s function studies, New Journal of Physics 15, 039501 (2013)

2013
[29]

Ovdat, Y

O. Ovdat, Y. Don, and E. Akkermans, Vacancies in graphene: Dirac physics and fractional vacuum charges, Phys. Rev. B102, 075109 (2020)

2020
[30]

M. M. Ugeda, I. Brihuega, F. Guinea, and J. M. Gómez- Rodríguez, Missing atom as a source of carbon mag- netism, Phys. Rev. Lett.104, 096804 (2010)

2010
[31]

Ovdat, J

O. Ovdat, J. Mao, Y. Jiang, E. Y. Andrei, and E. Akker- mans, Observing a scale anomaly and a universal quan- tum phase transition in graphene, Nature Communica- tions 8, 507 (2017)

2017
[32]

Mariani, L

E. Mariani, L. I. Glazman, A. Kamenev, and F. von Op- pen, Zero-bias anomaly in the tunneling density of states of graphene, Phys. Rev. B76, 165402 (2007)

2007
[33]

V. V. Cheianov and V. I. Fal’ko, Friedel oscillations, im- purity scattering, and temperature dependence of resis- tivity in graphene, Phys. Rev. Lett.97, 226801 (2006)

2006
[34]

Gade and F

R. Gade and F. Wegner, The n = 0 replica limit of u(n) and u(n)so(n) models, Nuc. Phys. B360, 213 (1991)

1991
[35]

Gade, Anderson localization for sublattice models, Nuclear Physics B398, 499 (1993)

R. Gade, Anderson localization for sublattice models, Nuclear Physics B398, 499 (1993)

1993
[36]

Fabrizio and C

M. Fabrizio and C. Castellani, Anderson localization in bipartite lattices, Nuc. Phys. B583, 542 (2000)

2000
[37]

X. Wu, X. Li, Z. Song, C. Berger, and W. A. de Heer, Weak antilocalization in epitaxial graphene: Evidence for chiral electrons, Phys. Rev. Lett.98, 136801 (2007)

2007
[38]

P.Ostrovsky, I.Gornyi,andA.Mirlin,Electrontransport in disordered graphene, Phys. Rev. B74, 235443 (2006)

2006
[39]

F. V. Tikhonenko, D. W. Horsell, R. V. Gorbachev, and A. K. Savchenko, Weak localization in graphene flakes, Phys. Rev. Lett.100, 056802 (2008)

2008
[40]

P. R. Wallace, The band theory of graphite, Phys. Rev. 16 71, 622 (1947)

1947
[41]

V. M. Pereira, A. H. Castro Neto, and N. M. R. Peres, Tight-binding approach to uniaxial strain in graphene, Phys. Rev. B80, 045401 (2009)

2009
[42]

Montambaux, F

G. Montambaux, F. Piéchon, J.-N. Fuchs, and M. O. Goerbig, Merging of dirac points in a two-dimensional crystal, Phys. Rev. B80, 153412 (2009)

2009
[43]

E. H. Lieb, Two theorems on the hubbard model, Phys. Rev. Lett. 62, 1201 (1989)

1989
[44]

C.A.CoulsonandG.S.Rushbrooke,Noteonthemethod of molecular orbitals, Mathematical Proceedings of the Cambridge Philosophical Society36, 193–200 (1940)

1940
[45]

Rigby and R

M. Rigby and R. Mallion, On the eigenvalues and eigenvectors of certain finite, vertex-weighted, bipartite graphs, Journal of Combinatorial Theory, Series B27, 122 (1979)

1979
[46]

Xiao and Y

R. Xiao and Y. X. Zhao, Revealing the spatial nature of sublattice symmetry, Nature Communications15, 3787 (2024)

2024
[47]

Nielsen and M

H. Nielsen and M. Ninomiya, The adler-bell-jackiw anomaly and weyl fermions in a crystal, Physics Letters B 130, 389 (1983)

1983
[48]

M. I. Katsnelson, Graphene carbon in two dimensions (Cambridge university press, 1979)

1979
[49]

Clifford, Preliminary sketch of biquaternions, Proceed- ings of the London Mathematical Society s1-4, 381 (1871)
[50]

P. A. M. Dirac and R. H. Fowler, The quantum theory of the electron, Proceedings of the Royal Society of Lon- don. Series A, Containing Papers of a Mathematical and Physical Character 117, 610 (1928)

1928
[51]

J. C. Slonczewski and P. R. Weiss, Band structure of graphite, Phys. Rev.109, 272 (1958)

1958
[52]

K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Electric field effect in atomically thin carbon films, Science 306, 666 (2004)

2004
[53]

A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, The electronic properties of graphene, Rev. Mod. Phys.81, 109 (2009)

2009
[54]

Cahangirov, M

S. Cahangirov, M. Topsakal, E. Aktürk, H. Şahin, and S. Ciraci, Two- and one-dimensional honeycomb struc- tures of silicon and germanium, Phys. Rev. Lett.102, 236804 (2009)

2009
[55]

J. Wang, S. Deng, Z. Liu, and Z. Liu, The rare two- dimensional materials with Dirac cones, National Science Review 2, 22 (2015)

2015
[56]

Malko, C

D. Malko, C. Neiss, F. Viñes, and A. Görling, Competi- tion for graphene: Graphynes with direction-dependent dirac cones, Phys. Rev. Lett.108, 086804 (2012)

2012
[57]

Huang, W

H. Huang, W. Duan, and Z. Liu, The existence/absence of dirac cones in graphynes, New Journal of Physics15, 023004 (2013)

2013
[58]

Dutreix, Impurity and boundary modes in the hon- eycomb lattice, Theses, Université Paris Sud - Paris XI (2014)

C. Dutreix, Impurity and boundary modes in the hon- eycomb lattice, Theses, Université Paris Sud - Paris XI (2014)

2014
[59]

Ando and T

T. Ando and T. Nakanishi, Impurity scattering in carbon nanotubes – absence of back scattering –, Journal of the Physical Society of Japan67, 1704 (1998)

1998
[60]

Yndurain, Effect of hole doping on the magnetism of point defects in graphene: A theoretical study, Phys

F. Yndurain, Effect of hole doping on the magnetism of point defects in graphene: A theoretical study, Phys. Rev. B 90, 245420 (2014)

2014
[61]

Dutreix, G.-H

C. Dutreix, G.-H. H., I. Brihuega, M. I. Katsnelson, C. Chapelier, and V. T. Renard, Measuring the berry phase of graphene from wavefront dislocations in friedel oscillations, Nature 574, 219 (2019)

2019
[62]

Bena, Effect of a single localized impurity on the lo- cal density of states in monolayer and bilayer graphene, Phys

C. Bena, Effect of a single localized impurity on the lo- cal density of states in monolayer and bilayer graphene, Phys. Rev. Lett.100, 076601 (2008)

2008
[63]

B. R. K. Nanda, M. Sherafati, Z. S. Popović, and S. Sat- pathy, Electronic structure of the substitutional vacancy in graphene: density-functional and Green’s function studies, New Journal of Physics14, 083004 (2012)

2012
[64]

Y.Abulafia, A.Goft, N.Orion,andE.Akkermans,Wave- fronts dislocations measure topology in graphene with defects (2023), arXiv:2307.05185

2023 arXiv
[65]

Dutreix, B

C. Dutreix, B. Matthieu, D. Pierre, and M. Fabrice, Wavefront dislocations reveal the topology of quasi-1d photonic insulators, Nature Communications 12, 3571 (2021)

2021
[66]

Nakahara,Geometry, topology and physics , Graduate student series in physics (Hilger, Bristol, 1990)

M. Nakahara,Geometry, topology and physics , Graduate student series in physics (Hilger, Bristol, 1990)

1990
[67]

Yankowsky, A

E. Yankowsky, A. Schwarz, and S. Levy,Quantum Field Theory and Topology, Grundlehren der mathematischen Wissenschaften (Springer Berlin Heidelberg, 2013)

2013
[68]

M. F. Atiyah and I. M. Singer, The index of elliptic op- erators on compact manifolds, Bulletin of the American Mathematical Society 69, 422 (1963)

1963
[69]

M. F. Atiyah and I. M. Singer, The index of elliptic op- erators: I, Annals of Mathematics87, 484 (1968)

1968
[70]

A. J. Niemi and G. W. Semenoff, Spectral asymmetry on an open space, Phys. Rev. D30, 809 (1984)

1984
[71]

Stone, Zero modes, boundary conditions and anoma- liesonthelatticeandinthecontinuum,AnnalsofPhysics 155, 56 (1984)

M. Stone, Zero modes, boundary conditions and anoma- liesonthelatticeandinthecontinuum,AnnalsofPhysics 155, 56 (1984)

1984
[72]

Hatsugai, Chern number and edge states in the integer quantum hall effect, Phys

Y. Hatsugai, Chern number and edge states in the integer quantum hall effect, Phys. Rev. Lett.71, 3697 (1993)

1993
[73]

N. Weik, J. Schindler, S. Bera, G. C. Solomon, and F. Evers, Graphene with vacancies: Supernumerary zero modes, Phys. Rev. B94, 064204 (2016)

2016
[74]

Akkermans and G

E. Akkermans and G. Montambaux,Mesoscopic Physics of Electrons and Photons (Cambridge University Press, Cambridge, 2007)

2007
[75]

C.-Y. Hou, C. Chamon, and C. Mudry, Electron frac- tionalization in two-dimensional graphenelike structures, Phys. Rev. Lett.98, 186809 (2007)

2007
[76]

Jackiw and S.-Y

R. Jackiw and S.-Y. Pi, Chiral gauge theory for graphene, Phys. Rev. Lett.98, 266402 (2007)

2007
[77]

O. V. Yazyev and L. Helm, Defect-induced magnetism in graphene, Phys. Rev. B75, 125408 (2007)

2007
[78]

P. O. Lehtinen, A. S. Foster, Y. Ma, A. V. Krashenin- nikov, and R. M. Nieminen, Irradiation-induced mag- netism in graphite: A density functional study, Phys. Rev. Lett. 93, 187202 (2004)

2004
[79]

Y. Ma, P. O. Lehtinen, A. S. Foster, and R. M. Nieminen, Magnetic properties of vacancies in graphene and single- walled carbon nanotubes, New Journal of Physics6, 68 (2004)

2004
[80]

M. A. H. Vozmediano, M. P. López-Sancho, T. Stauber, and F. Guinea, Local defects and ferromagnetism in graphene layers, Phys. Rev. B72, 155121 (2005)

2005
[81]

O. V. Yazyev, Magnetism in disordered graphene and irradiated graphite, Phys. Rev. Lett.101, 037203 (2008)

2008
[82]

O. V. Yazyev, Emergence of magnetism in graphene materials and nanostructures, Reports on Progress in Physics 73, 056501 (2010)

2010
[83]

Jiang, P

Y. Jiang, P. Lo, D. May, G. Li, G. Guo, F. Anders, T. Taniguchi, K. Watanabe, J. Mao, and E. Andrei, In- 17 ducing kondo screening of vacancy magnetic moments in graphene with gating and local curvature, Nature Com- munications 9, 2349 (2018)

2018
[84]

González-Herrero, J

H. González-Herrero, J. M. Gómez-Rodríguez, P. Mallet, M. Moaied, J. J. Palacios, C. Salgado, M. M. Ugeda, J.- Y. Veuillen, F. Yndurain, and I. Brihuega, Atomic-scale control of graphene magnetism by using hydrogen atoms, Science 352, 437 (2016)

2016
[85]

Nair, I.-L

R. Nair, I.-L. Tsai, M. Sepioni, O. Lehtinen, J. Keinonen, A. Krasheninnikov, A. C. Neto, M. Katsnelson, A. Geim, and I. Grigorieva, Dual origin of defect magnetism in graphene and its reversible switching by molecular dop- ing, Nature Communications4, 2010 (2013)

2013
[86]

Faccio, L

R. Faccio, L. Fernández-Werner, H. Pardo, C. Goyenola, O. N. Ventura, and A. W. Mombru, Electronic and struc- turaldistortionsingrapheneinducedbycarbonvacancies and boron doping, The Journal of Physical Chemistry C 114, 18961 (2010)

2010
[87]

M. D. Bhatt, H. Kim, and G. Kim, Various defects in graphene: a review, RSC Adv.12, 21520 (2022)

2022
[88]

A. M. Valencia and M. J. Caldas, Single vacancy defect in graphene: Insights into its magnetic properties from theoretical modeling, Phys. Rev. B96, 125431 (2017)

2017
[89]

R. R. Nair, M. Sepioni, I.-L. Tsai, O. Lehtinen, J. Keinonen, A. V. Krasheninnikov, T. Thomson, A. K. Geim, and I. V. Grigorieva, Spin-half paramagnetism in graphene induced by point defects, Nature Physics8, 199 (2012)

2012
[90]

Hubbard and B

J. Hubbard and B. H. Flowers, Electron correlations in narrow energy bands, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 276, 238 (1963)

1963

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.

[1] [1]

The adatom potentialˆVad in (27) is given by (37) with just one non-zero diagonal element Haa

Valley coupling electronic density: adatom Here, we apply the previously described method to compute the first-order termδρ(1) ad (r) resulting from val- ley coupling due to an adatom [59], highlighting essential computational steps and the outcome (further details in Appendix A). The adatom potentialˆVad in (27) is given by (37) with just one non-zero di...

[2] [2]

This is compared to the local density correc- tion (57) for an adatom potential to differentiate between topological vacancies and non-topological adatoms

Valley coupling electronic density: vacancy We use the previously mentioned method to calculate the first-order electronic density changeδρ(1) vac (r) due to a vacancy. This is compared to the local density correc- tion (57) for an adatom potential to differentiate between topological vacancies and non-topological adatoms. The vacancy potential ˆVvac in (...

[3] [3]

B. J. Wieder, B. Bradlyn, J. Cano, Z. Wang, M. G. Vergniory, L. Elcoro, A. A. Soluyanov, C. Felser, T. Ne- upert, N. Regnault, and A. Bernevig B., Topological ma- terials discovery from crystal symmetry, Nature Reviews Materials 7, 196 (2024)

2024

[4] [4]

M. Z. Hasan and C. L. Kane, Colloquium: Topological insulators, Rev. Mod. Phys.82, 3045 (2010)

2010

[5] [5]

A.I.Khan, A.Ramdas, E.Lindgren, H.-M.Kim, B.Won, X. Wu, K. Saraswat, C.-T.Chen, Y. Suzuki, F.H. da Jor- nada, I.-K. Oh, and E. Pop, Surface conduction and re- duced electrical resistivity in ultrathin noncrystalline nbp semimetal, Science 387, 62 (2025)

2025

[6] [6]

Kitaev, Fault-tolerant quantum computation by anyons, Annals of Physics303, 2 (2003)

A. Kitaev, Fault-tolerant quantum computation by anyons, Annals of Physics303, 2 (2003)

2003

[7] [7]

Fu and C

L. Fu and C. L. Kane, Superconducting proximity effect and majorana fermions at the surface of a topological insulator, Phys. Rev. Lett.100, 096407 (2008)

2008

[8] [8]

R.TakahashiandS.Murakami,Thermoelectrictransport in perfectly conducting channels in quantum spin hall systems, Phys. Rev. B81, 161302 (2010)

2010