{"id":"95d77bda-b2de-40bb-a86a-6763955c3e41","arxiv_id":"2608.11944","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In few-layer NiI2, increasing thickness rotates and lengthens the spin spiral wavevector, while wrinkles deflect it locally, establishing thickness and strain as tuning knobs for this 2D multiferroic.","lead":"Scientists measured how the magnetic spiral in ultra-thin nickel iodide films changes as the films grow thicker and as they bend. The results show that thickness and local strain are two practical knobs for controlling the material's magnetic and electric properties.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thickness-dependence model is calibrated to the data it claims to explain: interlayer exchange magnitudes are scaled to align with experiment, so Fig. 4d/e does not independently test the dimensional-crossover mechanism.","rationale":"The reader's weakest_assumption identified the same core issue: the analytical thickness model is partially calibrated against the experimental data it is used to explain. The paper's own Supplement II-1 admits both a J3 adjustment and a global scaling of interlayer interactions to align with experiment, so the quantitative agreement in Fig. 4d/e is not a parameter-free test of the dimensional-crossover mechanism. This is load-bearing because the paper's central claim is not merely that q changes with thickness, which the SP-STM data support, but that the change is 'primarily driven by enhanced interlayer exchange energy.' If the fitted interlayer scale absorbs thickness-dependent physics, the computed trajectory could reproduce the data without the proposed mechanism being correct. The curvature-induced deflection claim is less affected because it is supported by a DFT-trained machine-learning potential and shows a qualitative deflection, although its flat-monolayer wavevector differs from experiment; that is a secondary weakness. The reader's CONDITIONAL verdict already reflects this uncertainty, so no verdict change is needed.","tokens_in":22500,"tokens_out":5004,"duration_ms":53209,"concrete_test":"Recompute the q(N) trajectory for N=1-7 without fitting: take interlayer exchange parameters directly from DFT (Table S1, U=4 eV) with no global scaling, keep J3 fixed at 2.63 meV or at the underlying HSE/DFT value, and solve the full N-layer Hamiltonian (or the effective model with no adjustable factor). If the computed wavelength increase and rotation from [110] toward [1-10] disappear, the agreement in Fig. 4d/e is calibration-dependent and the dimensional-crossover mechanism is not independently supported; if the trajectory survives, the scaling concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Supplement II-1 explicitly states that J3 was enhanced from 2.25 to 2.63 meV to match the monolayer wavevector, and that after DFT calculation the interlayer interaction ratios were fixed while their magnitudes were scaled 'to align with experimental data.' This means the comparison in Fig. 4d/e is not an independent prediction: at least one scale factor is calibrated against the same measured q(N) points being explained. The effective (N-1)/N intralayer mapping in Part II-5 also assumes exactly antiparallel interlayer alignment and treats all N layers as identical; any layer-resolved variation (surface versus interior layers, stacking fluctuations, finite-size effects) can be absorbed into the fitted scale factor. The computed q(N) trajectory therefore cannot by itself establish that enhanced interlayer exchange is the driver of the observed crossover. The experimental thickness trend is solid, but the mechanistic attribution requires an uncalibrated test.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a combined SP-STM and multi-scale modeling study of few-layer NiI2 films (1-7 monolayers) grown by MBE. The authors observe that as film thickness increases, the spin-spiral wavelength increases and the wavevector rotates from near [110] toward [1-10], which they attribute to a dimensional crossover driven by enhanced interlayer exchange interactions. They also observe that local wrinkles and bent regions deflect the spin-spiral wavevector, and they use SpinGNN++ spin-lattice dynamics simulations to argue that this deflection arises from curvature-induced modification of intralayer J3/J1 and Kitaev interactions. The experimental work includes FFT-based wavevector extraction from independent terraces, a NbSe2 substrate control, and stacking characterization by STM and STEM. The modeling includes an analytical single-q spin Hamiltonian, DFT extraction of interlayer parameters, and a machine-learned spin-lattice potential for the curvature simulations.","tokens_in":22582,"tokens_out":3497,"duration_ms":41080,"significance":"If the central mechanistic claims hold, the paper establishes two experimentally accessible tuning knobs—film thickness and local curvature—for non-collinear helical magnetism and associated electric polarization in a van der Waals multiferroic. The experimental dataset is a valuable systematic map of q(N) for 1-7 ML NiI2, and the NbSe2 control strengthens the claim that substrate charge transfer does not dominate the thickness evolution. The SpinGNN++ simulations provide an independent route to connect local structural deformation with exchange modification, and the reported test-set MAE of 0.036 meV/atom and reproduced energy rankings support the reliability of that potential. However, the mechanistic attribution of the thickness crossover is weakened by admitted fitting of model parameters to the same experimental data used for comparison, so the paper's central claim requires an uncalibrated or out-of-sample validation to be fully convincing.","major_comments":[{"comment":"The comparison in Fig. 4d/e is not an independent prediction of the dimensional-crossover mechanism. Supplement Part II-1 states that J3 was enhanced from 2.25 meV to 2.63 meV to match the monolayer wavevector and that the interlayer interaction magnitudes were scaled \"to align with experimental data\" after fixing their DFT ratios. Because the same experimental q(N) points are used to set at least one free scale, the calculated thickness trajectory can reproduce the trend by construction. The authors should provide an uncalibrated prediction using the raw DFT interlayer magnitudes (or a band reflecting the Ueff dependence shown in Table S1), or an out-of-sample test such as predicting q for a thickness or stacking variant not used in the fitting.","section":"Supplement Part II-1 and Fig. 4d/e"},{"comment":"The effective intralayer mapping J1,eff = J1 - 2(N-1)/N J2_perp assumes strictly antiparallel interlayer alignment (S_R0 = -S_R0,perp) for every adjacent layer and represents the film as N identical layers with (N-1)/N averaged couplings. The SP-STM step-edge data in Fig. 2f show a half-period offset at one 3ML/4ML step, which is suggestive but does not establish a 180-degree interlayer phase for all interior interfaces of all thicknesses, especially as q rotates with N. If the interlayer phase deviates from pi or if layer-resolved variations (surface versus interior layers, local stacking fluctuations) are present, the fitted scale factor can absorb these effects and the q(N) trajectory in Fig. 4c-e would no longer be a valid test of the interlayer-exchange mechanism. The authors should state what evidence constrains the interlayer phase and quantify the sensitivity of the predicted trajectory to phi and to layer-dependent parameters.","section":"Supplement Part II-5 and Fig. 2f"},{"comment":"The attribution of the curvature-induced deflection to local reduction of J3/J1 and Kitaev interaction K rests on exchange tensors predicted by the SpinGNN++/Spin-Allegro++ potential on a curved lattice. While the potential's overall accuracy is documented for flat monolayer configurations, no direct DFT verification is provided for the bent supercell itself, where bond lengths and angles deviate substantially from the training distribution. A direct DFT calculation of exchange parameters for a few representative locally curved configurations (or at least a comparison of J3/J1 and K predicted by the ML potential against DFT for those configurations) would materially strengthen the claim that curvature-induced exchange modification is the microscopic cause of the observed deflection.","section":"Section 3.3 and Supplement Part III-3"}],"minor_comments":[{"comment":"The caption describes the phase diagram as being in \"J3-K space,\" while the text in Section 3.2 states that the diagram is in \"(J1, K) space\" with J1 set to -1 meV; please make the axes consistent.","section":"Fig. 4b caption"},{"comment":"The main text refers to \"Fig. S6 (online)\" for the bond-length and bond-angle changes in the curved lattice, but in the Supplement Fig. S6 is the NbSe2 substrate control data; the figure numbering for the curved-lattice structural analysis should be corrected.","section":"Section 3.3 and Supplement Part I-5"},{"comment":"The phrase \"turning knobs\" in the final paragraph should read \"tuning knobs\" to match the abstract.","section":"Discussion and Conclusion"},{"comment":"The main text consistently uses \"SpinGNN++,\" while the supplement introduces \"Spin-Allegro++\" as the implementation; please clarify the relationship between these two names in the methods or supplement.","section":"Supplement Part III"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically competent and the experimental core is solid, but the central mechanistic attribution for the thickness crossover is calibrated to the data it claims to explain. I would encourage the editor to request an uncalibrated prediction or an out-of-sample test rather than treating the current Fig. 4d/e comparison as sufficient. The curvature study is more independent and could be strengthened by direct DFT checks on bent configurations. No concerns about data integrity or citation practices beyond the calibration issue."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. First, this is the first systematic SP-STM thickness series, 1 to 7 ML, of the NiI2 spin spiral, and the central observations look solid: the wavevector rotates from near [110] toward [1-10] and the wavelength grows from 1.7 to 2.1 nm as the film thickens, and wrinkles deflect the spiral locally by 40-68 degrees. Second, the thickness-crossover mechanism is backed by a model that is partly calibrated against the data it claims to explain, so Fig. 4d/e should be read as a consistency check rather than an independent prediction.\n\nThe experimental core is the strength. The FFT wavevector extraction is repeated on independent terraces with error bars, the q-2q phase relation confirms spiral order, the half-period offset across a step edge supports antiferromagnetic interlayer coupling, and the NbSe2 control (Supplement I-5) argues convincingly that interfacial charge transfer is not driving the thickness trend. The saturating wavelength and the near-linear rotation of alpha are clean trends; I trust them.\n\nThe wrinkle result is the most novel and the most defensible. The SpinGNN++ potential is trained on DFT data, not on the deflection angles, so the simulated deflection of about 60 degrees is a genuine prediction. The local suppression of J3/J1 and K in the curved region is a concrete mechanistic step. Caveat: the potential's flat-monolayer ground state does not quite match the measured monolayer wavevector, which softens the quantitative claim, but the mechanism stands. The absence of robust topological charge at the deflection is a good control, distinguishing this from the earlier domain-wall textures.\n\nThe soft spot is the thickness model, and the authors deserve credit for admitting it: Supplement II-1 states that J3 was increased from 2.25 to 2.63 meV to match the monolayer, and that interlayer exchange magnitudes were scaled to align with the experimental data. A free scale factor means the computed q(N) trajectory is partly a fit, not a test. The (N-1)/N mapping also assumes strictly antiparallel interlayer alignment and identical layers; the supplement generalizes this to a uniform phase, but the main model does not. The experimental trend stands on its own; the attribution to enhanced interlayer exchange is plausible but not proven by the calculation.\n\nNo data or code are shipped, which matters for a paper this simulation-heavy. Who this is for: anyone working on 2D magnetism, vdW multiferroics, or SP-STM; the paper would make a good reading-group discussion on calibration versus prediction. I would cite the experimental thickness series. It deserves a serious referee: revise and resubmit, with data/code released and the model reframed as corroboration rather than derivation.","headline":"Solid new SP-STM thickness series and wrinkle-deflection observations in NiI2; the crossover mechanism is plausible but partly calibrated to the same data it explains, so treat the theory as corroboration, not proof.","tokens_in":23246,"tokens_out":5928,"would_cite":true,"duration_ms":54627,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["75.70.Ak","75.25.-m","68.37.Ef"],"model":"deepseek-v4-flash","headline":"Thickness and local bending continuously tune the spin-spiral wavevector in few-layer NiI2 films, giving two practical knobs for engineering helical magnetism and the ferroelectric polarization it carries.","keywords":["spin spiral","NiI2","dimensional crossover","spin-polarized STM","multiferroic","interlayer exchange","curvature-induced magnetism","van der Waals magnets"],"falsifier":"A layer-resolved spin-polarized STM experiment that images the spin phase on both sides of a step in a 3- or 4-ML NiI2 terrace would settle the interlayer assumption: finding a phase clearly different from 180 degrees, or finding that the measured 3-ML wavevector falls far from the model's q(N) curve when computed interlayer couplings are used without rescaling, would falsify the thickness-driven dimensional-crossover claim.","tokens_in":22217,"feed_emoji":"🧲","tokens_out":11012,"duration_ms":102649,"temperature":0.7,"pith_summary":"Using spin-polarized scanning tunneling microscopy on films grown one to seven atomic layers thick, this paper tracks how the helical magnetic order of NiI2 changes as the material approaches the two-dimensional limit. It finds that as thickness grows, the spin-spiral wavelength increases from about 1.7 nm toward 2.1 nm and the wavevector rotates continuously from near [110] toward [1-10], a behavior it attributes to antiferromagnetic interlayer exchange rebalancing the frustrated intralayer interactions. The same measurements show that wrinkles and bent regions deflect the spiral locally by 40–68 degrees, an effect reproduced in spin-lattice dynamics simulations and traced to curvature-induced changes in the exchange parameters. If the interpretation is right, film thickness and local strain become two independent tuning methods for non-collinear helical magnetism and the electric polarization that accompanies it in van der Waals multiferroics.","feed_headline":"Spin spiral in NiI2 rotates as film thickens or bends","feed_subtitle":"Spin-polarized STM tracks the magnetic wavevector from 1 to 7 layers, then shows wrinkles bending it by up to 68 degrees.","key_machinery":"The central object is a single-q spin-spiral energy functional for AA-stacked NiI2 built from the frustrated intralayer Heisenberg exchanges J1–J3, a Kitaev term K (a bond-directional anisotropic exchange), and a biquadratic term B, with interlayer Heisenberg exchanges included as effective intralayer modifications for an N-layer film. The load-bearing identity is the mapping of antiferromagnetic interlayer coupling onto effective ferromagnetic intralayer couplings, $J_{1,\\mathrm{eff}} = J_1 - \\frac{2(N-1)}{N} J_2^\\perp$ and $J_{2,\\mathrm{eff}} = J_2 - \\frac{N-1}{N} J_{3,2}^\\perp$, which lets the thickness trajectory be read off a monolayer phase diagram whose incommensurate ICα region supports continuous rotation of the wavevector between crystallographic directions. For curvature, the machinery is a machine-learned spin-lattice dynamics potential that relaxes a bent monolayer and predicts how bending redistributes J3/J1 and K spatially, producing the deflected spiral in simulation.","core_discovery":"The central claim is that the spin-spiral state in few-layer NiI2 is continuously tunable by thickness and by local curvature. In AA-stacked films from 1 to 7 monolayers, the authors observe a monotonic increase of the spiral wavelength and a rotation of its wavevector (the direction and period of the magnetic modulation) from about 7 degrees to about 26 degrees relative to [110], converging toward but not reaching the bulk value, and they attribute this to progressively stronger antiferromagnetic interlayer exchange acting through effective intralayer couplings that scale as (N−1)/N. At wrinkles with curvature around 10–20 inverse micrometers, the spiral direction locally reorients by 40–68 degrees and then recovers; spin-lattice dynamics simulations show that bending changes the J3/J1 ratio and the Kitaev interaction in the curved region. The paper concludes that thickness changes the global balance of competing interactions while local strain perturbs them at the nanoscale, each providing a distinct knob to reorient the helical magnetic order and to manipulate the associated ferroelectric polarization.","pith_inferences":["A testable extension would be to bend single-layer NiI2 by controlled amounts and measure the deflection angle as a function of curvature, checking whether it follows a monotonic relation set by the local change in J3/J1 and K.","The (N−1)/N effective-coupling argument is generic: other frustrated van der Waals helimagnets with antiferromagnetic interlayer stacking should show a similar thickness-driven rotation of the wavevector, and comparing materials such as CoI2 would separate the universal mechanism from NiI2-specific parameters.","Because the deflection is localized at the wrinkle, an array of wrinkles could act as a rewritable template: spin-polarized STM combined with local bias or field might write and erase the local spiral orientation, making the polarization pattern mechanically programmable.","The paper's wavevector extraction assumes a single-q spiral; if a thickness or strain regime with multi-q states or competing domains is found, the FFT-based analysis alone could misread the texture, so phase-resolved real-space imaging of deflected regions would sharpen the claim."],"forward_implications":["At 7 ML the spiral wavelength has nearly saturated while still differing from bulk NiI2, implying the AA stacking of the MBE-grown films defines a distinct thin-film limit whose crossover trajectory depends on stacking order.","Because the 2q charge modulation persists at every thickness, the spin-driven electric polarization is expected to persist in few-layer NiI2, so thickness tunes the magnetic and ferroelectric order parameters together.","Curvature-induced deflection is a local reorientation of a single-q spiral with no topological charge, distinguishing it from previously reported spin-spiral domain walls and offering a strain-based route to pattern multiferroic domains.","Substrate-induced charge transfer is not the controlling factor: NiI2 films on a strongly hole-doping NbSe2 substrate still show nearly the same wavevector as films on HOPG."],"supporting_citations":[{"why":"Supplies the intralayer Heisenberg, Kitaev, and biquadratic spin Hamiltonian and parameter values for NiI2 that the film model extends.","marker":"[37]"},{"why":"Prior analytical treatment of monolayer and few-layer NiI2 spiral phases that the multi-layer effective-coupling model builds on.","marker":"[54]"},{"why":"Establishes the monolayer spin-spiral wavevector, the 2q charge modulation, and the properties of spin-spiral domain walls used as the experimental baseline.","marker":"[34]"},{"why":"Independent spin-resolved imaging of helical order in mono- and bilayer NiI2 used for comparison of wavevector values.","marker":"[36]"},{"why":"Provides the bulk NiI2 magnetic structure and spiral wavelength that define the bulk limit toward which the film crossover moves.","marker":"[26]"},{"why":"Reports single-layer van der Waals multiferroicity in NiI2 and the monolayer-versus-bulk comparison that motivates the thickness study.","marker":"[29]"},{"why":"Provides the machine-learning spin-lattice dynamics framework used to simulate the curved lattice and to extract local exchange interaction maps.","marker":"[55]"},{"why":"Shows monolayer NiI2 on a different substrate with a similar spin spiral, supporting the claim that substrate charge transfer does not dominate the wavevector evolution.","marker":"[35]"}],"fun_headline_variants":["NiI2 spin spiral tuned by thickness and curvature","Thickness and strain bend NiI2 magnetic spiral","Spin spiral in NiI2 shifts with layers and wrinkles","Dimensional crossover reorients NiI2 spin helix","Wrinkles and layer count steer NiI2 spin spiral"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The thickness model assumes adjacent NiI2 layers are strictly antiparallel, treats every layer as identical with averaged (N−1)/N effective couplings, and rescales the computed interlayer exchange ratios once to match the same experimental data being explained; if the interlayer angle, stacking registry, or the scale factor changes with thickness, the predicted q(N) trajectory would not be a clean test of dimensional crossover.","fun_headline_variants_meta":{"raw":{"variants":["NiI2 spin spiral tuned by thickness and curvature","Thickness and strain bend NiI2 magnetic spiral","Spin spiral in NiI2 shifts with layers and wrinkles","Dimensional crossover reorients NiI2 spin helix","Wrinkles and layer count steer NiI2 spin spiral"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000172,"raw_usage":{"total_tokens":1285,"prompt_tokens":966,"completion_tokens":319,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":582,"completion_tokens_details":{"reasoning_tokens":239}},"tokens_in":582,"tokens_out":319,"duration_ms":3213,"temperature":1.0,"reasoning_tokens":239,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:21:35.791109+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A layer-resolved spin-polarized STM experiment that images the spin phase on both sides of a step in a 3- or 4-ML NiI2 terrace would settle the interlayer assumption: finding a phase clearly different from 180 degrees, or finding that the measured 3-ML wavevector falls far from the model's q(N) curve when computed interlayer couplings are used without rescaling, would falsify the thickness-driven dimensional-crossover claim.","supporting_citations":[{"cited_title":"Giant chiral magnetoelectric oscillations in a van der Waals multiferroic","cited_arxiv_id":null,"evidence_quote":"Supplies the intralayer Heisenberg, Kitaev, and biquadratic spin Hamiltonian and parameter values for NiI2 that the film model extends."},{"cited_title":"Molecular beam epitaxy of highly crystalline monolayer molybdenum disulfide on hexagonal boron nitride","cited_arxiv_id":null,"evidence_quote":"Prior analytical treatment of monolayer and few-layer NiI2 spiral phases that the multi-layer effective-coupling model builds on."},{"cited_title":"Possible persistence of multiferroic order down to bilayer limit of van der Waals material NiI2","cited_arxiv_id":null,"evidence_quote":"Establishes the monolayer spin-spiral wavevector, the 2q charge modulation, and the properties of spin-spiral domain walls used as the experimental baseline."},{"cited_title":"Atomic‐scale visualization of multiferroicity in monolayer NiI2","cited_arxiv_id":null,"evidence_quote":"Independent spin-resolved imaging of helical order in mono- and bilayer NiI2 used for comparison of wavevector values."},{"cited_title":"Strain gradient mediated magnetism and polarization in monolayer VSe2","cited_arxiv_id":null,"evidence_quote":"Provides the bulk NiI2 magnetic structure and spiral wavelength that define the bulk limit toward which the film crossover moves."},{"cited_title":"Evidence for a single -layer van der Waals multiferroic","cited_arxiv_id":null,"evidence_quote":"Reports single-layer van der Waals multiferroicity in NiI2 and the monolayer-versus-bulk comparison that motivates the thickness study."},{"cited_title":"Electric field control of moiré skyrmion phases in twisted multiferroic NiI2 bilayers","cited_arxiv_id":null,"evidence_quote":"Provides the machine-learning spin-lattice dynamics framework used to simulate the curved lattice and to extract local exchange interaction maps."},{"cited_title":"Orientation-selective spin-polarized edge states in monolayer NiI 2","cited_arxiv_id":null,"evidence_quote":"Shows monolayer NiI2 on a different substrate with a similar spin spiral, supporting the claim that substrate charge transfer does not dominate the wavevector evolution."}],"review_version":1}