{"id":"1dc2aa39-0aa5-49a3-be7c-cb730b8e332e","arxiv_id":"2412.13042","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Impurity-induced density ripples in a nodal superconductor carry a vortex whose integer vorticity equals the topological winding difference of the wavefunctions at the two connected nodes, under explicit conditions on the impurity and the STM tip.","lead":"This paper shows that the winding of electron wavefunctions around the zero-energy points of a superconductor can be read off from the ripples an impurity creates, as seen by a scanning tunneling microscope. For the material NbSe2 in a magnetic field, it predicts that a magnetic impurity and a spin-polarized tip reveal the winding of all twelve nodes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The vortex-winding correspondence hinges on projecting the Green's function onto the cone band after moving all winding into n=1; the full unprojected check that would validate this is deferred, so the measured vortex is not yet tied to the physical node winding.","rationale":"The reader's weakest_assumption is exactly the same load-bearing point I find: the projection of the Green's function onto the cone band after a gauge choice that places all winding in n=1. I agree. The full Green's function in Eq. (8) includes all bands and is gauge-invariant; the projected form (13)-(16) is not, and the paper explicitly defers the numerical check in Appendix C1. Nothing in the analytic derivation rules out the gapped band contributing to the vortex charge in the exact LDOS, nor does it show that the measured charge equals the total winding W=Σ_n W_n rather than the projected W_{n=1}. This is a correctness risk, not a stylistic issue. Because the reader already set CONDITIONAL with this same concern, no verdict movement is needed; the condition should be a full unprojected calculation, ideally with a gauge-rephasing check. The other points raised by the reader (main-text vs appendix chirality condition, missing code/data) are real but secondary; the projection/gauge issue is the one that, if wrong, invalidates the central claim.","tokens_in":25621,"tokens_out":14984,"duration_ms":157389,"concrete_test":"Perform a full, unprojected T-matrix calculation of the spin-LDOS for the NbSe2 model of Eq. (39) with impurity Vx and tip Mx,My (parameters of Fig. 2), using the exact 4-band Green's function in the original spin/Nambu basis (no linearization, no projection). Fourier transform and extract the phase winding around each ΔKij; compare with the projected n=1 prediction ΔW_n=1 and with the full-model total winding ΔW=Σ_n W_n. Repeat the exact calculation after a smooth rephasing of the gapped-band eigenvectors that shifts W_2 by ±1: the physical LDOS should be invariant, while the projected prediction changes. Matching ΔW with invariance would validate the projection; any deviation or gauge sensitivity would falsify the central vortex-winding claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III.B/III.E claims that a non-zero winding difference ΔW between two nodes forces a vortex/dislocation of charge ±ΔW in the LDOS, with ΔW defined through the phases θ_n of the Q-matrix eigenvalues (Eqs. (7), (12), (15), (29)). The load-bearing premise is that the phase of the single isolated cone band n=1, after fixing the chiral-basis gauge, carries the full topological information. But Eq. (7) defines the node winding only as a sum over all bands, W=Σ_n W_n. As Appendix C1 states, the Q-matrix basis is not unique: 'we may reshuffle the non-trivial winding between the different bands n. The projection to n = 1 may therefore remove some winding contributions from the Green's function.' The authors fix this by choosing a basis in which all winding sits in n=1 and then project the Green's function to that band (Eqs. (13)-(16)). Nothing in the paper proves that the exact, gauge-invariant spin-LDOS produced by the full Green's function (8) with physical impurity/tip operators has a vortex at ΔKij whose charge equals W_{n=1} in that chosen gauge, rather than some other combination of band windings and matrix-element weights. Since the phases of BdG eigenstates are arbitrary, the 'winding of a node' is gauge-dependent; the paper's canonical choice (gapped bands have W_{n>1}=0) is plausible but not demonstrated to be the quantity an STM measures. The authors themselves defer the decisive check: 'a full (numerical) calculation of the Green's function including the gapped bands, without projection, may be used to confirm the validity of the projection procedure' (Appendix C1). Until that check is performed, the central claim that the observed vortex equals the topological winding difference is conditional.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a general theory for extracting the difference of chiral winding numbers of two Dirac nodes in a 2D nodal chiral system from impurity-induced LDOS modulations. Using a chiral-basis Green's function and the T-matrix/Born approximation, it derives conditions on the chirality preferences of an STM tip and an impurity under which the LDOS at a scattering vector ΔKij displays a wavefront dislocation in real space and a phase vortex in momentum space with charge equal to the winding difference ΔW = WKi - WKj. The method is applied to the proposed nodal Ising superconducting phase of monolayer NbSe2 under an in-plane magnetic field, where the authors predict that a magnetic impurity combined with a spin-polarized STM tip can reveal the winding of all 12 nodes. The paper also discusses anisotropy corrections, the relative robustness of momentum-space vortices versus real-space dislocations, and the experimental feasibility including a quantitative Born-approximation estimate.","tokens_in":25868,"tokens_out":10366,"duration_ms":93250,"significance":"If the central correspondence is fully substantiated, the paper provides a concrete experimental route to measure a lower-dimensional topological invariant—the winding of superconducting nodes—via STM, which is currently lacking. The work extends the graphene wavefront-dislocation technique to superconductors and other chiral nodal systems, and it gives explicit, falsifiable predictions for a specific material (NbSe2). Strengths include the detailed analytical derivation of the conditions, the generalization to anisotropic cones and complex chirality parameters, the careful treatment of experimental resolution, and the honest enumeration of approximations. The authors also provide a quantitative estimate for the Fe impurity in NbSe2, showing that the Born approximation is marginally reasonable (G0V ≈ 0.3). These features make the paper a valuable contribution if the identified gaps are addressed.","major_comments":[{"comment":"The general conditions for a non-zero dislocation charge are misstated for complex chirality preferences. With IΔ = h1 χM and I−Δ = h1 χV (Eq. (30)), the requirement |IΔ| ≠ |I−Δ| in Eq. (28) is equivalent to |χM| ≠ |χV|, not to χM ≠ ±χV as claimed in Eq. (36). Equation (34b) as written only excludes χM = χV. For instance, χM = 1, χV = i satisfies both Eq. (34b) and Eq. (36) yet has equal magnitudes, so the dislocation is absent. The correct conditions for complex χ are Eqs. (A8a)-(A8b) of Appendix A 2; the main text should either adopt them or explicitly restrict the derivation to real-valued χ.","section":"III.C, Eqs. (27)-(36)"},{"comment":"The central identification of the measured LDOS vortex charge with the physical node winding difference ΔW is not yet established, because the derivation projects the Green's function onto the single cone band n = 1 after fixing a gauge in which all winding sits in that band. As the authors state in Appendix C1, the winding number is defined only as a total over bands (Eq. (7)) and \"we may reshuffle the non-trivial winding between the different bands n\"; the projection to n = 1 can therefore remove winding contributions. The exact Green's function (8) includes the gapped bands, and the impurity/tip operators act in the full Hilbert space, so the weight of each band in the LDOS need not match the n = 1 projection. The paper explicitly defers the decisive check to \"a full (numerical) calculation of the Green's function including the gapped bands, without projection.\" Since the claimed correspondence between the measured vortex and the node winding is the paper's central result, this gap must be closed, either by a rigorous argument that the gapped bands do not alter the winding of the LDOS phase at the relevant scattering vector, or by a numerical evaluation of the unprojected Green's function for a concrete model (e.g., the NbSe2 model of Eq. (39)).","section":"III.B and Appendix C1, Eqs. (7), (13)-(16)"}],"minor_comments":[{"comment":"The statement that \"their total in the Brillouin zone must vanish\" is used to argue that all windings follow from the pairwise differences; a brief justification or reference for this topological constraint would be helpful.","section":"I. Introduction"},{"comment":"The full T-matrix expression is written for a local impurity potential, but the text does not explicitly state that the impurity is local until later; a sentence clarifying this would avoid ambiguity.","section":"III.A"},{"comment":"The same symbol θKi is used for the momentum-space phase in Eq. (11) and the real-space phase in Eq. (15); consider using a different notation (e.g., θKi^(q) and θKi^(r)) to avoid confusion.","section":"III.B, Eqs. (11) and (15)"},{"comment":"Two references appear as unresolved placeholders \"[?]\" in the footnote about the Su-Schrieffer-Heeger chain and Theorem 1; these need to be completed.","section":"III.E, footnote 53"},{"comment":"The caption refers to green and purple vectors, but the figure may be rendered in grayscale; please ensure the color coding is clear for all readers or label the vectors directly.","section":"Fig. 2 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a serious theoretical contribution with a clear and testable proposal. The main concern is the projection/gauge issue, which the authors themselves flag in Appendix C1; I believe it is addressable through a numerical calculation on the NbSe2 model or a sharper analytic argument. The condition error in Section III.C is local and fixable. I recommend major revision with a request for the additional verification rather than rejection, as the central idea is sound and the limitations are transparently discussed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere is my take on arXiv:2412.13042. The paper takes Dutreix et al.'s graphene QPI-dislocation mechanism and extends it to general chiral nodal systems, including nodal superconductors. The genuinely new pieces are the chirality-preference conditions on impurity and tip (χM, χV), the d-wave no-go statement, and the NbSe2 spin-polarized STM recipe. These are derived analytically, not fitted, and the authors are candid about limitations: Born estimate, anisotropic cone corrections, and approximate chiral symmetry.\n\nThe derivation is internally consistent. The LDOS decomposition in Eq. (19) into terms with winding ±Δθij is clean, and the conditions for observing the dislocation/vortex follow straightforwardly. I checked the logic from the chiral-basis Green's function to the winding conditions and did not find an error. The NbSe2 section is concrete, and the estimate in Appendix D that V G0 ≈ 0.3 is a useful sanity check.\n\nThe main soft spot is the projection to the cone band. The node winding is defined as a sum over bands, and the authors fix a gauge where all winding sits in the n=1 band, then project the Green's function to that band. They explicitly defer the full unprojected check in Appendix C1. This means the central claim—that the observed vortex equals the physical node winding difference—is conditional. The gapped band is separated by 2h, so the approximation is plausible, but given that the paper's entire purpose is to read a topological quantity from LDOS, I would want that check done before accepting the correspondence as established. I don't think this is fatal; it's a gap that a referee should ask to close.\n\nA smaller issue: the main-text condition Eq. (34a) is written for real χ, while the full complex form is only in Appendix A2. That's a presentation problem, not a physics one.\n\nNo code or data is shipped, but that's excusable for a theory paper; the analytic expressions are the main product.\n\nWho is this for? Experiment groups doing STM on candidate nodal superconductors, and theorists working on QPI. The paper deserves a serious peer review—the method is timely and the limitations are honestly stated. I would recommend engaging with it and asking for the unprojected numerical check as a major revision condition.","headline":"A useful generalization of QPI dislocation detection to nodal superconductors, with a concrete NbSe2 proposal, but the key projection step is deferred and needs verification.","tokens_in":26577,"tokens_out":5598,"would_cite":true,"duration_ms":51877,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The topological winding of superconducting nodes can be read from impurity-induced LDOS ripples, whose Fourier transform carries a vortex of strength equal to the winding difference.","keywords":["topological winding","nodal superconductor","chiral symmetry","local density of states","quasiparticle interference","scanning tunneling microscopy","NbSe2","phase vortex"],"falsifier":"Compute the full Bogoliubov-de Gennes Green's function without cone-band projection for the NbSe2 model, generate the impurity LDOS, and take its Fourier transform; if the vortex vorticity at any ΔKij differs from the winding difference predicted by the projected theory, the central claim collapses.","tokens_in":25267,"feed_emoji":"🌀","tokens_out":6796,"duration_ms":57341,"temperature":0.7,"pith_summary":"This paper claims that the topological winding of nodal points in chiral two-dimensional systems, including nodal superconductors, can be read directly from the local density of states measured by a scanning tunneling microscope near a single impurity. A difference in winding between two nodes produces a phase vortex in the Fourier transform of the impurity-induced LDOS, with vorticity equal to the winding difference, provided the impurity and the STM tip prefer opposite chiralities of the Bogoliubov wavefunctions. The paper derives explicit conditions on the impurity and tip, expressed through their chirality-preference ratios, and shows they are met in monolayer NbSe2 under an in-plane magnetic field with a magnetic impurity and a spin-polarized tip, making all 12 nodal windings detectable. If the claim holds, a purely local STM measurement can determine lower-dimensional topological invariants that bulk transport or edge-state probes cannot reach.","feed_headline":"Superconducting node windings show up as vortices in STM data","feed_subtitle":"A magnetic impurity plus spin-polarized tip can reveal the winding of all 12 nodes in monolayer NbSe2.","key_machinery":"The central object is the real-space phase function $\\theta_{K_i}(r) = \\arctan[(v_\\parallel/v_\\perp) \\tan(W_{K_i} \\theta_r + \\phi_i)]$, which inherits the node's winding number $W_{K_i}$ from the linearized dispersion and appears in the off-diagonal propagator $g_1$. Writing the scattering-induced LDOS as $\\mathrm{Re}[\\varrho(r)e^{i\\varphi_r}]$ with $\\varphi_r = \\Delta K_{ij}\\cdot r + \\mathrm{Arg}[I_0 + I_\\Delta e^{i\\Delta\\theta} + I_{-\\Delta} e^{-i\\Delta\\theta}]$, the winding of $\\mathrm{Arg}[\\cdot]$ around the impurity yields the dislocation charge $W_{\\delta\\rho} = \\pm\\Delta W$. The conditions for a non-zero $W_{\\delta\\rho}$ reduce to two chirality-preference ratios: $\\chi_M = M_A/M_B$ for the tip and $\\chi_V = V_A/V_B$ for the impurity, with the requirement $|h_1(\\chi_M+\\chi_V)| > |h_0(\\chi_M\\chi_V+1)|$ and $\\chi_M \\neq \\pm\\chi_V$.","core_discovery":"In a chiral Hamiltonian written in a basis of A and B eigenstates, the complex phase θn(k) of the off-diagonal Q-matrix carries the winding. After linearizing around a node and projecting to the cone band, the real-space Green's function g1(r) acquires the same phase as a function of the real-space polar angle θr, so the cone's winding WKi appears in the phase of the propagator. Scattering between two nodes i and j by an impurity then contributes LDOS terms cos(ΔKij·r ± Δθij(θr)) to the wavefront, and when the tip and impurity have different chirality preferences, the phase Δθij winds ΔW = WKi - WKj times, creating a dislocation of strength ±ΔW in real space. Crucially, in the Fourier transform of the LDOS this becomes a single vortex at ΔKij whose vorticity equals the total enclosed dislocation charge, a quantity that is robust to finite field of view. The paper's application to the twelve-node nodal superconductor proposed for monolayer NbSe2 under an in-plane field shows that a magnetic impurity combined with a spin-polarized STM tip satisfies the winding conditions, so each node's winding is in principle measurable.","pith_inferences":["Beyond the paper: a decisive numerical test would be a full, unprojected Green's function calculation of the impurity LDOS in NbSe2; if the vortex vorticity at each ΔKij still equals the total winding difference, the projection to the cone band is validated.","Beyond the paper: the same scattering mechanism could be used on the surfaces of three-dimensional Weyl semimetals or nodal-line superconductors, where the projected surface LDOS may carry winding signatures from the bulk nodes.","Beyond the paper: the requirement that tip and impurity prefer opposite chiralities suggests that functionalized or spin-polarized STM tips could act as chirality filters, potentially enabling detection of higher-winding Dirac points or multi-band nodes with |W| > 1.","Beyond the paper: if confirmed, the technique gives a local, real-space probe of topological invariants that does not require edge states or global Berry curvature cancellation, which could be applied to disordered or small samples."],"forward_implications":["A single STM measurement, not quantized transport, can extract the winding numbers of superconducting nodes, which are otherwise inaccessible experimentally.","In a nodal superconductor such as NbSe2, the measured winding differences can discriminate pairing symmetries, since only certain pairing functions yield non-zero windings at the nodes.","The method extends beyond superconductors to any chiral system with linearly dispersing nodes, including graphene-like semimetals, and to nodal lines in three dimensions.","Momentum-space vortices remain observable even when real-space dislocations are masked by other defects, because the vorticity at ΔKij accumulates all dislocation charges inside the field of view.","The explicit chirality-preference conditions provide a recipe to choose an impurity and tip (or functionalized tip) for any candidate material."],"supporting_citations":[{"why":"Graphene experiment detecting Dirac cone winding via wavefront dislocations in Friedel oscillations; the method this paper generalizes to chiral systems and superconductors.","marker":"[22]"},{"why":"Proposal of magnetic-field-driven nodal topological superconductivity in monolayer transition metal dichalcogenides; supplies the physical context for NbSe2.","marker":"[37]"},{"why":"Provides the BdG Hamiltonian, basis and unitary transformation used for the NbSe2 node winding and scattering calculations.","marker":"[39]"},{"why":"Derives the Green's function expression for a chiral system used in Eq. (8), though in one dimension and without nodes.","marker":"[41]"},{"why":"Defines dislocations in wave trains, the framework behind writing the LDOS as a complex field with phase winding.","marker":"[43]"},{"why":"Shows a dislocation with constant charge around a vacancy in graphene, used as a contrasting case in the robustness discussion.","marker":"[47]"},{"why":"Supplies the Yu-Shiba-Rusinov impurity strength estimate for Fe in NbSe2 used to justify the Born approximation.","marker":"[56]"},{"why":"Reports experimental evidence of nodal and nematic superconducting phases in NbSe2 monolayers, supporting the twelve-node picture.","marker":"[49]"}],"fun_headline_variants":["STM vortices expose node winding in superconductors","Spin-polarized STM reads node winding from impurity vortices","Winding topology of nodes emerges in STM impurity patterns","Magnetic impurity plus spin tip reveals NbSe2 node windings"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that after a gauge choice putting all winding into the cone band, no gapped band carries any of the node's winding, so the projected propagator measures the true node winding; the authors state that a full numerical Green's function check without projection is still needed.","fun_headline_variants_meta":{"raw":{"variants":["STM vortices expose node winding in superconductors","Spin-polarized STM reads node winding from impurity vortices","Winding topology of nodes emerges in STM impurity patterns","Magnetic impurity plus spin tip reveals NbSe2 node windings"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000544,"raw_usage":{"total_tokens":2688,"prompt_tokens":1114,"completion_tokens":1574,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":730,"completion_tokens_details":{"reasoning_tokens":1506}},"tokens_in":730,"tokens_out":1574,"duration_ms":13156,"temperature":1.0,"reasoning_tokens":1506,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T13:32:42.801164+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full Bogoliubov-de Gennes Green's function without cone-band projection for the NbSe2 model, generate the impurity LDOS, and take its Fourier transform; if the vortex vorticity at any ΔKij differs from the winding difference predicted by the projected theory, the central claim collapses.","supporting_citations":[{"cited_title":"Volovik, Zeros in the fermion spectrum in superfluid systems as diabolical points, JETP Lett 46, 98 (1987)","cited_arxiv_id":null,"evidence_quote":"Graphene experiment detecting Dirac cone winding via wavefront dislocations in Friedel oscillations; the method this paper generalizes to chiral systems and superconductors."},{"cited_title":"Zhang, J","cited_arxiv_id":null,"evidence_quote":"Proposal of magnetic-field-driven nodal topological superconductivity in monolayer transition metal dichalcogenides; supplies the physical context for NbSe2."},{"cited_title":"Spin wave vortex as topological probe of magnetic texture","cited_arxiv_id":"2405.03566","evidence_quote":"Provides the BdG Hamiltonian, basis and unitary transformation used for the NbSe2 node winding and scattering calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the Green's function expression for a chiral system used in Eq. (8), though in one dimension and without nodes."},{"cited_title":"Uldemolins, A","cited_arxiv_id":null,"evidence_quote":"Shows a dislocation with constant charge around a vacancy in graphene, used as a contrasting case in the robustness discussion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Yu-Shiba-Rusinov impurity strength estimate for Fe in NbSe2 used to justify the Born approximation."},{"cited_title":"Shaffer, J","cited_arxiv_id":null,"evidence_quote":"Reports experimental evidence of nodal and nematic superconducting phases in NbSe2 monolayers, supporting the twelve-node picture."}],"review_version":1}