{"id":"85f9492d-3f9c-4bf3-84f0-5413d90c93fd","arxiv_id":"2508.19128","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"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.","lead":"The paper proposes reading off a material's topological winding number from the twists, called dislocations, in the electron density pattern around a defect, using scanning tunneling microscope images. The authors demonstrate the recipe on a vacancy and an adatom in graphene and argue it offers a standardized way to certify topological states.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Filtering protocol's hand-set choices are not shown to be robust, so the extracted winding may not be the topological index.","rationale":"The paper's central claim is that a standardized Fourier-filtering protocol extracts the topological winding number from δρ. The reader identified the weakest assumption as the filtering protocol's hand-set choices; I agree. The proposed concrete check directly varies those choices on the same data and on synthetic data with known winding. If the results are invariant, the concern is resolved. If not, the 'definitive' framing is unsupported. I do not see an internal inconsistency: the index-theory chain is inherited from prior work, and the vacancy example is plausible. The missing support is empirical robustness, not mathematical contradiction. The paper does provide tight-binding numerics and a demonstration on published STM data, which are real evidence; the concern is about generality, not fabrication. Therefore the reader's CONDITIONAL verdict should stand; no change.","tokens_in":8626,"tokens_out":8669,"duration_ms":98578,"concrete_test":"Recompute the winding (4) from the filtered inverse FT of the tight-binding vacancy data (Fig. 1) with: (i) filter windows 4×4, 8×8, and 16×16 pixels; (ii) each of the three inequivalent satellite-pair choices; (iii) r-prefactor varied (none, r, r^2). Also run the identical protocol on synthetic data generated from Eq. (7) with added noise. If the extracted winding changes or the dislocation disappears for any choice, the readout is not protocol-independent. As a second check, apply the same protocol to a chiral lattice model with a defect of known index ±2; the extracted winding should be ±2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Paper's central chain (2) is operationalized by reading ν in (4) from a dislocation in δρ after a Fourier filter. The load-bearing step is the hand-set protocol in Supplementary 'Filtering intervalley scattering' and Fig. 1b: six satellite peaks appear, 'to avoid overlapping signals of the three duplicated pairs, we filter around two out of the six'; an '8-pixel window' keeps only selected peaks; the inverse FT is taken after 'phase information being reinstated by incorporating the matching phase at each pixel.' For the extracted winding to equal Index D, the intervalley phase χ(r,n) of Eq. (7) must survive filtering without contamination. No test is shown varying the peak pair, window size, or r-prefactor; no leakage estimate from the excluded four peaks, dc component, or Friedel rings; no uncertainty on the experimental readout (Fig. 4b). Since the paper claims a 'definitive' and 'general' method, and concedes only simple examples are shown with broader cases 'planned,' this missing robustness leaves the central claim conditional.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a general method to extract the topological winding number ν of chiral-symmetric Hamiltonians from wavefront dislocations in the local electronic density δρ(r). The central chain is Nzm = |Index D| = |ν| (Eq. 2), with ν read from a contour integral of the phase gradient of the intervalley LDOS (Eq. 4). The method is demonstrated on a graphene vacancy (topological, ν = ±1) versus an adatom (non-topological), using tight-binding simulations and a published STM image (Fig. 4). The authors argue that the dislocation pattern in δρ after a standardized Fourier filtering step directly measures the topological invariant.","tokens_in":8913,"tokens_out":2823,"duration_ms":32375,"significance":"If the proposed protocol is robust, it would provide a rare observable—real-space imaging of the winding number—for chiral topological phases, and would potentially extend to other platforms with intervalley/defect interference. The manuscript is commendably explicit about numerical parameters (dE = 10^-4 t, i0+ = 10^-5 t, PBC, V0 = -3t, ε_F = 0.6t) and the vacancy/adatom contrast is physically plausible. However, the two load-bearing steps—the adopted form of the intervalley phase and the hand-tuned Fourier filter—are not independently validated, so the central claim of a general and definitive method remains conditional.","major_comments":[{"comment":"The extracted winding number is not shown to be robust to the filter's hand-set choices. The protocol retains two of the six satellite peaks, uses an 8-pixel window, and reinstates phases after inverse filtering; no test varies the peak pair, window size, or the r (or r²) prefactor, and no leakage estimate is given for the excluded peaks, dc component, or Friedel rings. Since Eq. (4) requires contour-independent phase data, a demonstration that the readout ν is invariant under these choices is essential to support the method's generality.","section":"Supplementary 'Filtering intervalley scattering' and Fig. 1b"},{"comment":"The phase χ(r,n) = ΔK·r + nθ with n = 2 is imported from the same group's ref. 14. This already fixes an angular winding of 2, so the subsequent identification ν = ±1 in Fig. 3b and Fig. 4b is substantially predetermined by the assumed analytic form rather than independently measured. To support the claim that δρ(r) 'allows computing' ν, the manuscript should either derive Eq. (7) within this paper, or test the protocol against a system whose winding is known from an independent method and verify that the extracted ν does not depend on the assumed n.","section":"Eqs. (7)–(8)"},{"comment":"The equality Nzm = |Index D| = |ν| is asserted, but the paper does not independently compute Index D for the vacancy potential; it cites ref. 14 for Index D = ν = ±1. For a genuine measurement, the numerical tight-binding model should allow a direct count of zero modes (or a direct computation of Index D) to be compared with the wavefront-dislocation readout. Without such a check, the measured ν and the zero-mode count are connected only by an external assumption, leaving the central chain untested.","section":"Section 'Interferences Measure Topology', text near Eq. (2)"}],"minor_comments":[{"comment":"\"T echnion\" contains a spurious space; should read \"Technion.\"","section":"Author affiliation"},{"comment":"The sentence \"variations in local electronic density δρ(r) in (1) as a result of external potentials\" should refer to Eq. (3), not Eq. (1).","section":"Main text, paragraph after Eq. (5)"},{"comment":"The statement that the FT filter \"was done on the original 36 × 57 numerical results\" is unclear: the figure shows a cropped region, so the description should specify whether filtering and inverse FT are performed on the full lattice before cropping.","section":"Fig. 1 caption"},{"comment":"The adatom case is described as having \"no dislocations near the adatom\" after filtering, but Fig. 6 shows alternating ±1 dislocations at large distances with zero net winding. This distinction should be explicitly reconciled in the caption or text to avoid confusion.","section":"Figs. 5 and 6"}],"recommendation":"major_revision","confidential_remarks":"The manuscript builds heavily on the same group's earlier work (refs. 12–14), and the 'measurement' is in effect a filter-based visualization of the phase expression derived there. The missing robustness analysis is the main blocker; it can be fixed in revision but requires substantial additional numerical experiments rather than a simple textual change."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real new thing here is the concrete protocol: multiply δρ by r, take an FFT, keep one diametrically opposed pair of the six intervalley satellite peaks, filter with an 8-pixel window, inverse-transform with the phase reinstated, and read the dislocation winding. Applied to a graphene vacancy it gives ±1, and to an adatom it gives canceling ±1 pairs for net zero. That contrast is physically sensible and the tight-binding numerics are honestly documented (dE, i0+, PBC, V0, εF all stated). Re-analyzing the published STM image is a nice demonstration, even if the winding there is read by eye.\n\nWhere I part company with the paper is the frame. Equations (7)-(9) and the index chain come from the group's own earlier work (refs. 12-14), so the paper's genuine contribution is the filtering recipe, not the underlying identification. The recipe itself has hand-set choices—which two of the six peaks, the 8-pixel window, the r-prefactor—and no sensitivity analysis is shown. Nothing rules out leakage from the excluded peaks or the dc/Friedel components contributing comparable phase gradients. The experimental readout has no uncertainty, and the orbital-smoothing decay length a from the supplementary is not given a value. The claims in the abstract about a 'general', 'definitive' method and quantum-technology building blocks are not supported by two simple examples; the authors themselves concede broader cases are planned.\n\nI don't think this is a fatal problem. The core demonstration is coherent and the protocol could be genuinely useful if it survives robustness checks. The circularity worry is real but not disqualifying: yes, the n=2 angular index is built into the imported phase χ, so the dislocation winding is predetermined once you assume (7). But the paper's value is in showing that this phase appears in realistic local-density maps after filtering, which is an empirical claim worth checking. Give the authors a serious referee, not a desk reject. The referee should demand: filter-window and peak-pair variation, a leakage estimate, an error bar on the experimental dislocation, the value of a, and a tone-down of the general claims. If those are fixed, this becomes a solid methods paper for the AIII/defect community.","headline":"A clearly specified STM-image protocol for reading a chiral winding number from local density dislocations, demonstrated on a graphene vacancy and an adatom; the specific filter choices are new but the underlying formulas are inherited, and the robustness and scope claims outrun the evidence.","tokens_in":9459,"tokens_out":2610,"would_cite":false,"duration_ms":27488,"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":"A wavefront dislocation in the local electronic density carries the winding number of a chiral topological state, readable directly from STM images.","keywords":["topological winding number","wavefront dislocation","local electronic density","scanning tunneling microscopy","chiral symmetry","index theorem","graphene vacancy","intervalley interference"],"falsifier":"Use the same STM image of a graphene vacancy and filter around each of the other two pairs of satellite peaks, with different filter-window sizes; if the winding number from the contour integral changes, the read-out is not a protected invariant. Alternatively, simulate the vacancy density with two different radial prefactors (r vs r²) and check that the winding stays ±1.","tokens_in":8458,"feed_emoji":"🌀","tokens_out":12079,"duration_ms":102002,"temperature":0.7,"pith_summary":"The paper argues that for chiral-symmetric Hamiltonians, the integer topological invariant can be extracted from an STM-like image of the local electronic density rather than from edge-state transport or band-structure calculations. The load-bearing identity is Nzm=|Index D|=|ν|, which identifies the number of topological zero modes with both the analytic index of the off-diagonal block D and the bulk winding number. The paper shows that, for a graphene vacancy, the intervalley interference term makes the phase of the density wind by ±1 around the defect, producing a visible wavefront dislocation; integrating the phase gradient along any contour around the defect returns the winding number. The same filtering protocol applied to a chiral-symmetry-breaking adatom yields Friedel oscillations and no stable dislocation, showing the method separates topological from non-topological defects. This gives a direct, imaging-based route to identify and characterize topological quantum states.","feed_headline":"Wavefront dislocations in STM maps reveal a material's winding number","feed_subtitle":"A contour-independent integral over the density phase counts zero modes and separates a vacancy from an adatom.","key_machinery":"The central object is the phase χ(r,n)=ΔK·r+nθ living in the intervalley interference term of the density. The argument turns on the identity Nzm=|Index D|=|ν|, where D is the off-diagonal block of a chiral Hamiltonian, and on the contour integral (1/2π)∮∇χ·dr=ν. The wavefront dislocation is the real-space manifestation of this identity: a single extra wavefront terminating at the defect, whose winding number counts the topological zero modes. The Fourier-filtering protocol extracts exactly this intervalley phase from an STM image, making the integer observable.","core_discovery":"The central claim is that the local electronic density contains a direct signature of the topological winding number. For a chiral Hamiltonian, the paper uses the index-theorem chain Nzm=|Index D|=|ν| to link the number of zero modes to the topological invariant ν. For a vacancy in graphene, the intervalley contribution to the density variation takes the form δρΔK(r)=F(r)[cos(ΔK·r)−cos(ΔK·r+2θ)], whose phase χ(r,n)=ΔK·r+nθ has a winding (1/2π)∮∇χ·dr=±1 around the defect. Solving Dψ=0 gives Index D=ν=±1 and Nzm=1, so the observed dislocation pattern is interpreted as a topological zero mode. The paper demonstrates that this pattern appears both in tight-binding simulations and in existing STM","pith_inferences":["A built-in consistency test follows from the six satellite peaks: filtering around each of the three diametrically opposed pairs should return the same ν; the paper does not report this cross-check, so it is a natural next validation.","The method should extend to other local chiral defects and to classical wave systems (photonic, acoustic, cold-atom) where a local field is imaged, since the dislocation mechanism relies only on interference and chirality.","The protocol's practical limit is likely phase reconstruction near the dislocation core, where the intervalley amplitude is small; injecting controlled noise into synthetic density maps would show how much filtering robustness the method has."],"forward_implications":["STM images of a single defect become a direct probe of the winding number, eliminating the need for edge-state transport or Hall measurements.","A standardized Fourier-filtering step separates topological defects (vacancy, ν=±1) from non-topological ones (adatom, no stable winding).","The protocol works on already published STM data, so it can be applied to existing images without new instrumentation.","The same index-theorem chain, applied to other chiral-symmetric Hamiltonians, would give a local real-space route to their winding numbers."],"supporting_citations":[{"why":"Supplies the index-theorem statement Index D = ν for chiral Hamiltonians and the topological description of graphene defects that the density readout is built on.","marker":"12"},{"why":"Gives the intervalley density formula δρΔK(r)=F(r)[cos χ(r,0)-cos χ(r,2)] and the calculation Dψ=0 that fixes Index D=±1 for the vacancy and describes the adatom case.","marker":"14"},{"why":"Provides the published STM data of a graphene vacancy that the method is demonstrated on.","marker":"42"},{"why":"Establishes the wavefront-dislocation technique for reading Berry phase from Friedel oscillations, which the present method adapts to intervalley interference and uses as the non-topological adatom comparison.","marker":"43"},{"why":"Foundational index theorem connecting analytic and topological indices, the mathematical basis of the chain Nzm=|Index D|=|ν|.","marker":"5"},{"why":"Nielsen-Ninomiya no-go theorem, cited to guarantee the fermion-doubling two-valley structure that the intervalley interference pattern requires.","marker":"40"}],"fun_headline_variants":["Density dislocations expose a material's winding number","STM density phase knots reveal topological topology","Wavefront defects in local density count winding","A material's topology read from density wavefronts","Contour-free dislocations in STM yield winding number"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The method assumes that the filtering step pulls out one clean interference pattern and that nothing else (background oscillations, the other satellite signals, or numerical noise near low-amplitude points) adds comparable phase gradients; if that fails, the winding number read from the image is not protected.","fun_headline_variants_meta":{"raw":{"variants":["Density dislocations expose a material's winding number","STM density phase knots reveal topological topology","Wavefront defects in local density count winding","A material's topology read from density wavefronts","Contour-free dislocations in STM yield winding number"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000274,"raw_usage":{"total_tokens":1474,"prompt_tokens":738,"completion_tokens":736,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":482,"completion_tokens_details":{"reasoning_tokens":674}},"tokens_in":482,"tokens_out":736,"duration_ms":8306,"temperature":1.0,"reasoning_tokens":674,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T15:57:23.201379+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Use the same STM image of a graphene vacancy and filter around each of the other two pairs of satellite peaks, with different filter-window sizes; if the winding number from the contour integral changes, the read-out is not a protected invariant. Alternatively, simulate the vacancy density with two different radial prefactors (r vs r²) and check that the winding stays ±1.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the index-theorem statement Index D = ν for chiral Hamiltonians and the topological description of graphene defects that the density readout is built on."},{"cited_title":"Defects Potentials for Two-Dimensional Topological Materials","cited_arxiv_id":"2507.01530","evidence_quote":"Gives the intervalley density formula δρΔK(r)=F(r)[cos χ(r,0)-cos χ(r,2)] and the calculation Dψ=0 that fixes Index D=±1 for the vacancy and describes the adatom case."},{"cited_title":"M., Brihuega, I., Guinea, F","cited_arxiv_id":null,"evidence_quote":"Provides the published STM data of a graphene vacancy that the method is demonstrated on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the wavefront-dislocation technique for reading Berry phase from Friedel oscillations, which the present method adapts to intervalley interference and uses as the non-topological adatom comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Foundational index theorem connecting analytic and topological indices, the mathematical basis of the chain Nzm=|Index D|=|ν|."},{"cited_title":"& Ninomiya, M","cited_arxiv_id":null,"evidence_quote":"Nielsen-Ninomiya no-go theorem, cited to guarantee the fermion-doubling two-valley structure that the intervalley interference pattern requires."}],"review_version":1}