{"id":"6a3da28a-66fa-4385-83de-c69d1d4f2a6e","arxiv_id":"2510.08253","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Spatially resolved IETS maps of monolayer NiI2 show stripe-modulated inelastic signals at about 4 meV that the authors identify as electromagnons, providing atomic-scale visualization of these multiferroic modes.","lead":"Low-temperature scanning tunneling microscopy on monolayer NiI2 reveals in-gap excitations whose spatial maps show a stripe modulation attributed to electromagnons, collective modes that mix magnetic and electric order. The work demonstrates a way to visualize these quasiparticles at the atomic scale, relevant for electrically controlled spintronics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mean-field P^2 term's λ and χ_e are unspecified, so the predicted stripe modulation—the electromagnon fingerprint—can be tuned to match experiment.","rationale":"Reading in good faith, the paper delivers a plausible and technically ambitious experiment: temperature-dependent STM shows multiferroic stripes vanishing at Tc ≈ 10.4 K; IETS reveals in-gap features; spatial maps show a stripe modulation at energies below 5 meV that are consistent with the spin-spiral periodicity. The theory side uses a J1-J3 spin model with single-ion anisotropy, DFT-based phonons, and an HPBdG linear-spin-wave treatment. There is genuine independent support: the measured q-vector and Tc are used to fit J1 and J3, and the phonon spectrum is ab initio. The electromagnon identification, however, depends on the spatial modulation fingerprint. That fingerprint is produced by an effective polarization-energy term whose parameters λ and χ_e are not reported. The reader's weakest assumption pinpoints this correctly. I agree that this is the load-bearing soft spot. Because the term is quartic and treated by a mean-field decoupling that effectively renormalizes the spin-spin couplings based on the ground-state spiral, the resulting spatial modulation is sensitive to the decoupling scheme and to the chosen λ/χ_e. Without independent values, the match between theory (Fig. 3f) and experiment (Fig. 3e) is not strong evidence. This does not invalidate the observation of low-energy inelastic modes with a stripe modulation, but it does prevent the claim of 'unambiguous evidence' for electromagnons. A conditional verdict is appropriate, with the requested additional information (λ, χ_e, error bars, and λ→0 test). Hence no change to the reader's CONDITIONAL verdict.","tokens_in":12895,"tokens_out":4476,"duration_ms":57663,"concrete_test":"Perform a first-principles determination of λ and χ_e for monolayer NiI2: compute the ground-state spin spiral, evaluate the local polarization P_i = λ Σ_j S_i × r_ij × (S_j - S_i) and the electronic dielectric susceptibility χ_e (e.g., from DFPT or finite-field). Fix λ and χ_e to these values in Eqs. 6-8, recompute the HPBdG local spectral function, and compare the predicted stripe modulation and its bias dependence to the experimental dI/dV maps in Fig. 3e. In addition, with the same parameters, compute the magnon gap and compare to the THz electromagnon resonance reported in refs. 8 and 13. If the resulting modulation is absent or an order of magnitude weaker than observed, the fingerprint is an artifact of parameter choice rather than robust electromagnon physics.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central identification of electromagnons rests on the spatial stripe modulation in the local spectral function (Fig. 3e,f). The theoretical fingerprint is generated solely by the polarization-fluctuation term H_pol = Σ_i P_i^2/(2ε0χ_e), treated by a mean-field decoupling in Methods Eq. 6. Neither λ (the inverse-DMI coupling) nor χ_e (the dielectric susceptibility) is given a numerical value anywhere in the text or Methods. Since the modulation is a direct consequence of this term, and since the term is quartic and decoupled in a specific way (Eqs. 7-8), the amplitude and even the existence of the modulation in the mean-field spectrum depend on these free parameters. If λ/χ_e is effectively chosen to reproduce the observed stripe contrast, the comparison is circular. The manuscript also states that the second-order term is 'negligible' for the ground state, so it is not constrained by the q-vector or Tc. Moreover, the mean-field decoupling neglects the dynamical polarization fluctuations that are the essence of an electromagnon; a static renormalization of exchange may produce a spatial modulation that is not representative of the true coupled spin-polarization mode. The authors do not show the λ→0 limit, and no independent DFT estimate of λ and χ_e is provided. Thus the distinctive fingerprint—the main new result—is not yet shown to follow from independently fixed physics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports low-temperature STM/IETS measurements on monolayer NiI2 grown on HOPG, combined with DFT, Monte Carlo, and spin-wave calculations. The authors identify low-energy in-gap IETS features as electromagnons and claim that spatially resolved dI/dV maps at these energies reveal a stripe-like modulation that matches theoretical predictions for the electromagnon local spectral function. Temperature-dependent STM imaging is used to locate the multiferroic transition at Tc≈10.4 K and to estimate the exchange parameters J1=-1.25 meV and J3=0.42 meV. The central claim is that this constitutes 'unambiguous evidence of the internal structure of electromagnons.'","tokens_in":13288,"tokens_out":5095,"duration_ms":51318,"significance":"If the identification holds, this would be the first atomic-scale real-space visualization of electromagnons and would establish a potentially general STM-based methodology for probing magnetoelectric collective modes in van der Waals multiferroics. The paper combines a substantial experimental dataset (temperature-dependent imaging, spectroscopic mapping, IETS histograms) with a multi-level theoretical framework (DFT phonons, spin model, Monte Carlo, HPBdG spin waves). The experimental observations themselves are interesting and likely reproducible. However, the strength of the claim depends on the theoretical fingerprint being derived from independently fixed physics, and this is currently not fully established.","major_comments":[{"comment":"The spatial modulation that constitutes the paper's central electromagnon fingerprint is generated by the polarization-fluctuation term H_pol = Σ P_i^2/(2ε0χ_e), treated by the mean-field decoupling of Eq. (6). The coupling λ and the susceptibility χ_e are never assigned numerical values anywhere in the text or Methods. The paper explicitly states that this term has a negligible effect on the ground-state spin configuration, so it is not constrained by the measured q-vector or Tc. Consequently, the amplitude and even the existence of the modulation in the HPBdG spectrum depend on parameters that could in principle be tuned to match the observed stripe contrast. The λ→0 limit is not shown, and no independent DFT estimate of λ or χ_e is provided. The authors should specify these parameters (or their ratio), show how the predicted spectral-function modulation depends on them, and demonstrat","section":"Methods, Eqs. (6)–(8)"},{"comment":"The theoretical maps predict a phonon-related stripe modulation at ~12 mV with half the spin-spiral periodicity, but the experiment shows no visible phonon modulation. The text explains this post hoc by a stronger electron-electromagnon coupling relative to electron-phonon coupling, without specifying these couplings or providing any estimate. Moreover, the I_TDOS used to assign the IETS peaks in Fig. 3(c,d) is computed with equal weights for electromagnon and phonon contributions, an arbitrary choice that the authors acknowledge. Since the peak assignment and the real-space fingerprint are validated by matching to this theory, the discrepancy at 12 mV and the arbitrary weighting need to be addressed quantitatively (e.g., via matrix-element estimates, selection rules, or a parameter scan) rather than explained post hoc.","section":"Fig. 3(e,f) and 'Observation of electromagnon excitations'"},{"comment":"The exchange parameters J1=-1.25 meV and J3=0.42 meV, as well as Tc≈10.4 K, are extracted by fitting data from the same sample: the stripe periodicity gives J3/J1≈-0.3 and the thermal evolution is fitted with 1/σ=A(Tc-T)^{1/4}. No error bars or uncertainties are reported for these quantities. Because these parameters set the energy scale of the computed electromagnon modes that are then compared with the IETS peaks, the lack of uncertainty propagation is a gap in the central argument. Independent determination (e.g., DFT or literature values) or a propagation-of-errors analysis is needed to avoid circularity.","section":"Fig. 2(g,h) and 'Temperature-dependent observation of multiferroicity'"},{"comment":"The mean-field decoupling of the quartic P_i^2 term replaces one spin operator by its ground-state expectation value, resulting in a static renormalization of the exchange and anisotropy tensors. The eigenmodes are thus magnons of an effective static spin Hamiltonian; their interpretation as electromagnons relies on the semiclassical identification p_i = λ/6a Σ m_i × r_ij × (m_j - m_i). The paper does not demonstrate that this captures the coupled spin-polarization dynamics that define an electromagnon, nor that dynamical polarization fluctuations are negligible. A benchmark against a full coupled spin-polarization or spin-lattice calculation, or at least a discussion of the regime of validity of the mean-field decoupling, would strengthen the claim that the observed modulation is a fingerprint of electromagnons rather than of the static mean-field magnons.","section":"Methods, Eq. (6) and following semiclassical interpretation"}],"minor_comments":[{"comment":"The fit 1/σ=A(Tc-T)^{1/4} is described but no fit residuals, goodness-of-fit, or confidence interval for Tc is provided; adding these would help the reader judge the robustness of Tc≈10.4 K.","section":"Fig. 2(g,h)"},{"comment":"The Methods state a peak-to-peak bias modulation of 0.25 mV for 'short-range spectra' and 10 mV for 'long-range spectra'. The distinction is not defined in the main text; please clarify which spectra are shown in Figs. 3(a,b) and how the modulation amplitude affects the IETS peak resolution.","section":"Fig. 3(b) and Methods"},{"comment":"The notation for the polarization term is inconsistent: the main text writes P_i^2/(2ε0χ_e), while Eq. (6) uses λ² multiplied by the spin-spin cross terms. Please define the relationship between λ, χ_e, and the prefactor explicitly.","section":"Methods, Eq. (6)"},{"comment":"The definition of the unitary matrix U_i^(jk) is unclear; 'the two indices in parenthesis (jk) to the spatial indices of the matrix' is difficult to parse. Please write out the transformation explicitly or reference a standard notation.","section":"Methods, HPBdG section"},{"comment":"The phrase 'unambiguous evidence of the internal structure of electromagnons' is stronger than the presented comparison supports, given the theoretical caveats above. A more cautious wording would better match the current evidence.","section":"Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be of interest to the mtrl-sci/cond-mat community, and the experimental dataset is valuable. However, the central claim of 'unambiguous' electromagnon visualization hinges on the P^2-term mean-field parameters λ and χ_e, which are not given, and on a post hoc explanation for the missing phonon modulation. These issues are fixable within the manuscript's scope (by adding parameter values, robustness checks, and an independent benchmark), but they need to be addressed before publication. I would not recommend rejection, as the experimental observations and the theoretical framework are largely sound; the matching between experiment and theory is not yet sufficiently constrained."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: this is the first attempt at real-space STM/IETS visualization of electromagnons in a monolayer multiferroic, and the experimental dataset is substantial. The temperature-dependent stripe order, the in-gap IETS features appearing below Tc, and the spatial maps at −4 mV showing a stripe modulation are all new and well presented. The authors also compute the phonon spectrum in the ferroelectric Cm phase with DFT, which is a nontrivial piece of solid work.\n\nThe soft spot is the theory used to label the stripes as electromagnons. The spatial modulation in the local spectral function comes entirely from the polarization-fluctuation term H_pol = Σ P_i^2/(2ε0χ_e), treated in mean field in Methods Eq. 6. No numerical values are given for λ or χ_e, and the paper states the term barely affects the ground state, so nothing else in the paper constrains it. That means the amplitude and even the existence of the stripe modulation in the computed maps are free choices. The λ→0 limit is not shown. This is a real gap, and it makes the claim of unambiguous fingerprint overdrawn. The mean-field decoupling also replaces the polarization field by its expectation value, so it is arguably not capturing the dynamical hybridized mode one wants to call an electromagnon — it's a static renormalization of exchange.\n\nSecondary issues: J1, J3, and Tc are reported without error bars. J1 and J3 are derived from the same sample's stripe periodicity and Tc and then used to predict the electromagnon energies measured on that same sample, so the comparison is not fully independent. The IETS histogram is compared to the theoretical integrated DOS with equal weights for phonons and electromagnons; the inflection points don't depend on weights, but the experimental peak heights and widths are not modeled. And the predicted phonon spatial modulation at 12 mV is absent, explained afterward by unspecified stronger electron–electromagnon coupling.\n\nThat said, the experimental core looks solid and the paper is a serious attempt. It deserves a thorough referee and probably major revision on the theory side — values for λ and χ_e, a sensitivity analysis, an independent check against the known THz electromagnon energies (refs 8, 13), and some account of why the phonon modulation is invisible. If those are fixed, this could be a strong paper.\n\nRecommendation: send to review. I would bring it to a reading group now for the method discussion.","headline":"First real-space IETS maps of purported electromagnons in monolayer NiI2 — promising, but the theoretical fingerprint depends on unspecified parameters.","tokens_in":13777,"tokens_out":2756,"would_cite":false,"duration_ms":23341,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["75.85.+t","68.37.Ef"],"model":"deepseek-v4-flash","headline":"Inelastic tunneling maps reveal electromagnons in monolayer NiI2.","keywords":["electromagnons","multiferroics","NiI2","scanning tunneling microscopy","inelastic tunneling spectroscopy","spin-spiral order","magnetoelectric coupling","two-dimensional materials"],"falsifier":"Compute the low-energy spin-wave spectrum and spatial dI/dV maps for the same J1-J3 model with the polarization term omitted: if the stripe modulation and the magnon gap persist, they cannot be electromagnon fingerprints. Alternatively, measure the same IETS maps in a NiI2 monolayer on a substrate that strongly screens the electric polarization; a suppressed or shifted stripe modulation would confirm the magnetoelectric origin, while an unchanged stripe would refute it.","tokens_in":12784,"feed_emoji":"🧲","tokens_out":4782,"duration_ms":36436,"temperature":0.7,"pith_summary":"This paper reports the first atomic-scale observation of electromagnons—collective excitations that oscillate in both spin and electric polarization—in a single layer of the multiferroic NiI2. Using low-temperature scanning tunneling microscopy, the authors find sharp inelastic tunneling features below 5 meV that appear only below the multiferroic transition temperature, and spatially resolved maps at those energies show a stripe modulation matching the spin-spiral periodicity. Comparing with spin-model and first-principles calculations, they identify these modes as electromagnons and argue that the stripe-modulated local spectral function is a distinctive real-space fingerprint of the coupled spin-polarization dynamics. If correct, this gives a general method to probe electromagnons at the atomic scale in two-dimensional multiferroics and connects directly to efforts toward electrically tunable, low-energy spintronics.","feed_headline":"Atomic-scale maps expose electromagnons in monolayer NiI2","feed_subtitle":"A stripe-modulated spectral fingerprint identifies these spin-polarization waves, opening a route to electrically tunable spintronics.","key_machinery":"The central object is the electromagnon, a collective mode of the spin-spiral multiferroic in which spin fluctuations are coupled to electric-polarization fluctuations through the inverse Dzyaloshinskii-Moriya interaction. The crucial theoretical step is adding a polarization-fluctuation energy cost, H_pol = sum_i P_i^2/(2 eps0 chi_e), to the J1-J3 Heisenberg spin Hamiltonian and decoupling it in a mean-field (Tyablikov-like) approximation. This term opens a gap in the magnon spectrum and produces a spatial modulation in the local spectral function with half the spin-spiral periodicity—the fingerprint the experiment targets. The calculations use linear spin-wave theory on a rotated Holstein-","core_discovery":"The paper's central claim is that the low-energy inelastic excitations observed in monolayer NiI2 by scanning tunneling microscopy are electromagnons, not ordinary magnons or phonons. The evidence is twofold: spectroscopically, the in-gap features below 5 meV appear and sharpen as the system is cooled through the multiferroic transition, and their energy positions match the inflection points of a calculated density of states that sums electromagnon and phonon contributions; spatially, dI/dV maps at electromagnon energies show a stripe modulation with the periodicity of the spin spiral, which the authors' theory says arises from the coupling between magnetic fluctuations and the emergent elec","pith_inferences":["Editorial inference: the same technique should work in other van der Waals multiferroics and in twisted bilayers, where the spin-spiral periodicity can be tuned; the stripe period in the IETS maps should track the magnetic q-vector if the electromagnon assignment is correct.","Editorial inference: because the modulation is attributed to the polarization-fluctuation term, changing the dielectric environment (e.g., by substrate choice or capping layers) should shift the gap and modulation amplitude—a controlled experiment that would test the mechanism independently.","Editorial inference: if the electromagnon fingerprint is robust, then the strong visibility of the stripe at biases above the electromagnon energy range implies a large electron-electromagnon coupling relative to electron-phonon coupling; quantifying that coupling would be a natural next step."],"forward_implications":["Electromagnons can be detected and spatially resolved with STM-based inelastic tunneling spectroscopy, not only with terahertz optical probes.","The stripe-modulated local spectral function provides a real-space fingerprint for identifying electromagnons in other spin-spiral multiferroics.","The polarization-fluctuation energy term that produces the modulation also stabilizes long-range magnetic order in two dimensions, offering a concrete mechanism to evade the Mermin-Wagner restriction.","Atomic-scale access to electromagnon modes makes it possible to study their interaction with defects, domains, and heterostructure interfaces in a way bulk probes cannot.","The energy scale set by the multiferroic transition (Tc ~ 10.4 K) and the extracted exchange couplings tie the observed modes to the same physics that governs the static multiferroic order."],"fun_headline_variants":["Atomic maps capture electromagnons in single-layer NiI2","Stripe patterns reveal electromagnons in monolayer NiI2","Direct STM evidence for electromagnons in NiI2 monolayer","In-gap modes shown to be electromagnons in monolayer NiI2","Spin-polarization waves spotted at atomic scale in NiI2"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the predicted stripe-modulated spectral function follows from adding a polarization-fluctuation energy cost to the spin Hamiltonian and decoupling it in mean field, without independently fixing the strength of that term (the paper gives no numerical values for lambda or chi_e), so the fingerprint could in principle be an artifact of the model.","fun_headline_variants_meta":{"raw":{"variants":["Atomic maps capture electromagnons in single-layer NiI2","Stripe patterns reveal electromagnons in monolayer NiI2","Direct STM evidence for electromagnons in NiI2 monolayer","In-gap modes shown to be electromagnons in monolayer NiI2","Spin-polarization waves spotted at atomic scale in NiI2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000226,"raw_usage":{"total_tokens":1324,"prompt_tokens":782,"completion_tokens":542,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":452}},"tokens_in":526,"tokens_out":542,"duration_ms":4939,"temperature":1.0,"reasoning_tokens":452,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T10:46:39.983080+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the low-energy spin-wave spectrum and spatial dI/dV maps for the same J1-J3 model with the polarization term omitted: if the stripe modulation and the magnon gap persist, they cannot be electromagnon fingerprints. Alternatively, measure the same IETS maps in a NiI2 monolayer on a substrate that strongly screens the electric polarization; a suppressed or shifted stripe modulation would confirm the magnetoelectric origin, while an unchanged stripe would refute it.","supporting_citations":[],"review_version":1}