{"id":"eb63d5cd-70ef-4d40-92f0-d9ef4f2a9df5","arxiv_id":"2512.16211","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Periodic particle-hole-symmetric gap modulations in bulk FeSe are attributed to pair-breaking scattering interference from magnetic impurities, not pair density wave order.","lead":"This experiment finds small periodic wiggles in the superconducting gap of bulk FeSe near magnetic impurities, with a spacing that matches the theory of pair-breaking scattering interference. The work suggests such gap ripples can be explained without invoking exotic pair density wave order.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gap maps may be an artifact of a spatially varying Nb-tip gap; the in-phase Δ+ / Δ− pattern is exactly the common-mode signature, and the phase-referenced QPI cannot exclude it.","rationale":"The reader's weakest assumption is the right one: the primary quantitative output is the two-gap fitting pipeline, and no error bars or controls are given. I refined it by identifying a concrete artifact mechanism: a spatially varying Nb-tip gap or energy offset produces a common-mode shift of both coherence peaks, which would appear exactly as the observed in-phase Δ_+/Δ− modulation and would survive the symmetrized/antisymmetrized decomposition. The phase-referenced QPI is not an independent rescue because any particle-hole-symmetric modulation yields the same sign pattern. The paper has genuine strengths—the superconductive tip, YSR identification, and the Josephson modulation—but these are not sufficient to rule out the artifact. A targeted test, refitting with a free spatial Δ_Nb, would settle the concern. Pending that, the CONDITIONAL verdict is appropriate and unchanged.","tokens_in":10100,"tokens_out":6993,"duration_ms":75302,"concrete_test":"Refit the raw dI/dV stack from Fig. 3 with Δ_Nb allowed to vary pixel-by-pixel as a free parameter, while keeping the FeSe two-gap model fixed. If the resulting Δ_Nb(r) map shows a Fourier peak at 1.3 nm^-1 comparable to the originally reported Δ_α± modulation, and the FeSe gap maps lose the Q0 peak, the artifact hypothesis is confirmed. If Δ_Nb(r) is spatially uniform and the Q0 peak persists in the FeSe gap maps, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central evidence for PBSI is the Q0 modulation in the two-gap fit outputs (Fig. 3(a–h)). The fitting pipeline (Supplemental Sections II–III, ref. [36]) deconvolves the Nb-tip DOS and then extracts Δ_α±(r) and Δ_ε±(r) with a two-gap model. The reported in-phase modulation of Δ_α+ and Δ_α− (enhanced in the symmetrized sum, suppressed in the difference) is precisely the signature of a spatially varying common-mode energy shift. A spatially varying Nb superconducting gap Δ_Nb(r), or an uncontrolled energy offset, would shift both coherence peaks in the same direction and survive the symmetrized/antisymmetrized analysis. No error bars, control measurements, or spatial maps of the fitted Δ_Nb are provided. The phase-referenced QPI in Fig. 4 is not a clean independent arbiter: any particle-hole-symmetric modulation (including a tip-induced common shift) produces positive g~(q, +E, −E), and the predicted sign pattern follows from the coherence-peak structure rather than from PBSI uniquely. Thus the claim rests on the fidelity of the deconvolution, which has not been demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports STM/SJTM measurements at 300 mK on bulk FeSe using superconducting Nb tips. The authors identify subsurface magnetic impurities through spatially dispersing in-gap Yu-Shiba-Rusinov states and reduced Josephson current. Around these impurities, they extract spatial maps of the two superconducting gaps Δ_{α±}(r) and Δ_{ε±}(r) from deconvolved dI/dV spectra and find particle-hole symmetric modulations with wavevectors Q0 ≈ 1.3 nm⁻¹ and Q1 that match intra-pocket scattering vectors of the α and ε Fermi pockets. They also report a faint modulation of the maximum Josephson current at Q0 and propose a phase-referenced QPI analysis whose sign pattern is consistent with the s=1 pair-breaking scattering interference (PBSI) formalism of ref. [28]. The central claim is that these observations establish PBSI as a viable origin of superconducting gap modulations in a superconductor without preexisting density-wave order.","tokens_in":10415,"tokens_out":3573,"duration_ms":42830,"significance":"If the central claim holds, the paper makes an important contribution by providing experimental support for a recently proposed alternative to pair-density-wave interpretations of periodic gap modulations. The use of superconducting tips for improved energy resolution, the identification of subsurface magnetic scatterers, and the proposal of a phase-referenced QPI diagnostic are valuable methodological steps. The paper explicitly connects its observations to falsifiable predictions from external theory (refs. [28,30]) and includes a raw LDOS coherence-peak trace (Fig. 3(j)) that shows a real 0.12 meV oscillation. The significance is therefore potentially high for the STM community studying gap modulations, provided the fit-derived gap maps are shown to be robust against instrumental artifacts.","major_comments":[{"comment":"The central evidence for PBSI is the Q0 modulation in the gap maps, but these maps are outputs of a deconvolution and two-gap fitting pipeline (Supplemental Sections II–III, ref. [36]). No error bars or uncertainties are given for Δ_{α±}(r), no statistical significance is assigned to the FFT peak at Q0, and no control is shown that the fitting procedure itself does not produce a spurious 1.3 nm⁻¹ modulation. The raw LDOS trace in Fig. 3(j) is helpful but is one line cut; it does not establish statistical significance over the full field of view. The authors should provide quantitative controls: for example, FFT of fit residuals, simulated maps with a spatially uniform gap passed through the same pipeline, or a comparison of the extracted Δ_{α±}(r) with an independent coherence-peak extraction method. Without such controls, a tip-induced common-mode artifact cannot be excluded.","section":"§3, Fig. 3(a–j)"},{"comment":"The phase-referenced QPI signal g̃(q, E1, E2) is presented as an independent and direct probe of PBSI. However, the sign pattern predicted in Fig. 1(g) for s=1 is also what one would expect from any particle-hole symmetric modulation of the coherence peaks, including a spatially varying common-mode energy shift or a spatially varying Nb tip gap. Since the g(r,E) images used in Fig. 4 are the same dI/dV maps from which the gap maps are derived, the phase-referenced QPI is an internal consistency check, not an independent arbiter. To strengthen the claim, the authors should show that the observed sign pattern is inconsistent with a common-mode shift, for example by mapping the fitted Nb gap or an energy offset and demonstrating that it does not oscillate at Q0, or by repeating the measurement with a different tip condition.","section":"§4, Fig. 4"},{"comment":"The Josephson current modulation is admitted in the text to be only 'faintly observed', with crest lines placed by hand. No FFT or quantitative line profile of I_J(r) is given, and the claimed period and phase relative to the gap modulation are not extracted. If this measurement is used as supporting evidence for simultaneous Δ(r) and superfluid-density modulation, it needs quantitative analysis with errors; otherwise it should be clearly labeled as suggestive and not load-bearing.","section":"Josephson imaging, Fig. 3(k–l)"}],"minor_comments":[{"comment":"The notation Δ_{0,2±} and Δ_{1,2±} is confusing: the comma-separated subscripts are easily misread. Consistent notation such as Δ_{α±} and Δ_{ε±} would improve readability.","section":"Notation throughout"},{"comment":"In Eq. (1), the definition 'Δ_{),(}=Δ_0±V' appears to contain a typographical artifact. Please clarify the subscript and define all quantities, including the delta-function broadening used in the plotted LDOS.","section":"Eq. (1)"},{"comment":"The statement that 'no spin or charge order exists in bulk FeSe' is categorical; it would benefit from a qualification and a citation, since magnetic fluctuations and possible short-range order have been discussed in the FeSe literature.","section":"§2, first paragraph"},{"comment":"The sign symmetry argument for reducing the six combinations to three is correct, but the color scale and sign convention of g̃ should be stated explicitly in the figure or caption to avoid ambiguity.","section":"Fig. 4(f)"},{"comment":"Ref. [36] is a previous work by the same group and is used for the fitting pipeline. The paper would be stronger if the Supplemental Material cross-referenced the specific equations or procedures from ref. [36] in more detail, since the current text simply says 'Supplemental Material Sections II, III'.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This is a technically sophisticated and timely study, and the proposed phase-referenced QPI approach is a useful contribution. The main risk is that the fit-derived gap maps are not sufficiently validated against instrumental or analysis artifacts, and the phase-referenced QPI does not fully exclude a common-mode coherence-peak shift. The authors should be asked for quantitative controls and error analysis. If those are provided, the paper could become suitable for publication; in the current form the central claim rests on a single unverified analysis pipeline."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should read this one. The group has done a careful STM study of bulk FeSe with superconducting Nb tips and found gap modulations at the expected intra-pocket wavevectors, along with YSR states around subsurface magnetic defects and a faint Josephson-current modulation. The attribution of these oscillations to pair-breaking scattering interference rather than a PDW is new for a density-wave-free superconductor, and the phase-referenced QPI diagnostic they introduce is a genuinely useful tool. The raw coherence-peak traces in Fig. 3(j) show a real 0.12 meV oscillation, which is more direct evidence than most gap maps in this literature. Credit is earned for energy resolution and for choosing a clean system where composite PDW explanations are unlikely.\n\nThe soft spots are real but not disqualifying. The gap maps are outputs of a fitting pipeline with no error bars, and the stress-test worry holds up: a spatially varying Nb tip gap would produce exactly the in-phase Δ+ / Δ− pattern they report, and the symmetrized/antisymmetrized analysis cannot distinguish that from PBSI. The raw traces are the best counter, but they are still after tip-gap deconvolution, and no control measurement or spatial map of the fitted Nb gap is shown. The phase-referenced QPI sign pattern is consistent with PBSI but not unique to it—any particle-hole symmetric modulation yields the same structure, so it is an internal consistency check rather than an independent arbiter. The Josephson modulation is admitted as faint, and there is no quantitative comparison to the PBSI model or any released code/data. All of these are addressable with more analysis and controls.\n\nVerdict: this deserves a serious referee. It is not a finished case, but it is an important experimental result with a plausible alternative interpretation that the authors have not yet excluded. I would bring it to a reading group specifically to discuss the common-mode tip-gap pitfall—it is a nice example of why fitting outputs need systematic-error controls. I'd cite it cautiously in any PDW discussion, and I'd want the revision to include error bars, a control for tip-gap variation, and ideally a direct fit to the PBSI lineshape.","headline":"A serious, well-executed STM study that makes the strongest case yet for PBSI in a clean superconductor, but the central gap maps rest on fit outputs without error bars, and the stress-test concern about a spatially varying tip gap is legitimately unresolved.","tokens_in":10930,"tokens_out":2174,"would_cite":true,"duration_ms":25998,"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":"Bulk FeSe's periodic superconducting-gap modulations arise from magnetic-impurity pair-breaking scattering, not from pair density wave order.","keywords":["pair-breaking scattering interference","FeSe superconductor","superconducting gap modulation","scanning tunneling microscopy","Yu-Shiba-Rusinov states","phase-referenced quasiparticle interference","pair density wave","Josephson tunneling"],"falsifier":"Re-measure the same field of view with a normal metallic tip: if coherence-peak positions extracted directly from the raw spectra (without deconvolution and two-gap fitting) do not reproduce the 1.3 nm^-1 particle-hole-symmetric oscillation, or if the oscillation vanishes in regions far from subsurface magnetic impurities, the PBSI attribution fails.","tokens_in":9983,"feed_emoji":"🧲","tokens_out":3828,"duration_ms":41019,"temperature":0.7,"pith_summary":"This paper tries to establish that the periodic, particle-hole-symmetric modulations of the superconducting gap observed in bulk FeSe are caused by pair-breaking scattering interference (PBSI) around magnetic impurities, not by pair density wave order. Using scanning tunneling microscopy with superconducting Nb tips for enhanced energy resolution, the authors resolve two distinct superconducting gaps and find both modulate in space with wavevectors matching intra-pocket scattering vectors of the Fermi surface. Around subsurface magnetic scatterers identified by Yu-Shiba-Rusinov states, the gap modulations have an amplitude of at least 0.12 meV and are accompanied by spatial modulation of the Josephson current. A phase-referenced quasiparticle interference analysis independently reproduces the expected PBSI sign pattern. If correct, the work establishes PBSI as a viable alternative explanation for gap modulations in superconductors that lack preexisting charge or spin density waves, urging caution before attributing such modulations to finite-momentum pairing.","feed_headline":"FeSe gap ripples come from magnetic scatterers, not pair density waves","feed_subtitle":"Superconducting-tip STM shows 0.12 meV gap ripples around magnetic impurities — no pair density wave needed.","key_machinery":"The pair-breaking scattering interference (PBSI) formalism, which describes how Bogoliubov quasiparticles scattering between banana-shaped constant-energy contours connected by a wavevector Q modulate the local density of states and hence the coherence-peak positions: for sign-preserving scattering (s=1) in the presence of magnetic scatterers, the two coherence peaks at ±Δ± shift in a particle-hole symmetric manner with amplitude set by the scattering potential V. The key experimental tool is phase-referenced quasiparticle interference, g̃(q,E1,E2)=|g(q,E1)|cos[θ(q,E1)−θ(q,E2)], which reveals whether modulations at two energies are in or out of phase; combined with superconductive Nb tips th","core_discovery":"The paper reports that in bulk FeSe, subsurface magnetic impurities produce spatially dispersing Yu-Shiba-Rusinov states, and around these impurities the superconducting gaps on both the alpha and epsilon Fermi-surface pockets modulate periodically with wavevectors equal to the intra-pocket scattering vectors Q0 and Q1. The modulations are particle-hole symmetric and in-phase between positive and negative energies, with amplitude at least 0.12 meV, and the maximum Josephson current modulates with the same periodicity. Phase-referenced quasiparticle interference images show a hot spot at Q0 whose sign pattern across four energies matches the s=1 PBSI prediction. The authors conclude these obs","pith_inferences":["If PBSI is generic, then some previously reported 'pair density wave' gap modulations in superconductors without coexisting density waves may actually be PBSI; re-examining those materials with phase-referenced QPI could settle the interpretation.","Because the lower bound on the scattering potential is only 0.12 meV, even very weak magnetic impurities can generate detectable gap modulations, implying PBSI-induced ripples may be ubiquitous in any superconductor with magnetic disorder and would show up at Fermi-surface nesting wavevectors.","A testable extension: deliberately introducing magnetic impurities (e.g., by electron irradiation or controlled doping) should increase the density of PBSI modulations, while mapping around purely non-magnetic impurities should show no s=1 PBSI gap modulation.","The same phase-referenced QPI technique could be applied to extract the sign of the order parameter between different pockets in multiband superconductors, potentially distinguishing s± from s++ pairing even in the presence of strong disorder."],"forward_implications":["Gap modulations of order 0.1 meV or larger can arise without any finite-momentum pairing, so claims of pair density wave order in superconductors lacking preexisting density waves must rule out PBSI.","The modulation amplitude is set by the impurity scattering potential and is independent of Zeeman energy, giving a quantitative prediction that can be checked in other materials.","Phase-referenced quasiparticle interference provides a practical, independent method to identify PBSI-induced gap modulations, applicable even when direct gap mapping is difficult.","In FeSe specifically, the intra-alpha-pocket scattering vector Q0 is the dominant channel, while the epsilon-pocket signal is weaker but consistent, constraining models of orbital-selective pairing.","The observed Josephson current modulation around magnetic impurities confirms that PBSI affects superfluid density as well as the single-particle gap, linking the two measurable signatures."],"fun_headline_variants":["FeSe ripples: magnetic scatterers, not pair density waves","Josephson STM reveals FeSe gap ripples from magnetic YSR states","Gap modulations in FeSe now tied to pair-breaking scattering interference","Magnetic impurities shape FeSe gaps - no PDW needed"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the deconvolution of the Nb-tip spectra and the two-gap fitting procedure faithfully recover the true local density of states of FeSe; if that pipeline itself produces a spurious 1.3 nm^-1 oscillation, the central gap-map evidence collapses.","fun_headline_variants_meta":{"raw":{"variants":["FeSe ripples: magnetic scatterers, not pair density waves","Josephson STM reveals FeSe gap ripples from magnetic YSR states","Gap modulations in FeSe now tied to pair-breaking scattering interference","Magnetic impurities shape FeSe gaps - no PDW needed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000586,"raw_usage":{"total_tokens":2569,"prompt_tokens":701,"completion_tokens":1868,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":445,"completion_tokens_details":{"reasoning_tokens":1803}},"tokens_in":445,"tokens_out":1868,"duration_ms":15525,"temperature":1.0,"reasoning_tokens":1803,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T15:36:14.911391+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-measure the same field of view with a normal metallic tip: if coherence-peak positions extracted directly from the raw spectra (without deconvolution and two-gap fitting) do not reproduce the 1.3 nm^-1 particle-hole-symmetric oscillation, or if the oscillation vanishes in regions far from subsurface magnetic impurities, the PBSI attribution fails.","supporting_citations":[],"review_version":1}