{"id":"587ddc64-ddd5-4a41-b167-29d0b3633967","arxiv_id":"1908.02955","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A reduced response function analysis of quasiparticle interference shows that the number of hot arcs in a spin-resolved channel equals twice the topological charge of a Dirac point.","lead":"This paper develops a reduced response function framework for analyzing quasiparticle interference in scanning tunneling microscopy. It shows how counting bright arcs in spin-resolved interference patterns can reveal the topological charge of a Dirac point.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Admitted false GJDOS features at ±2K overlap the intervalley arc indicators, so the claimed exclusion of finite-lifetime discrepancies fails exactly where the topological-number indicator is read.","rationale":"Good faith: paper develops RRF, derives Eq.16, gives numerical evidence that S and J share singularities in ideal convex cases, and shows arc patterns in Bi2Te3, BiTeI, and graphene. The topological rule 2|n| arcs for ideal H_n is analytically plausible and the numerics support it for S in several models. Credit: parameter-free derivations and comparisons with experimental QPI patterns. However, the central claim has two linked dependencies: (1) Eq.14 is asserted after 'tedious derivation' and not shown, so the singular equivalence is not fully demonstrated; (2) the finite-lifetime passage from S to J is justified only by an exclusion that the text itself undermines at ±2K. I focus on (2) because it directly threatens the applicability of the topological indicator to real QPI data, which is the advertised payoff. The reader's weakest assumption identifies the same finite-lifetime exclusion; my concern makes it concrete by pointing to the paper's own admission of false J hot spots at the indicator momenta. I do not view this as fatal: the materials-section results are computed with S, and a revised statement restricting the indicator to the stable part of the response, or proving that false J features do not shift arc counts, would preserve the core contribution. Hence the verdict stays conditional.","tokens_in":13814,"tokens_out":10389,"duration_ms":116714,"concrete_test":"Compute S00(p,ω) and J00(p,ω) for monolayer graphene at ω=0.3 with a finite quasiparticle lifetime η in the range used for the paper's figures, and compare the maps in windows around p=2K, p=−2K, and p=2(K−K′). Determine whether J00 has local maxima ('false hot spots') at the same momentum points where S00 shows the intervalley arcs, and whether the number of disconnected bright segments in J00 matches 2N=2. Repeat for bilayer and trilayer (N=2,3). If J00's false features overlap the arcs or change the segment count, the finite-lifetime exclusion fails at the indicator momenta, and the central claim that GJDOS is justified for studying QPI and for topological-number extraction requires revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing concern is the finite-lifetime exclusion that separates GJDOS from the full FT-LDOS. The paper asserts in Sec. III (intro) that discrepancies between S and J are limited to recognizable patterns such as a p=0 hot spot or asterisk and can be excluded 'at first hand' when using J. This assertion is not proved, and it is contradicted by the paper's own graphene section: Sec. IV C closes with 'joint density of states give false hot spots at p=0, ±2K, ±2K, and ±2(K−K′).' These momenta are exactly where the intervalley topological indicators of Fig. 9 are located: the 2N disconnected hot arcs arise from intervalley scattering at ±2K, ±2K′, and ±2(K−K′). Thus the admitted false features overlap the readout region of the central indicator. If the stated justification is followed and J is used directly, the false hot spots will be counted as part of or will obscure the hot arcs. The numerical demonstrations of the materials sections use S (the real part), not J, which sidesteps the problem, but the central claim as stated is that J is justified for studying QPI and that the arc-counting rule is an indicator of the topological number. The unproven and partially contradicted exclusion is therefore load-bearing; without it, the rule is valid for ideal zero-lifetime models but not established for experimental finite-lifetime QPI.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a reduced response function (RRF) framework for quasiparticle interference, decomposing the Fourier-transformed LDOS into real and imaginary parts S and J (the generalized joint density of states, GJDOS). It claims that S and J share the same singularities in ideal zero-lifetime systems and that, under finite quasiparticle lifetime, discrepancies between them are limited to recognizable false features that can be 'excluded at first hand,' thereby justifying the use of GJDOS for QPI analysis. From an analytical expression for GJDOS, the paper derives topological-number indicators: for the scattering/probe channel alpha = beta = x, a Dirac point of topological charge n produces 2|n| hot arcs, while the alpha = beta = 0 channel gives an odd-even distinction. The indicators are demonstrated numerically on Bi2Te3, BiTeI, and ABC-stacked N-layer graphene.","tokens_in":1633,"tokens_out":3103,"duration_ms":97372,"significance":"The paper offers a clean decomposition of the FT-LDOS response into real and imaginary parts, a transparent spin-coherent-factor algebra (Sec. II B), and a closed-form GJDOS expression (Eq. (16)) that enables a geometric reading of QPI patterns. The proposed topological-number indicators for the ideal Hamiltonian Eq. (20) are elegant and produce falsifiable predictions: 2|n| hot arcs for the x-x channel and an odd-even selection for the 0-0 channel. The numerical simulations for Bi2Te3, BiTeI, and ABC-stacked graphene are internally consistent and show that global patterns survive complicated constant-energy contours. If the derivation behind Eq. (14) is supplied and the status of J relative to S is clarified, the framework would be a useful contribution to the QPI literature. The current manuscript, however, does not establish the central claim that GJDOS is a generally justified substitute for FT-LDOS.","major_comments":[{"comment":"The central justification for using J (GJDOS) in place of the full FT-LDOS rests on the assertion in the Sec. III introduction that finite-lifetime discrepancies between S and J are limited to recognizable patterns (a p=0 hot spot or an asterisk-like feature) that can be excluded 'at first hand' when using J. No proof of this classification is provided. The assertion is directly contradicted by the paper's own closing statement of Sec. IV C, which admits that the joint density of states gives false hot spots at p=0, plus or minus 2K, plus or minus 2K prime, and plus or minus 2(K - K prime). These momenta are precisely the regions where the intervalley topological indicator is read: the 2N hot arcs in Fig. 9 arise from intervalley scattering at those same wave vectors. Thus the admitted false features overlap the readout region of the central indicator. Since the numerical demonstrations of the indicators in Sec. IV are all performed with S (the real part) rather than J, the paper does not actually demonstrate that J can be used directly for topological-number inference in finite-lifetime systems. This is a load-bearing gap: without a proof of the exclusion or a restriction of the indicator claim to S, the central claim that GJDOS is a justified tool for QPI analysis is not established.","section":"Sec. III (intro) and Sec. IV C"},{"comment":"The asymptotic formula Eq. (14) is introduced with the phrase 'after a tedious derivation we can show,' but no derivation is supplied in the main text or in an appendix. This formula is the mathematical foundation for the statement that the real and imaginary parts of the reduced response function share the same singularities and for the subsequent joint-curvature analysis used throughout the paper, including the classification of higher-order singularities in Fig. 2 and the stability discussion in Sec. II D. Without a presented derivation or a reference to a complete proof, the central analytic result of the paper is not independently verifiable from the manuscript. The authors should include the derivation in an appendix or provide a detailed outline of the stationary-phase or contour-integral argument that leads to Eq. (14).","section":"Sec. II C, Eq. (14)"}],"minor_comments":[{"comment":"The sentence 'It is justified that the generalized joint density of states, which is the imaginary part of RRF, for studying QPI' lacks a main verb; it should be rewritten, for example, as 'the generalized joint density of states, which is the imaginary part of the RRF, is justified for studying QPI.'","section":"Abstract"},{"comment":"The text states that S is 'an autocorrelation of A(k,omega) and A(k,omega)'; from Eq. (10), S is a cross-correlation of A and B (A B + B A), not an autocorrelation of A with itself. The subsequent discussion of the sign change of B confirms that this is a typographical error that should be corrected.","section":"Sec. II D"},{"comment":"The sentence 'distinct 2N pieces of disconnected hot arcs appear for intravalley scattering' should read 'intervalley scattering,' consistent with the figure caption and the surrounding discussion that places the hot arcs at wave vectors such as plus or minus 2K and plus or minus 2(K - K prime).","section":"Sec. IV C"},{"comment":"The list of false hot spots 'p=0, plus or minus 2K, plus or minus 2K, and plus or minus 2(K - K prime)' appears to contain a duplicated 'plus or minus 2K'; presumably the second one should be 'plus or minus 2K prime.'","section":"Sec. IV C"},{"comment":"The caption 'Behavior of dot-E_{k_omega(t0)}+p' is unclear notation; it likely refers to the energy at k_omega(t0)+p and should be typeset accordingly for readability.","section":"Fig. 1 caption"},{"comment":"The 'Positive-negative indicator' is defined by rotating the scattering and probe channels 'a little,' but the required rotation angle (presumably pi/(2n)) is not specified, which makes the prescription ambiguous.","section":"Sec. III B"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal scope and the reference list is appropriate. The main concern is the gap between the paper's stated justification of GJDOS and its own numerical demonstrations, which use S rather than J. I recommend that the revision address this directly, possibly by separating two claims: (i) S and J share singularities in ideal systems (which requires the missing derivation for Eq. (14)), and (ii) J may be used in practice with recognized false features, while the topological indicators are best read from S or from J only after a controlled analysis. The numerical demonstrations already support the indicator for S, so the paper's core contribution to topological inference could be preserved with a more modest claim about J."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth a referee's time, but the central justification for reading QPI patterns with GJDOS has a hole, and the authors themselves dig it.\n\nWhat's new: the reduced response function R = S + iJ, the joint-curvature classification of singular orders (Eq. 15), and the 2|n| hot-arc rule for the α=β=x channel. That rule—counting disconnected hot arcs as a readout of the winding number of a Dirac point—is concrete, checkable, and a genuinely useful diagnostic. The spin coherent factor algebra in Sec. II B checks out, and the numerical demonstrations for Bi2Te3, BiTeI, and ABC graphene are internally consistent and reproduce the expected patterns.\n\nThe soft spot is exactly where the reader's stress-test lands. The paper argues (Sec. III intro) that discrepancies between S and J under finite lifetime are limited to recognizable patterns like a p=0 hot spot or asterisk and can be 'excluded at first hand' when using J. Then in Sec. IV C, the authors write that 'joint density of states give false hot spots at p=0, ±2K, ±2K, and ±2(K−K′).' Those are the momenta where their own intervalley topological indicators live—the 2N hot arcs of Fig. 9. So the admitted false features overlap the readout region of the central indicator. The numerical figures use S, not J, which sidesteps the problem, but the paper's stated claim is that J is justified for QPI. Without a proof of the exclusion, or a qualification that the rule applies to ideal zero-lifetime models only, the arc-counting indicator is not established for experimental finite-lifetime data. This is load-bearing, not a cosmetic issue.\n\nThe math itself is mostly sound; Eq. (14) is stated as 'after a tedious derivation' without proof, which is a legitimate referee request, but the stress-test concern is more serious. The citation to Ref. 32 for stationary phase is appropriate, and the paper is honest about the false features—which makes the oversight more surprising.\n\nBottom line: this is a clever framework with a plausible central rule, but the paper needs to either prove the exclusion or restrict the claim before the rule is ready for experimental use. Send it out; a good referee will catch this and the revision can strengthen it.","headline":"Clever framework and a concrete arc-counting rule, but the paper's own graphene section admits false GJDOS hot spots at exactly the momenta where the rule is read, so the central justification for using J is not established.","tokens_in":14632,"tokens_out":3807,"would_cite":false,"duration_ms":35405,"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":"For a Dirac point of topological charge $n$, an $x$-scattering and $x$-probe quasiparticle interference image shows exactly $2|n|$ disconnected bright arcs, making arc counts topological-number indicators.","keywords":["quasiparticle interference","reduced response function","generalized joint density of states","scanning tunneling spectroscopy","topological charge","Dirac point","spin-momentum locking","ABC-stacked graphene"],"falsifier":"Take the model $H_n(\\mathbf k)=k^n(\\cos n\\theta\\,\\sigma_x+\\sin n\\theta\\,\\sigma_y)$ with a small warping term added to flatten a segment of a constant-energy contour, and compute both the real part $S_{xx}$ and the imaginary part $J_{xx}$ at finite lifetime; if any bright feature of $J_{xx}$ that has no counterpart in $S_{xx}$ lies at the same momentum as one of the expected arcs, the count is ambiguous. Experimentally, a spin-resolved FT-STS measurement of ABC-stacked trilayer graphene counting the bright arcs in the intravalley and intervalley regions would settle whether the predicted $2N$ arc counts appear exactly as stated.","tokens_in":13633,"feed_emoji":"🔬","tokens_out":16725,"duration_ms":164637,"temperature":0.7,"pith_summary":"Scanning tunneling microscopy sees quasiparticle interference (QPI) as standing-wave ripples around impurities, and its Fourier transform gives patterns that depend on the band structure. This paper develops a reduced response function whose imaginary part, the generalized joint density of states, is shown to share all singular features of the full Fourier-transformed local density of states except a few recognizable artifacts. That justifies using this simpler object to read QPI patterns directly from constant-energy contours. For gapless systems the paper derives global indicators of topology: in the channel where both scattering and probing are $x$-polarized, a Dirac point of topological charge $n$ produces exactly $2|n|$ disconnected bright arcs, so counting arcs reads the topological number. The same indicators are shown numerically to survive complicated band geometry in Bi$_2$Te$_3$, BiTeI, and ABC-stacked $N$-layer graphene.","feed_headline":"Count the arcs to read a Dirac point's topological charge","feed_subtitle":"A reduced response function shows 2|n| disconnected bright arcs for a Dirac point of charge n in spin-resolved FT-STS.","key_machinery":"The load-bearing object is the reduced response function $R_{\\alpha\\beta}(\\mathbf p,\\omega)=S_{\\alpha\\beta}(\\mathbf p,\\omega)+iJ_{\\alpha\\beta}(\\mathbf p,\\omega)$, assembled from an autocorrelation of spectral functions $A_s(\\mathbf k,\\omega)=\\delta(\\omega-E^s_{\\mathbf k})$ and their Hilbert partners $B_s(\\mathbf k,\\omega)$. Its imaginary part $J_{\\alpha\\beta}$, the generalized joint density of states, is the practical imaging tool: it integrates the spin coherent factor $F^{ss'}_{\\alpha\\beta}(\\mathbf k+\\mathbf p,\\mathbf k)$ over pairs of momenta on constant-energy contours that are separated by $\\mathbf p$. The singularity condition $v^{s'}_{\\mathbf k_0+\\mathbf p}\\times v^s_{\\mathbf k_0}=0$ reduces QPI pattern formation to contour geometry, and the channel-dependent structure of the coherent factor encodes how each scattering and probing channel suppresses or restores those singularities: for $\\alpha=\\beta=x$ the factor reflects spins across the $x$-axis, which selects the directions that form $2|n|$ arcs. Because the real part is the Hilbert transform of the imaginary part, the two have the same singular behavior, which is what licenses reading FT-LDOS from the simpler $J$ once its recognizable artifacts are discarded.","core_discovery":"The paper's central claim is that for a two-band Hamiltonian $H(\\mathbf k)=E_0(\\mathbf k)+\\mathbf d(\\mathbf k)\\cdot\\boldsymbol\\sigma$, the reduced response function $R_{\\alpha\\beta}(\\mathbf p,\\omega)=S_{\\alpha\\beta}(\\mathbf p,\\omega)+iJ_{\\alpha\\beta}(\\mathbf p,\\omega)$ carries the full singular information of quasiparticle interference: its real part $S$ and imaginary part $J$ share the same singularities, and $J$ is the generalized joint density of states, an autocorrelation of spectral functions along constant-energy contours weighted by a spin coherent factor $F^{ss'}_{\\alpha\\beta}$. The singularities occur when a pair of points on two contours linked by $\\mathbf p$ have group velocities satisfying $v^{s'}_{\\mathbf k_0+\\mathbf p}\\times v^s_{\\mathbf k_0}=0$, i.e. parallel or antiparallel velocities, with an inverse-square-root divergence that becomes stronger where the joint curvature vanishes. For an ideal Dirac point $H_n(\\mathbf k)=k^n(\\cos n\\theta\\,\\sigma_x+\\sin n\\theta\\,\\sigma_y)$ of topological charge $n$, the paper proves that the channel $\\alpha=\\beta=x$ produces exactly $2|n|$ disconnected hot arcs, while the charge channel $\\alpha=\\beta=0$ suppresses backscattering for odd $n$; the number of arcs is therefore a global topological-number indicator. Numerical evaluation for Bi$_2$Te$_3$, BiTeI, and ABC-stacked $N$-layer graphene shows that the indicator remains legible when constant-energy contours are nonconvex and when local features such as cusp-scattering hot spots and false joint-density peaks distort the images.","pith_inferences":["Inference: because the mechanism uses only contour geometry and a channel-dependent coherent factor, the same $2|n|$ arc-counting should extend to multi-band crossings, nodal lines, or non-Hermitian band touchings as long as generalized constant-energy contours and an effective spin/pseudospin factor exist; the paper's two-band restriction is a convenience, not a prerequisite.","Inference: the 'artifacts are excludable at first hand' rule could be turned into a quantitative protocol by comparing $S_{\\alpha\\beta}$ and $J_{\\alpha\\beta}$ over a range of finite lifetimes; any feature of $J$ that disappears from $S$ as a contour flattens is spurious, so the topological indicator could be automated rather than judged by eye.","Inference: the sign of the topological charge is read from the direction the arcs rotate when the $x$-channel is rotated slightly, which predicts that reversing the chirality of a sample reverses the rotation direction in a spin-resolved FT-STS experiment; this is a direct consequence the paper does not itself test."],"forward_implications":["A spin-resolved FT-STS experiment with an $x$-polarized impurity and $x$-sensitive probe can read the topological charge of a Dirac point by counting bright arcs: $2|n|$ arcs means charge $n$.","For the graphene family, two valleys with opposite charges $\\pm N$ produce $2N$ disconnected hot arcs in intervalley QPI, and the charge channel distinguishes whether the charge difference is odd or even.","The generalized joint density of states is a legitimate stand-in for full FT-LDOS in QPI analysis, provided its known artifacts, the $\\mathbf p=0$ hot spot and asterisk-like features from flat or near-nesting contours, are recognized and discarded.","Zero joint curvature at cusp-to-cusp scattering produces higher-order singularities whose hotter spots should dominate experimental QPI images of Bi$_2$Te$_3$ when the constant-energy contour becomes nonconvex; the paper reproduces the measured FT-STS patterns.","The topological arc indicators survive nonconvex contours and complicated local geometry, so they need not wait for idealized band structures to be useful in real materials."],"supporting_citations":[{"why":"It introduced the joint density of states as the origin of hot spots in QPI, the picture the reduced response function generalizes.","marker":"[4]"},{"why":"It generalized JDOS to include the internal spin structure of quasiparticles, the quantity whose validity the paper justifies.","marker":"[5]"},{"why":"It cautioned that GJDOS can give false QPI features, which the paper explains as recognizable artifacts.","marker":"[22]"},{"why":"It showed discrepancies between GJDOS and full FT-LDOS, motivating the reduced response function construction.","marker":"[31]"},{"why":"It defined the topological charge of Dirac points used to formulate the arc-counting indicator.","marker":"[33]"},{"why":"It supplied the effective Hamiltonians for ABC-stacked N-layer graphene, the testbed with Dirac points of charges +N and -N.","marker":"[34]"},{"why":"It supplied the hexagonal-warping Hamiltonian for Bi2Te3 surface states used to test whether the indicators survive nonconvex contours.","marker":"[19]"},{"why":"It supplied the low-energy Bi2Te3 model with topological number Q=1, the base model for the numerical QPI patterns.","marker":"[38]"},{"why":"It provided experimental FT-STS patterns of Bi2Te3 against which the numerical QPI images are checked.","marker":"[20]"}],"fun_headline_variants":["Count arcs to read Dirac point charge","2|n| arcs = topological charge in QPI","Quasiparticle arcs reveal topological numbers","Arc count tells Dirac point's charge","Topological charge by arc number in QPI"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The arc-counting indicator rests on the assumption, stated rather than proved, that with finite quasiparticle lifetime the only differences between the easily computed imaginary part and the real QPI signal are recognizable artifacts such as a $\\mathbf p=0$ hot spot or an asterisk-like pattern, and that these never appear at the momenta where the $2|n|$ arcs are being counted.","fun_headline_variants_meta":{"raw":{"variants":["Count arcs to read Dirac point charge","2|n| arcs = topological charge in QPI","Quasiparticle arcs reveal topological numbers","Arc count tells Dirac point's charge","Topological charge by arc number in QPI"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00043,"raw_usage":{"total_tokens":2254,"prompt_tokens":1057,"completion_tokens":1197,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":673,"completion_tokens_details":{"reasoning_tokens":1129}},"tokens_in":673,"tokens_out":1197,"duration_ms":12779,"temperature":1.0,"reasoning_tokens":1129,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:28:54.067955+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the model $H_n(\\mathbf k)=k^n(\\cos n\\theta\\,\\sigma_x+\\sin n\\theta\\,\\sigma_y)$ with a small warping term added to flatten a segment of a constant-energy contour, and compute both the real part $S_{xx}$ and the imaginary part $J_{xx}$ at finite lifetime; if any bright feature of $J_{xx}$ that has no counterpart in $S_{xx}$ lies at the same momentum as one of the expected arcs, the count is ambiguous. Experimentally, a spin-resolved FT-STS measurement of ABC-stacked trilayer graphene counting the bright arcs in the intravalley and intervalley regions would settle whether the predicted $2N$ arc counts appear exactly as stated.","supporting_citations":[{"cited_title":"Capriotti , author D","cited_arxiv_id":null,"evidence_quote":"It introduced the joint density of states as the origin of hot spots in QPI, the picture the reduced response function generalizes."},{"cited_title":"\\ Wang \\ and\\ author D.-H","cited_arxiv_id":null,"evidence_quote":"It generalized JDOS to include the internal spin structure of quasiparticles, the quantity whose validity the paper justifies."},{"cited_title":"Kohsaka , author T","cited_arxiv_id":null,"evidence_quote":"It cautioned that GJDOS can give false QPI features, which the paper explains as recognizable artifacts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It showed discrepancies between GJDOS and full FT-LDOS, motivating the reduced response function construction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It defined the topological charge of Dirac points used to formulate the arc-counting indicator."},{"cited_title":"Min \\ and\\ author A","cited_arxiv_id":null,"evidence_quote":"It supplied the effective Hamiltonians for ABC-stacked N-layer graphene, the testbed with Dirac points of charges +N and -N."},{"cited_title":"\\ Lee , author C","cited_arxiv_id":null,"evidence_quote":"It supplied the hexagonal-warping Hamiltonian for Bi2Te3 surface states used to test whether the indicators survive nonconvex contours."},{"cited_title":"Beidenkopf , author P","cited_arxiv_id":null,"evidence_quote":"It provided experimental FT-STS patterns of Bi2Te3 against which the numerical QPI images are checked."}],"review_version":1}