{"id":"804a3c06-e42f-47c6-accd-c7aa3143fb39","arxiv_id":"2508.04927","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"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.","lead":"The paper proposes a strategy for turning normal, non-topological materials into topological ones by inserting localized defects such as vacancies or adatoms, which move the system to a different row of the tenfold classification. It shows lattice examples where a defect produces a nonzero winding number and even manufactures an effective spin-like quantum number from a spinless model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The BDI/CII guarantee rests on an unproved replacement of the vacancy's delta-derivative symbol by a smooth winding field; the bilayer's effective 2θ winding is asserted, not derived.","rationale":"The paper's central claim is a constructive recipe to turn a non-topological Dirac material into a prescribed topological class (BDI or CII) by engineering defects. For this to hold, the topological invariant computed from the continuum symbol (the winding of φ) must equal the index of the actual lattice Hamiltonian with the physical defect. The authors derive that the vacancy potential symbol is a derivative of a delta function, then replace it with a smooth φ(r) of the same angular winding, asserting that the radial profile is irrelevant. This is a homotopy assumption: the index is invariant under continuous deformations preserving ellipticity and symmetry, but no such deformation is exhibited, and the distributional nature of the delta-derivative is not addressed. The bilayer case is even less secure: the effective 2θ winding emerges from Löwdin partitioning of a matrix whose entries contain these distributions; the calculation is omitted, and the phase bookkeeping for vacancies on opposite sublattices is unclear. The numerical zero modes are supporting evidence but do not uniquely determine the winding number; a direct re-derivation or a real-space invariant computation would settle the issue. The AA-stacking inconsistency between main text and supplementary is a separate flaw affecting a secondary example, not the central BDI/CII construction. The reader's conditional verdict is appropriate: the concern is a missing proof, not a demonstrated counterexample, and a concrete re-derivation can resolve it. Therefore, I do not change the verdict.","tokens_in":21243,"tokens_out":9739,"duration_ms":111356,"concrete_test":"Independently re-derive the effective low-energy symbol for the AB-stacked bilayer with two aligned vacancies directly from the full tight-binding Hamiltonian (supplementary Eq. (46)) via the discrete Weyl transform, keeping the vacancy potential as -i e^{iθ}∂_r δ(r). Then apply Löwdin partitioning to the exact symbol and extract the phase winding of the off-diagonal block around the vacancy. If the winding is not exactly 2θ, or if the symbol is not elliptic after regularization, the continuum replacement and the predicted ν=±4 CII phase fail.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the construction 'guarantees the realization of a topological material in class BDI or CII' (main text after Eq. (8)) depends on replacing the physical vacancy potential, whose symbol is a derivative of a delta function (supplementary Eq. (35): V ∝ e^{iθ}(-i∂_r δ(r))), with an arbitrary smooth localized complex field φ(r)=|φ|e^{imθ}, declaring the radial profile irrelevant (supplementary after Eq. (37)). No proof is given that the full lattice symbol, including subleading terms, is homotopic to the Dirac form with the same winding for all models. For the monolayer, the computed winding ν=∓1 assumes the symbol is elliptic with an isolated zero set, but the delta-derivative is a distribution and its regularization is not justified. For the AB-stacked bilayer, the effective field with winding 2θ is obtained by Löwdin partitioning of the vacancy coupling matrix Φ (supplementary Eqs. (54)-(55)), which involves products of these distributions; the derivation is not shown, and the phase bookkeeping for vacancies on opposite sublattices is opaque. If the actual effective off-diagonal block has a different phase winding, the predicted CII phase with ν=∓4 and T²=−1 does not follow. The numerical zero modes are consistent but do not uniquely fix the invariant.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a constructive method for engineering topological phases by embedding point defects or textures into chiral Dirac materials. The method is framed through the authors' 'symbol of a Hamiltonian' approach, with the Atiyah–Singer index theorem used to relate the analytical index of an elliptic operator to a topological winding number. After laying out general rules for navigating the tenfold classification, the authors apply the construction to two-dimensional chiral materials: a monolayer honeycomb/brick-wall lattice with a vacancy is claimed to yield class BDI with winding number ν = ∓1; an AB-stacked bilayer with aligned vacancies is claimed to yield class CII, with an emergent valley-based T² = −1 symmetry and winding number ν = ∓4; AA stacking is claimed to preserve BDI but move the topological modes to finite energy. Exact-diagonalization results showing zero modes are presented in support.","tokens_in":21643,"tokens_out":11101,"duration_ms":138000,"significance":"If the construction is correct, it provides a concrete, design-oriented bridge between the tenfold classification and material engineering: starting from a non-topological Dirac material, a prescribed defect converts the system into a prescribed topological class. The paper's strengths are its explicit low-energy symbols, the concrete lattice realizations (monolayer vacancy, AB- and AA-stacked bilayers), and the numerical zero modes that are consistent with the claimed invariants. The monolayer vacancy result reproduces the authors' earlier PRB result, and the AB-stacked bilayer result is a new prediction. However, the central 'guarantee' rests on an unproved replacement of a distributional vacancy potential by a smooth complex scalar field, and the bilayer CII derivation is asserted rather than fully shown. The quantitative predictions ν = ∓1 and ν = ∓4 are therefore plausible but not established to the standard claimed by the paper.","major_comments":[{"comment":"The winding-number formula in Eq. (8) is written as an integral over dθ only, with the result ν = ∓1. The general formula in Supplementary Eq. (3) contains an integration over the full phase space, ∫ d^d k d^D r, normalized by S_{d+D}. As printed, Eq. (8) omits the momentum integration and the normalization factor. If the momentum integration is exactly trivial, that should be shown; if it contributes a constant factor, that factor must appear. The same issue affects the statement 'after a simple integration over the momentum variables' in Supplementary Eq. (40). This is load-bearing because the quantitative values ν = ∓1 and ν = ∓4 are the central predictions.","section":"Main text, Eq. (8); Supplementary Eq. (3)"},{"comment":"The vacancy potential is shown to have symbol V ∝ e^{iθ}(-i∂_r δ(r)). The authors regularize the delta function and then replace −i∂_rδ(r) by an arbitrary localized function φ(r), declaring 'only its angular dependence matters' and, in the bilayer, that the exact radial profile 'does not play a crucial role'. No proof or controlled asymptotic argument is given that this replacement preserves the topological index for the full lattice symbol, including subleading terms. For a distributional symbol, ellipticity and the index can depend on the regularization. This assumption is the mechanism behind the claimed 'guarantee' of class BDI or CII, so it needs either a proof of homotopy invariance under the replacement or a direct lattice-level check of the invariant for an actual vacancy.","section":"Supplementary Eqs. (35)-(37) and main text after Eq. (8)"},{"comment":"The derivation of the effective bilayer vacancy coupling Φ_eff is not shown. The statement 'Using Löwdin partitioning on Φ, which has the same structure as F, we find Φ_eff(r) = φ(r)(0 e^{2iθ}; e^{2iθ} 0)' is a single assertion. This step is responsible for the factor 2θ that converts the monolayer winding ∓1 into the bilayer winding ∓4, and it involves products of distributional vacancy fields on opposite sublattices. The phase bookkeeping is opaque. Please provide the explicit Löwdin computation for the 8×8 symbol, or otherwise demonstrate the 2θ phase from the lattice Hamiltonian.","section":"Supplementary Eqs. (54)-(55) and main text around Eq. (9)"},{"comment":"The claimed invariant ν = ∓4 for the quadratic-Dirac bilayer is not derived. The statement 'Similarly to (8), the associated winding number now becomes ν = ∓4' is insufficient: one factor of 2 is attributed to the quadratic dispersion and one to the combined phases, but the Jacobian integral for H_AB(k,r) is not displayed. An explicit calculation of the invariant from the symbol in Eq. (56) is required, especially to confirm that the two factors of 2 multiply rather than cancel.","section":"Main text after Eq. (9); Supplementary Eq. (56)"}],"minor_comments":[{"comment":"The Weyl transform convention differs between the main text and the supplement (R(j) vs. R(j)/2). The supplement notes the redefinition, but the two displayed formulas are not equivalent as written. Please use one convention consistently or state the relation explicitly.","section":"Main Eq. (1) vs. Supplementary Eq. (1)"},{"comment":"The symbol φ(r) is used both for the full complex field φ1 + iφ2 and for the radial amplitude, and Eq. (38) is not manifestly Hermitian as printed. Please clarify the notation and verify the conjugation of the off-diagonal blocks.","section":"Supplementary Eq. (38) and surrounding notation"},{"comment":"The relation NZM = |Index Q| = (1/2)|ν_Hf| is introduced without derivation and is not universal: if applied to the monolayer case |ν| = 1 it would give a half-integer. Please state the symmetry class and operator Q for which this relation holds.","section":"Main text after Eq. (10)"},{"comment":"The manuscript contains duplicated introductory blocks ('TOPAZ: the strategy!') and refers to figures that are not included in the text. Please ensure all figures and captions are present and remove extraneous material.","section":"General presentation"},{"comment":"The pairing of Dirac points is attributed to the Nielsen–Ninomiya theorem; this is a loose use of the theorem, which is traditionally about lattice fermion doubling. Please clarify the precise statement being used.","section":"Section on Nielsen–Ninomiya"}],"recommendation":"major_revision","confidential_remarks":"The paper's central claim is stated as a 'guarantee' but the proof level is below that standard: the profile-independence assumption and the bilayer 2θ phase are asserted rather than derived. The numerical zero modes support but do not uniquely fix the invariants. I would require the authors to either supply the explicit derivations and a lattice-level check of the index, or substantially soften the 'guarantee' language. The paper also leans heavily on the authors' own previous framework (refs. 9 and 53), which is not independently verified here; this is not a reason to reject, but it increases the importance of making the present derivation self-contained."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Chris, quick take on arXiv:2508.04927. The genuinely new piece is the AB-stacked bilayer with vacancies: from a spinless BDI material they get an effective CII phase with T^2=−1 and a predicted winding ±4, and their numerics show the two zero modes. That is a real constructive example, and if it holds up it is a nice route to spin-1/2 behavior from a spinless lattice. But the paper packages it inside a 'navigating the tenfold table' framework that is largely Teo–Kane plus their own symbol formalism; the monolayer vacancy result is already in their 2023 PRB. So the new result is narrower than the abstract suggests.\n\nThe strengths: the symbol-of-a-Hamiltonian method is used honestly on tight-binding models, the numerics are consistent with the index predictions, and the bilayer example is not in the cited literature. The idea that valley pseudospin can mimic T^2=−1 is worth taking seriously.\n\nThe soft spots are real. Eq. (8) and supplementary Eq. (3) write the winding as a single contour integral over θ; the momentum-space part of the phase-space integral is never shown. For the bilayer, the factor of 4 (two from the quadratic dispersion, two from the vacancy phases) is asserted, and the Löwdin partitioning step that gives the effective 2θ winding is a one-liner; the phase bookkeeping is opaque. There is also a direct contradiction between the main text and the supplementary on AA stacking: the main text says AA stacking breaks particle-hole symmetry and goes to class AI, but the supplementary derives a BDI symbol with complex field e^{iθ} and zero modes at ±t. One of those is wrong, or the two are talking about different configurations, and the paper does not say which.\n\nThe replacement of the vacancy's delta-derivative potential by a smooth localized φ(r)=|φ|e^{imθ} is plausible for the index, but it is not proved. Given that the guarantee claim rests on it, that needs a real argument.\n\nBottom line: the paper is not ready as is, but the AB bilayer example deserves referee time. A serious referee should ask for a complete phase-space winding calculation, a resolution of the AA-stacking contradiction, and a justification of the profile-independence of the index. With those, it could be a solid addition. I'd send it to review, but not accept in current form.","headline":"A genuinely new bilayer-vacancy example (BDI→CII, T^2=−1) is buried under a framework that mostly restates Teo–Kane and the authors' own previous work, and the key derivation is asserted rather than shown.","tokens_in":22083,"tokens_out":10128,"would_cite":false,"duration_ms":105237,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Adding engineered localized defects that connect the two valleys of a chiral Dirac material turns a non-topological phase into a topological one in a prescribed symmetry class, with the defect field's phase winding serving as the topologica","keywords":["topological phases","tenfold classification","defect engineering","winding number","chiral symmetry","graphene","Nielsen–Ninomiya theorem","Hamiltonian symbol"],"falsifier":"Exact diagonalization of a tight-binding honeycomb lattice with a vacancy whose potential is a finite-range function with the same $\\pm\\theta$ angular dependence but different radial decay: if the number of zero-energy modes (or the computed winding number) deviates from $\\mp1$ for monolayer, or from $\\mp4$ for aligned bilayer vacancies, while the phase winding is held fixed, the profile-independence assumption fails. An STM measurement of zero-mode counts on differently terminated vacancies would settle the same question experimentally.","tokens_in":21156,"feed_emoji":"🌀","tokens_out":15939,"duration_ms":144913,"temperature":0.7,"pith_summary":"A method to design topological phases on demand: embed a defect that couples the two valleys of a chiral Dirac material, generating a phase-winding complex scalar field that preserves chiral symmetry. The winding of $\\phi(r)=|\\phi|e^{im\\theta}$ around the defect becomes the topological invariant, moving the material from a nontopological BDI class into a topological class in BDI or CII, with zero modes localized at the defect. Two explicit constructions are given: a vacancy in monolayer graphene gives winding $\\mp1$ in BDI; aligned vacancies in AB-stacked bilayer graphene give winding $\\mp4$ in CII and generate an effective spin-1/2 degree of freedom. If correct, this turns the tenfold classification from a diagnostic table into an engineering recipe, applicable to any lattice whose low-energy spectrum obeys the Nielsen–Ninomiya theorem.","feed_headline":"One defect can make graphene topological","feed_subtitle":"A single vacancy gives graphene a protected winding; two aligned vacancies create an artificial spin.","key_machinery":"The Hamiltonian symbol $H(k,r)$ obtained from a discrete Weyl transform, together with the defect-induced complex scalar field $\\phi(r)=|\\phi|e^{im\\theta}$. The field's phase winding around the defect is the topological charge: a vacancy adds $\\phi_1\\sigma_x\\otimes\\tau_x+\\phi_2\\sigma_x\\otimes\\tau_y$ to the Dirac-point symbol, and the winding number $\\nu=\\frac{1}{2\\pi}\\oint d\\theta\\,\\partial_\\theta\\arg\\phi$ follows from the Atiyah–Singer index theorem as the count of localized zero modes. The accessible symmetry classes are fixed by $s=p-q\\bmod 8$, the parity condition $\\delta-s\\equiv 0 \\pmod 4$ selecting integer ($\\mathbb{Z}$) topology.","core_discovery":"The central claim: a topological phase with prescribed symmetry class can be engineered by localized defects that couple low-energy modes around Dirac points. The defect adds a complex scalar field $\\phi(r)=|\\phi|e^{\\mp i\\theta}$ (or $e^{\\mp 2i\\theta}$ for quadratic points), preserving chiral symmetry and connecting valleys, so that the symbol $\\mathbf{h}(k)\\cdot\\boldsymbol{\\gamma}$ becomes elliptic with winding number $\\nu_{H_f}=\\mp1$ (monolayer) or $\\mp4$ (bilayer). By the Atiyah–Singer index theorem these windings are tied to zero modes localized at the vacancies—two for the bilayer's winding 4—realizing BDI→BDI and BDI→CII transitions. The latter creates an effective $T^2=-1$ spin degree","pith_inferences":["A direct test of the profile-independence assumption is available: compute the winding number in a tight-binding honeycomb lattice for several vacancy models (removed site, renormalized hoppings, hydrogen-passivated) with identical $\\pm\\theta$ phase structure; the claim predicts the index stays $\\pm1$, and a deviating index would refute the continuum replacement.","The same mechanism should transfer to artificial Dirac platforms—photonic lattices, acoustic metamaterials, and cold-atom optical lattices—where 'defects' can be carved with controlled phase windings, offering a tabletop implementation of tenfold-navigation.","The two zero modes in the CII bilayer example, protected by chiral symmetry, suggest a natural topologically protected pseudospin qubit; whether the two localized modes can be addressed and braided without breaking the symmetry is an open extension.","Because the paper links the tenfold classification to two-qubit entanglement, the engineered effective spin might also serve as a physical resource for entanglement-based quantum gates, although the specific circuit-level construction is not given."],"forward_implications":["A non-topological chiral Dirac material with paired valleys (graphene, brickwall, Kagome, Mielke lattices) becomes topological under a single vacancy: monolayer lands in class BDI with winding $\\mp1$, bilayer in class CII with winding $\\mp4$.","Effective spin-1/2 degrees of freedom are generated on demand: the AB-stacked bilayer construction yields time-reversal symmetry with $T^2=-1$ starting from a system with $T^2=+1$, without any intrinsic spin-orbit coupling.","Topological zero modes localize at the vacancy sites themselves, so the bulk–edge correspondence holds with defects playing the role of quasi-boundaries inside the material.","The method's scope extends to non-lattice systems: the same defect-field construction works for quantum graphs, where the relevant operators are elliptic and the index theorem still applies.","Chiral materials satisfying the Nielsen–Ninomiya theorem form a general family for this 'tenfold navigation'; the table of possible $(d,s,D)$ combinations up to $d=3$ specifies which defect type (point, line, surface) realizes a $\\mathbb{Z}$ phase."],"supporting_citations":[{"why":"Defines the ten symmetry classes (T, P, C) that the method navigates; the classification framework's starting point.","marker":"[5]"},{"why":"Provides the periodic table and the Z/Z2 entries that determine which phases are reachable in each (δ,s).","marker":"[6]"},{"why":"Independent classification of 3D topological insulators and superconductors that completes the periodic table.","marker":"[7]"},{"why":"Introduces codimension δ = d−D for defects, the dimensionality reduction that lets a material change class.","marker":"[8]"},{"why":"Establishes the vacancy-in-graphene complex scalar field and its winding number, the template this paper generalizes to a systematic method.","marker":"[9]"},{"why":"The Atiyah–Singer index theorem reference used to link the topological index to the number of localized zero modes.","marker":"[14]"},{"why":"Nielsen–Ninomiya theorem guarantees Dirac points come in pairs of opposite chirality, the structural premise of the defect-coupling construction.","marker":"[32]"},{"why":"Löwdin partitioning used to derive the effective bilayer Hamiltonian whose quadratic dispersion yields the CII winding ±4.","marker":"[54]"}],"fun_headline_variants":["A single vacancy flips graphene's topology","Defects steer graphene into topological phases","One hole in graphene creates a topological band","Graphene becomes topological with a vacancy","Local defect toggles graphene's topological class"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"Replacing the physical vacancy potential (whose symbol behaves like a derivative of a delta function) with an arbitrary localized complex scalar field $\\phi(r)=|\\phi|e^{im\\theta}$, assuming only the phase winding matters while the exact radial profile does not alter the topological index; if this substitution fails for realistic vacancy shapes, the computed winding numbers do not describe the actual lattice.","fun_headline_variants_meta":{"raw":{"variants":["A single vacancy flips graphene's topology","Defects steer graphene into topological phases","One hole in graphene creates a topological band","Graphene becomes topological with a vacancy","Local defect toggles graphene's topological class"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000186,"raw_usage":{"total_tokens":1101,"prompt_tokens":625,"completion_tokens":476,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":369,"completion_tokens_details":{"reasoning_tokens":420}},"tokens_in":369,"tokens_out":476,"duration_ms":5733,"temperature":1.0,"reasoning_tokens":420,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T23:42:36.939683+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Exact diagonalization of a tight-binding honeycomb lattice with a vacancy whose potential is a finite-range function with the same $\\pm\\theta$ angular dependence but different radial decay: if the number of zero-energy modes (or the computed winding number) deviates from $\\mp1$ for monolayer, or from $\\mp4$ for aligned bilayer vacancies, while the phase winding is held fixed, the profile-independence assumption fails. An STM measurement of zero-mode counts on differently terminated vacancies would settle the same question experimentally.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the vacancy-in-graphene complex scalar field and its winding number, the template this paper generalizes to a systematic method."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Löwdin partitioning used to derive the effective bilayer Hamiltonian whose quadratic dispersion yields the CII winding ±4."}],"review_version":1}