{"id":"cf05a20c-71f9-4e70-869f-169e336c318e","arxiv_id":"2412.21003","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Inelastic neutron scattering on VI3 reveals anisotropic spin excitations that the authors fit with a dominant Kitaev interaction, K ≈ -7.8 meV.","lead":"Neutron scattering on the van der Waals ferromagnet VI3 reveals strongly anisotropic magnetic excitations across a wide region of reciprocal space. The authors fit these spectra with a spin model containing a large bond-dependent Kitaev interaction, proposing VI3 as a new platform for Kitaev physics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Kitaev term is not uniquely constrained: a five-parameter LSWT fit with no error bars and no demonstrated exclusion of longer-range or symmetric anisotropic models leaves K as an effective parameter that could absorb omitted physics.","rationale":"The reader's conditional verdict already captures the main risk: the five-parameter fit has no error analysis and the model space is incomplete, so K may be effective rather than microscopic. My stress-test identifies the more precise mechanism, parameter non-identifiability and absence of a quantitative model comparison, which is the load-bearing point. The paper has real strengths: a new high-resolution INS dataset, the explicit symmetry-allowed Hamiltonian, and the successful reproduction of unusual zone-boundary patterns that simpler Heisenberg/XXZ models fail to produce. However, the published text does not provide the quantitative evidence needed to distinguish a genuine Kitaev term from an effective anisotropic exchange, so the reader's CONDITIONAL verdict should stand. I do not see an internal inconsistency that would justify rejection; the concern is about evidential strength, not a demonstrated error.","tokens_in":10086,"tokens_out":5729,"duration_ms":63990,"concrete_test":"Re-analyze the raw or published S(Q,ω) data with a global χ² fit scanning J ∈ [-2,0], K ∈ [-12,0], Γ ∈ [-2,2], Γ′ ∈ [-2,2], A ∈ [-2,0] meV in the same SpinW model, and repeat with Hamiltonian extensions (add J2 or interlayer Jc; drop Γ′; drop K). Report the best-fit parameter covariance and AIC/BIC for each model. If the 95% confidence interval of K excludes |K| < 3 meV and the extended models fit significantly worse, the Kitaev claim is supported; if K is comparable to its uncertainty or an extended Heisenberg/anisotropic model reaches equal χ², the claim is not.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the fitted value K = -7.8 meV in Eq. 1. For this to be a genuine Kitaev interaction, the five parameters J, K, Γ, Γ′, A must be identifiable from the low-energy magnon dispersion and constant-energy slices, and the model space must be complete enough that K is not a stand-in for neglected terms. The paper reports a best fit but gives no χ² map, parameter covariance, or error bars, and the comparative models are only summarized, not quantified. This matters because linear spin-wave dispersions of bond-direction-dependent exchanges are often degenerate in combinations of J, K, Γ, Γ′; the distinctive six-pointed-star and V-shape features may be generated by several different parameter sets. The authors also approximate the distorted lattice as a perfect honeycomb and fit the high-energy orbital branch with parameters they state 'cannot be unambiguously determined,' signaling that the low-energy fit is similarly underconstrained. If, for example, second-neighbor Heisenberg exchange, interlayer coupling, or a symmetric Γ-only model fits the same data within the instrumental resolution, the fitted K would be an effective parameter and the claim of a large Kitaev interaction would not be established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports high-resolution inelastic neutron scattering measurements on the van der Waals ferromagnet VI3, covering a wide range of reciprocal space at 5 K. The measured spin excitation spectra display strongly anisotropic features, including a six-pointed star pattern around the Γ1 zone center and a triangular pattern around the inequivalent Γ1* zone center, as well as a V-shaped dispersion with a bend anomaly. The authors fit these data with a linear spin-wave model described by a nearest-neighbor J-K-Γ-Γ'-A Hamiltonian on a honeycomb lattice, obtaining J = -1.04 meV, K = -7.8 meV, Γ = 0.3 meV, Γ' = -1.2 meV, and A = -0.6 meV, and claim that the low-energy branch is well captured only when a large bond-dependent Kitaev interaction is included. The high-energy branch above ~6 meV is separately simulated with a J-A model for the orbital-quenched state. The paper interprets the results as evidence for a dominant ferromagnetic Kitaev interaction in VI3, discusses the proximity of the fitted parameters to a predicted S=1 Kitaev spin liquid regime (ξ = 0.7289), and connects the findings to the large anomalous thermal Hall effect and the anomalous monolayer Tc enhancement.","tokens_in":10342,"tokens_out":5665,"duration_ms":59729,"significance":"If the central claim is correct, VI3 would be a rare S = 1 van der Waals ferromagnetic Kitaev candidate with a dominant bond-dependent exchange, a result that would substantially broaden the materials platform for Kitaev physics and offer a new framework for understanding the anomalous thermal Hall effect, the pressure/structure sensitivity of magnetism, and the enhanced monolayer Tc. The paper's strengths include the wide momentum coverage of the INS data, the explicit comparison of several symmetry-allowed interaction models, and the identification of texture in constant-energy slices (six-pointed star and triangular patterns) that are qualitatively reproduced by the Kitaev model but not by simpler Heisenberg, XXZ, or DM models. However, the significance is currently tempered by the absence of quantitative fit metrics, the lack of demonstrated parameter identifiability, and the incomplete treatment of model-space alternatives, all of which are needed to establish that the fitted K is a genuine microscopic Kitaev interaction rather than an effective parameter that absorbs neglected physics.","major_comments":[{"comment":"The five parameters J, K, Γ, Γ′, A are fitted to the very same low-energy magnon data that are then used to claim agreement, but the paper provides no goodness-of-fit measure, no error bars, and no parameter covariance or confidence intervals. Without such information, it is impossible to judge whether K = -7.8 meV is uniquely determined by the data or whether other parameter combinations (e.g., different K with compensating Γ and Γ′) reproduce the spectra equally well. This is load-bearing because the central claim of a large Kitaev interaction depends on the identifiability of K from the fit. Please provide a χ² landscape or equivalent analysis showing that the fitted parameter set is a well-isolated minimum and that K is constrained by the data.","section":"Main text, fitting paragraph (Eq. (1) and Fig. 3)"},{"comment":"The Hamiltonian in Eq. (1) includes only nearest-neighbor J-K-Γ-Γ′ and single-ion anisotropy on an ideal honeycomb lattice. The manuscript states that the lattice is approximated as honeycomb because distortions are 'minimal', but no quantitative justification is given, and longer-range exchange (e.g., second-neighbor J2), interlayer coupling, or further-neighbor anisotropic terms are not tested. Since the fitted K could absorb the effect of these omitted terms, the claim of a genuine microscopic Kitaev interaction requires a quantitative demonstration that such terms are negligible or that their inclusion does not change the extracted K. Please provide fits with representative longer-range or interlayer couplings, or a theoretical estimate of their scale.","section":"Model space and lattice approximation (Eq. (1); paragraph 'Considering the minimal lattice distortions...')"},{"comment":"The paper asserts that alternative models (Heisenberg, XXZ, DM) cannot capture the observed anisotropic dispersion, but no quantitative comparison is provided in the main text; the reader is referred to Supplemental Figures without any residual analysis, χ² values, or confidence statements. The central claim is that Kitaev interactions are necessary, not merely sufficient, so the exclusion of simpler models must be demonstrated quantitatively. Please include explicit goodness-of-fit comparisons or at least representative fits of the alternative models with the same data and fitting procedure.","section":"Main text, 'In comparison, alternative models fail...' and Figs. S3–S5"},{"comment":"The parameter ξ is defined via K = sin(2πξ) and J = cos(2πξ) using the already fitted J and K values. Therefore, the statement that VI3 falls within the predicted S=1 Kitaev spin liquid range (0.7150 ≲ ξ ≲ 0.7775) is a restatement of the fitted K/J ratio, not an independent validation from the phase diagram. This is a circularity in the interpretation that should be removed or explicitly reframed as a derived consequence of the fit, rather than presented as corroborating evidence.","section":"Discussion, paragraph on S=1 Kitaev spin liquid phase (ξ = 0.7289)"}],"minor_comments":[{"comment":"The caption states that the high-energy mode parameters (J = -2.9 meV, A = -3.1 meV) 'cannot be unambiguously determined based on the available data.' This limitation is acknowledged only in the caption; it should be stated in the main text as well, and the orbital-state assignment of the high-energy branch should not be presented as being strongly constrained by the fit.","section":"Fig. 3 caption"},{"comment":"There are several typographical artifacts, including 'Kit aev' in the abstract, 'We not e that the bend' in the main text, and a duplicated phrase in the acknowledgements ('supported by the was supported by'). These should be corrected during production.","section":"Abstract and main text"},{"comment":"The text first describes the high-energy mode as 'above 6 meV' and later refers to 'the high-energy mode above 6.5 meV'; please make the energy threshold consistent.","section":"Main text, description of high-energy mode"},{"comment":"The statement 'Considering the minimal lattice distortions observed, we approximate the lattice as a honeycomb structure at low temperatures' would benefit from a quantitative justification, such as the magnitude of the distortion from diffraction data, to allow the reader to assess the approximation.","section":"Main text, lattice approximation paragraph"},{"comment":"The color scales of the experimental and simulated constant-energy slices are not specified; adding a common intensity scale (or stating how intensities are normalized) would make the visual agreement more transparent.","section":"Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript presents an interesting candidate Kitaev material, and the data appear rich. However, the central claim rests on a fit whose identifiability and completeness are not demonstrated; the paper would be substantially strengthened by adding quantitative fit diagnostics, parameter confidence regions, and explicit tests of model-space extensions. I would also encourage the authors to make the fitting code and reduced data available to support reproducibility, given the strong reliance on SpinW simulations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe real news here is the data. The paper reports the first high-resolution inelastic neutron scattering maps of VI3 covering the zone-boundary regions near K, M, and the inequivalent zone center Γ*—previously only the zone-center response was reported. The six-pointed star pattern near Γ1 and the triangular pattern near Γ* are new, striking, and well resolved. That alone makes the paper worth engaging.\n\nWhat the paper does well: the measurements look careful, with constant-energy slices and constant-Q cuts at multiple incident energies, and the qualitative failures of Heisenberg, XXZ, and DM models are plausible from the quoted comparisons. The authors also openly state that the high-energy orbital branch parameters cannot be unambiguously determined, and they note residual intensity at the zone center that linear spin-wave theory does not capture. That honesty is real.\n\nThe soft spots, in proportion. The central claim—dominant ferromagnetic Kitaev interaction, K ≈ −7.8 meV—rests on a five-parameter fit to the same data used to claim agreement. There are no error bars, no covariance estimates, no χ² maps, and no quantified comparison against alternative anisotropic models. The statement that alternative models fail is qualitative, and the supplemental comparisons are summarized rather than shown as numbers. Linear spin-wave dispersions are notoriously degenerate in combinations of J, K, Γ, and Γ′, so K could easily be absorbing second-neighbor exchange, interlayer coupling, or the known low-temperature distortions. The honeycomb approximation is a real simplification, acknowledged in one sentence but not tested. The ξ = 0.7289 argument is also circular: ξ is computed from the fitted J and K, so invoking the Fukui phase diagram is not an independent check. These problems do not kill the paper, but they mean the title's \"signatures\" is the right level of claim, while the abstract's \"large bond-dependent Kitaev interaction\" is a hypothesis rather than an established fact.\n\nI would take this into peer review rather than desk reject. The data are valuable, the contrast with simpler models is meaningful, and the field needs more 3d Kitaev candidates—but a responsible referee must push for error analysis, quantitative model discrimination, and at least one independent test, such as field dependence or continuum scattering. With that revision, the paper could be a useful contribution.\n\nFor your reading group, it is worth discussing both the physics and the methodology, especially the danger of overinterpreting fitted exchange parameters in a frustrated magnet.","headline":"New INS data on VI3 show genuinely anisotropic zone-boundary excitations, but the large fitted Kitaev term is one plausible parametrization rather than a uniquely established interaction.","tokens_in":10926,"tokens_out":1687,"would_cite":true,"duration_ms":20489,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper reports that a large ferromagnetic bond-dependent Kitaev interaction, with $K=-7.8$ meV, dominates the low-energy spin excitations of VI3 and accounts for their anisotropic patterns.","keywords":["Kitaev interaction","inelastic neutron scattering","VI3","van der Waals ferromagnet","honeycomb lattice","spin wave","bond-dependent exchange","S=1 Kitaev"],"falsifier":"Ab initio computation of the exchange constants for the low-temperature distorted structure could show that the microscopic nearest-neighbor Kitaev coupling is much smaller than 7.8 meV, meaning the fitted K is an effective parameter rather than a real bond-dependent exchange. Alternatively, a neutron scan across the 5.5 meV bend at $\\Gamma_1^*$ with better momentum and energy resolution would distinguish a sharp magnon crossing from a broad continuum, the latter indicating physics beyond linear spin-wave theory.","tokens_in":9875,"feed_emoji":"🧲","tokens_out":13098,"duration_ms":121181,"temperature":0.7,"pith_summary":"VI3 is a van der Waals ferromagnet whose V3+ ions form a honeycomb lattice, and it has long resisted explanation by conventional Heisenberg or Dzyaloshinskii-Moriya models. The paper presents high-resolution inelastic neutron scattering maps of its magnetic excitations over a wide region of reciprocal space and shows that the low-energy spin-wave branch is reproduced by a nearest-neighbor Hamiltonian with a large bond-dependent Kitaev term, $K S_i^\\gamma S_j^\\gamma$, in addition to Heisenberg, off-diagonal, and single-ion-anisotropy terms. The fitted parameters are $J=-1.04$ meV, $K=-7.8$ meV, $\\Gamma=0.3$ meV, $\\Gamma'=-1.2$ meV, and $A=-0.6$ meV, so the ferromagnetic Kitaev exchange dominates the low-energy magnetism. If the claim holds, VI3 becomes a rare S=1 van der Waals Kitaev magnet, and its anomalous thermal Hall effect, monolayer ordering-temperature increase, and strong pressure sensitivity would be natural consequences of bond-dependent exchange. This would extend the search for Kitaev physics from 4d/5d $S=1/2$ magnets to a 3d $S=1$ ferromagnet.","feed_headline":"VI3's spin waves signal Kitaev-dominated magnetism","feed_subtitle":"Neutron scattering puts ferromagnetic Kitaev coupling K ≈ -7.8 meV in charge of VI3's magnetism.","key_machinery":"The load-bearing object is the bond-dependent Kitaev interaction $K S_i^\\gamma S_j^\\gamma$ on the honeycomb lattice, where $\\gamma = x, y, z$ labels the three nearest-neighbor bond directions of the edge-sharing iodine octahedra. The paper puts this term inside the symmetry-allowed J-K-$\\Gamma$-$\\Gamma'$-A Hamiltonian and compares linear spin-wave simulations, averaged over the three 120-degree magnetic domains and including the V3+ magnetic form factor, against the measured constant-energy slices and dispersions. It is the Kitaev term that produces the asymmetric, non-ring-like momentum structure; the alternative models fail where the data are most distinctive.","core_discovery":"The paper's central claim is that the low-energy magnon branch of VI3, measured at 5 K across several Brillouin zones, is quantitatively described by the J-K-$\\Gamma$-$\\Gamma'$-A model on a honeycomb lattice, with a ferromagnetic Kitaev coupling nearly eight times larger in magnitude than the Heisenberg coupling. This one model reproduces the six-pointed star pattern near $\\Gamma_1$, the triangular pattern near $\\Gamma_1^*$, the anisotropic V-shaped dispersion with the ~5.5 meV bend along $M$-$\\Gamma_1^*$-$M$, and the evolution of single peaks into two peaks with increasing energy. Simpler models, including Heisenberg, XXZ, and Dzyaloshinskii-Moriya terms, produce ring-like or otherwise incompatible patterns in the same region. The high-energy branch above 6 meV is assigned to the orbital-quenched excited state and simulated separately with J-A parameters that the paper explicitly states cannot be unambiguously determined.","pith_inferences":["One step the paper leaves open is an ab initio computation of the individual exchange paths in the actual low-temperature distorted structure; such a calculation could show whether the fitted $K=-7.8$ meV is a genuine microscopic Kitaev coupling or an effective parameter absorbing other terms.","The separate J-A treatment of the high-energy orbital-quenched branch leaves the connection between the two orbital sectors unexplored; a two-orbital model linking the branches could be tested against the data above 6 meV.","If the Kitaev picture holds, monolayer VI3 becomes a direct place to search for field- or strain-induced Kitaev spin-liquid signatures, a prediction the paper hints at but does not develop in detail."],"forward_implications":["VI3 sits at $\\xi = 0.7289$, inside the range $0.7150 \\lesssim \\xi \\lesssim 0.7775$ proposed for an S=1 Kitaev spin liquid, so moderate tuning by pressure, strain, field, or iodine substitution could push it toward a spin-liquid state.","The dominance of $K$ over $J$ gives a natural explanation for the remarkable sensitivity of VI3's magnetism to pressure and structural distortion: bond-dependent exchange weights change quickly when the lattice distorts.","The large anomalous thermal Hall effect, previously attributed to DM interactions, is more plausibly tied to the Kitaev coupling, since the fitted Kitaev term is much larger than any DM-like term in the model.","The six-pointed star pattern with filled intensity near the zone center resembles the scattering pattern reported for $\\alpha$-RuCl3, so the paper's result suggests that continuum or fractionalized excitations, rather than simple single-magnon scattering, may contribute in VI3 as well."],"supporting_citations":[{"why":"Provides the microscopic edge-shared-octahedra mechanism that motivates bond-dependent Kitaev exchange in honeycomb magnets.","marker":"[16]"},{"why":"XMCD evidence of a substantial orbital moment with different orbital occupations at the two V sites, underpinning the orbital assignment of the low-energy branch.","marker":"[33]"},{"why":"DFT calculation showing the large-orbital-moment state is energetically favorable, used to assign the low- and high-energy magnon branches.","marker":"[35]"},{"why":"Earlier inelastic neutron scattering on VI3 that resolved two dispersive branches near $\\Gamma_1$ and set the baseline that the present measurements extend to zone boundaries.","marker":"[40]"},{"why":"The linear spin-wave simulation code used to compute the J-K-$\\Gamma$-$\\Gamma'$-A spectra with magnetic domains and the V3+ form factor.","marker":"[44]"},{"why":"Reports of the six-pointed star-shaped scattering pattern with filled center in $\\alpha$-RuCl3, the comparison that suggests continuum or fractionalized excitations in VI3.","marker":"[13,14]"},{"why":"Theoretical phase diagram giving the S=1 Kitaev spin liquid parameter range into which the paper places VI3 at $\\xi = 0.7289$.","marker":"[26]"},{"why":"Reported large anomalous thermal Hall effect in VI3, which the paper reinterprets as a possible consequence of the dominant Kitaev interaction.","marker":"[37]"}],"fun_headline_variants":["Kitaev coupling rules VI3's magnetism, neutron data shows","Dominant ferromagnetic Kitaev coupling found in VI3","Neutron scattering shows Kitaev interaction steers VI3's spin waves","VI3 magnon anisotropy points to strong Kitaev interaction","Kitaev physics in VI3: neutron data reveals large coupling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim collapses if the low-energy magnon branch is not accurately described by a single nearest-neighbor J-K-$\\Gamma$-$\\Gamma'$-A model on an ideal honeycomb lattice with the ground-state orbital configuration, because then the fitted K would be an effective parameter absorbing longer-range exchange, interlayer coupling, orbital mixing, and distortion effects rather than a genuine Kitaev interaction.","fun_headline_variants_meta":{"raw":{"variants":["Kitaev coupling rules VI3's magnetism, neutron data shows","Dominant ferromagnetic Kitaev coupling found in VI3","Neutron scattering shows Kitaev interaction steers VI3's spin waves","VI3 magnon anisotropy points to strong Kitaev interaction","Kitaev physics in VI3: neutron data reveals large coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000462,"raw_usage":{"total_tokens":2288,"prompt_tokens":898,"completion_tokens":1390,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":514,"completion_tokens_details":{"reasoning_tokens":1303}},"tokens_in":514,"tokens_out":1390,"duration_ms":10555,"temperature":1.0,"reasoning_tokens":1303,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T23:03:59.527942+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Ab initio computation of the exchange constants for the low-temperature distorted structure could show that the microscopic nearest-neighbor Kitaev coupling is much smaller than 7.8 meV, meaning the fitted K is an effective parameter rather than a real bond-dependent exchange. Alternatively, a neutron scan across the 5.5 meV bend at $\\Gamma_1^*$ with better momentum and energy resolution would distinguish a sharp magnon crossing from a broad continuum, the latter indicating physics beyond linear spin-wave theory.","supporting_citations":[{"cited_title":"Hovančík, J","cited_arxiv_id":null,"evidence_quote":"XMCD evidence of a substantial orbital moment with different orbital occupations at the two V sites, underpinning the orbital assignment of the low-energy branch."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"DFT calculation showing the large-orbital-moment state is energetically favorable, used to assign the low- and high-energy magnon branches."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier inelastic neutron scattering on VI3 that resolved two dispersive branches near $\\Gamma_1$ and set the baseline that the present measurements extend to zone boundaries."},{"cited_title":"Fukui, Y","cited_arxiv_id":null,"evidence_quote":"Theoretical phase diagram giving the S=1 Kitaev spin liquid parameter range into which the paper places VI3 at $\\xi = 0.7289$."}],"review_version":1}