{"id":"7d3d4244-4ae0-4d6c-9ca7-a1a3dc321ea1","arxiv_id":"2505.23502","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"KCeSe2 is characterized as a triangular-lattice magnet with dominant ferromagnetic Kitaev interaction and field-induced phase transitions at 1.67 T and 3.8 T.","lead":"This paper studies a new magnetic material, KCeSe2, whose atoms form a triangular lattice and behave like tiny quantum magnets. The authors find that the magnetic interactions are dominated by a bond-dependent Kitaev coupling, and that magnetic fields cause two sharp phase transitions, making the material a promising testbed for exotic quantum spin physics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dominant-K claim rests on a two-parameter INS fit to two spectral features with no uncertainty or uniqueness check; alternate (J±±, Jz±) families may change K drastically.","rationale":"The reader's weakest assumption focused on the completeness of the four-parameter nearest-neighbor Hamiltonian and the representativeness of finite-size clusters. My concern is more specific and more directly tied to the Kitaev-dominance claim: the off-diagonal parameters that determine K are extracted from an underconstrained fit to two spectral peaks, with no uncertainty quantification. This is a load-bearing soft spot because Eq. 2 is a linear map: any spread in J±± or Jz± translates directly into spread in K. If alternate parameter sets fit the INS equally well but give K ≈ 0 or K > 0, the paper's headline conclusion would not survive. The DMRG consistency with the field-induced transitions does not rescue the claim because the same Hamiltonian is used for both fitting and prediction; it is a useful internal check but not an independent test of the model. I therefore partially agree with the reader: the identified weakness is real, but the most specific pressure point is the INS parameter indeterminacy rather than the general model-completeness concern. The reader's CONDITIONAL verdict remains appropriate: the paper's central claim should be accepted only after a demonstrated uniqueness/robustness analysis of the extracted Kitaev coupling, or with a clearly stated caveat that K is model-dependent. I do not see grounds to move to REJECT or UNVERDICTED because the experimental data and the broad physics (stripe-yz order, strong anisotropy, field-induced transitions) are credible and the model is at least a viable working description. A focused sensitivity check could settle the matter without redoing the experiments.","tokens_in":14503,"tokens_out":4525,"duration_ms":44718,"concrete_test":"Perform a systematic parameter sensitivity analysis for the INS fit: compute ED (or LSW) powder-averaged S(Q,ω) for a dense grid of (J±±, Jz±) covering the full phase-diagram range (e.g., −1 to 1 K) while keeping J± = 0.40 K and Jzz = 0.61 K fixed, and evaluate a weighted least-squares residual against the measured 0.3–0.4 meV intensity window and the full |Q| dependence. Then bootstrap the observed intensities from counting statistics and propagate to a posterior for K = −2J±± − √2 Jz±. If the 95% credible interval for K crosses zero or includes K > 0, the 'dominant ferromagnetic Kitaev' claim fails; if it excludes zero with |K| > J, the claim is robust. A secondary check: refit the same datasets with the symmetric Γ-only model (set J±± = Jz± = 0) and report the increase in χ²; if the increase is within noise, the Kitaev terms are not required by the data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numbers K = −1.82 K and J = 1.34 K are not directly measured; they are linear combinations (Eq. 2) of four fitted exchange parameters. The diagonal part (J± = 0.40 K, Jzz = 0.61 K) is fitted to bulk magnetization, but the crucial off-diagonal terms are fixed by matching a single powder INS difference spectrum in the narrow range 0.3–0.4 meV (Fig. 3b). Only two features—a weak peak near 0.28 meV and a stronger peak near 0.35 meV—are used to determine two parameters (J±± = 0.58 K, Jz± = 0.47 K) on a 4×6 ED cluster. With two data features and two parameters, the fit has effectively zero overconstraint: no confidence intervals, no χ² landscape, no powder-averaged Q-dependence beyond a single |Q| window. Equation 2 gives K = −2J±± − √2 Jz±, so the sign and magnitude of K are algebraically controlled by exactly these two parameters. A correlated shift, e.g. (J±±, Jz±) = (0.2, 0.8), yields K ≈ −1.53 K; (0.9, 0.1) yields K ≈ −1.94 K; a family with J±± negative could even give positive K if the INS peaks can be reproduced by compensating changes in Jz± and J±. No uniqueness scan or bootstrap is shown. The subsequent DMRG 'predictions' of the 1.67 T and 3.8 T transitions use the same fitted Hamiltonian, so agreement is a self-consistency check, not an independent validation. Therefore the central claim of Kitaev dominance is not established beyond the specific and possibly non-unique parameter point.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports a multi-probe experimental and numerical study of the triangular-lattice antiferromagnet KCeSe2. Magnetization, ESR, powder neutron diffraction, powder INS, specific heat, and AC susceptibility data are combined with full-diagonalization, thermal-pure-quantum, and DMRG calculations. The authors fit the four nearest-neighbor anisotropic exchanges of Eq. (1) to magnetization and INS data, obtaining J± = 0.40 K, Jzz = 0.61 K, J±± = 0.58 K, and Jz± = 0.47 K. Using the linear transformation in Eq. (2), they convert these into J = 1.34 K, K = −1.82 K, Γ = −0.62 K, and Γ′ = 0.21 K, which they interpret as dominant ferromagnetic Kitaev and antiferromagnetic Heisenberg interactions stabilizing stripe-yz order. They further report field-induced transitions at 1.67 T and ~3.8 T along the a-axis, reproduced by DMRG, and interpret the lower-field transition as a quantum critical point with enhanced quantum fluctuations.","tokens_in":14890,"tokens_out":6795,"duration_ms":64257,"significance":"If the fitted Hamiltonian is correct, KCeSe2 would be a rare experimental realization of a triangular-lattice Kitaev system, and the observation of a field-tuned stripe-yz to canted transition at 1.67 T would be an interesting addition to the Kitaev-materials family. The paper's strengths are the combination of several bulk probes on a newly synthesized material and the explicit cross-checks between ED/TPQ/DMRG calculations and thermodynamic data. The DMRG reproduction of the 1.67 T transition is supportive because this field value was not used in the parameter fit. However, as detailed below, the central quantitative claim depends on a fit with limited overconstraint and no uncertainty analysis; the significance of the result therefore hinges on whether the parameter set is unique.","major_comments":[{"comment":"The off-diagonal couplings J±± and Jz± are determined by matching two features (a weak peak near 0.28 meV and a stronger peak near 0.35 meV) of a powder-averaged INS difference spectrum using a 4×6 ED cluster. This is effectively a two-parameter fit to two spectral features; no chi-square landscape, confidence intervals, or alternative-parameter scans are shown. Because Eq. (2) gives K = −2J±± − √2 Jz± and Γ = (−2J± + Jzz − 4J±± + √2 Jz±)/3, the dominant-K and Γ conclusions are algebraically controlled by exactly these two parameters. Correlated shifts of J±± and Jz± can change K substantially and even its sign if the INS peaks can be compensated by other changes in the Hamiltonian. The authors should provide a uniqueness scan over (J±±, Jz±), report the goodness of fit, and give propagated uncertainties for K, J, Γ, and Γ′. Without this, the abstract-level claim of dominant ferromagnetic Kitaev interaction is not quantitatively established.","section":"Low-energy Spin Excitations, Fig. 3(b), Eq. (2)"},{"comment":"The DMRG magnetization curves that reproduce the 1.67 T and 3.76 T features are computed from the same four-parameter Hamiltonian that was fitted to the magnetization and INS data. Agreement is therefore a self-consistency check rather than an independent prediction, although the 1.67 T value itself was not used in the fit. The authors should either describe the DMRG comparison as a consistency check or demonstrate robustness of the transition fields across the fitted-parameter uncertainty region; this distinction matters for the claim that DMRG predicted the transition.","section":"Magnetic field induced phase transition, Fig. 4(a)"},{"comment":"The fit assumes that the four nearest-neighbor exchanges in Eq. (1) form a complete low-energy model. Further-neighbor exchanges, multi-ion anisotropies, or small Dzyaloshinskii-Moriya terms are not considered, and the powder-averaged INS data provide only limited Q-space constraints. This is a correctness risk for the derived K and J values, which are linear combinations of the fitted parameters. The authors should explicitly state this limitation and, if possible, bound further-neighbor couplings, for example from the magnetic ordering wavevector or from a single-crystal dispersion once available.","section":"Model Hamiltonian and Low-energy Spin Excitations"}],"minor_comments":[{"comment":"The phrase “consistent with density matrix renormalization group (DMRG) calculations predictions” is ungrammatical and should read “DMRG calculation predictions” or “predictions of DMRG calculations.”","section":"Abstract"},{"comment":"The symbol i is used for both the site index and the imaginary unit in the phase factors γij and ladder operators S±; this is confusing and should be disambiguated.","section":"Eq. (1)"},{"comment":"The statement that this is “the first observation of field-induced phase transitions in a stripe-yz ordered triangular magnet” should be qualified in light of the field-induced quantum criticality reported in CsCeSe2 (Refs. 41 and 42).","section":"Conclusions"},{"comment":"The order-by-disorder interpretation in the text and abstract is plausible but is supported only by a citation and the model phase diagram; either provide a direct calculation showing the selection of stripe-yz order for these parameters or soften the wording.","section":"Magnetic anisotropy and stripe-yz order"},{"comment":"The INS difference spectrum is shown without error bars; please add them or state the statistical uncertainty so the significance of the 0.28 meV and 0.35 meV peaks can be assessed.","section":"Fig. 3(b)"},{"comment":"The legend for the dM/dH curves is not fully specified; please label the light-blue and dark-blue lines explicitly as experiment and DMRG calculation, respectively.","section":"Fig. 4(a)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of cond-mat.str-el and will be of interest to the Kitaev-materials community. My main concern is parameter uniqueness: the central numbers K and J come from a fit with limited overconstraint and no uncertainty quantification. I recommend requiring a robustness analysis before publication, and I also suggest the authors compare their fitted parameters with the published KCeS2 and CsCeSe2 values, since the bond-dependent anisotropy is expected to be similar across this family."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nYou should know two things: this is a genuinely new material, KCeSe2, with a careful multi-probe experimental study; and the headline number—dominant Kitaev K = −1.82 K—is softer than the abstract suggests.\n\nWhat the paper does well: it determines the stripe-yz ground state with neutron diffraction (with a clean refinement), documents strong easy-plane anisotropy via ESR and magnetization, and maps out field-induced transitions via AC susceptibility, with a candidate QCP near 1.67 T. The specific heat and magnetization calculations are consistent with the fitted model, and the DMRG phase diagram places the compound in the stripe-yz region.\n\nThe soft spot is the parameter extraction. J± and Jzz come from magnetization, but the off-diagonal J±± and Jz± are fixed by matching two features in a powder INS difference spectrum on a 4×6 ED cluster. Two features, two parameters, no error bars, no uniqueness scan. Since K = −2J±± − √2 Jz±, a correlated shift of those two parameters can change K by a few tenths of a kelvin, possibly even its sign. So the central claim of Kitaev dominance is not established at the level the text claims.\n\nThe field-induced transitions themselves are real; the DMRG agreement at 1.67 T is a self-consistency check, not an independent prediction. And the claim of being the first to see field-induced transitions in a stripe-yz triangular magnet conflicts with the cited CsCeSe2 work (Refs 41,42) — an oversight.\n\nThis is a solid experimental paper with an overinterpreted central claim. A serious referee should ask for a full parameter-uncertainty analysis, including simultaneous fits and a demonstration that the INS data rule out alternative (J±±, Jz±) points. Until that's done, the Kitaev numbers should be treated as provisional.\n\nWho gets value? Anyone tracking Kitaev candidates or rare-earth triangular magnets. It deserves a serious referee, but with the expectation of heavy revision. My recommendation: engage, but do not take the Kitaev numbers at face value.","headline":"A new triangular Kitaev candidate with careful bulk characterization, but the dominant-K claim rests on an underdetermined INS fit and should be treated as provisional.","tokens_in":15463,"tokens_out":3304,"would_cite":true,"duration_ms":30738,"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":"KCeSe2 is a triangular-lattice Kitaev magnet with dominant ferromagnetic Kitaev exchange K = -1.82 K, whose stripe-yz order and two field-induced transitions at 1.67 T and 3.8 T are reproduced by DMRG.","keywords":["triangular-lattice magnet","Kitaev interaction","Kitaev-Heisenberg model","stripe-yz order","order by disorder","field-induced phase transition","cerium chalcogenide","DMRG"],"falsifier":"A single-crystal inelastic neutron scattering measurement that resolves the full spin-wave dispersion would falsify the central claim if its mode energies and intensities cannot be reproduced by Eq. 1 with the fitted parameters, or if an equally good fit requires adding a second-neighbor exchange that changes $K$ by more than its statistical uncertainty.","tokens_in":14304,"feed_emoji":"🧲","tokens_out":8166,"duration_ms":85769,"temperature":0.7,"pith_summary":"This paper argues that the triangular-lattice antiferromagnet KCeSe2 is a rare experimental platform for Kitaev physics beyond honeycomb magnets. From magnetization, neutron diffraction, inelastic neutron scattering, and specific heat, the authors fit a four-parameter anisotropic spin Hamiltonian and convert it to an extended Kitaev-Heisenberg form, obtaining a dominant ferromagnetic Kitaev interaction $K = -1.82$ K alongside antiferromagnetic Heisenberg $J = 1.34$ K. They show that this combination stabilizes stripe-yz magnetic order through order-by-disorder, even though both interactions individually promote frustration. Applying a magnetic field along a Kitaev bond, they find phase transitions at 1.67 T and 3.8 T that density-matrix renormalization-group calculations reproduce as stripe-yz to stripe-canted to polarized. If the model is right, KCeSe2 offers a way to study Kitaev interactions on a triangular lattice and field-driven quantum critical fluctuations near 1.67 T.","feed_headline":"Kitaev exchange dominates in triangular-lattice KCeSe2","feed_subtitle":"Data pin ferromagnetic Kitaev K = -1.82 K; fields then drive transitions at 1.67 T and 3.8 T.","key_machinery":"The load-bearing object is the four-parameter nearest-neighbor anisotropic spin Hamiltonian of Eq. 1, with couplings $J_{zz}$, $J_{\\pm}$, $J_{\\pm\\pm}$, $J_{z\\pm}$, plus the exact transformation of Eq. 2 to the extended Kitaev-Heisenberg parameters $J$, $K$, $\\Gamma$, $\\Gamma'$. Since $K = -2J_{\\pm\\pm} - \\sqrt{2}J_{z\\pm}$, the fitted positive values $J_{\\pm\\pm}=0.58$ K and $J_{z\\pm}=0.47$ K are what turn the empirical fit into a dominant ferromagnetic Kitaev interaction. The same parameters feed exact-diagonalization and thermal-pure-quantum specific-heat checks, linear spin-wave comparisons, and a $6\\times 30$ DMRG cylinder whose entanglement entropy and static spin structure factors mark the stripe-yz phase and the field-induced transitions. Order-by-disorder -- quantum fluctuations lifting the classical degeneracy among stripe states -- is invoked to explain why frustration plus Kitaev exchange still selects stripe-yz order.","core_discovery":"The core discovery is that KCeSe2 is a triangular-lattice material where a ferromagnetic Kitaev term is the largest exchange scale. Using the anisotropic spin model of Eq. 1, the authors determine four exchange couplings ($J_{\\pm}=0.40$ K, $J_{zz}=0.61$ K, $J_{\\pm\\pm}=0.58$ K, $J_{z\\pm}=0.47$ K) by matching ESR g-factors, magnetization, inelastic neutron scattering, and specific heat, then map them via Eq. 2 onto an extended Kitaev-Heisenberg model to get $J=1.34$ K, $K=-1.82$ K, $\\Gamma=-0.62$ K, and $\\Gamma'=0.21$ K. This parameter set places KCeSe2 in the stripe-yz ordered phase in both their DMRG phase diagram and the theoretical $J$-$K$-$\\Gamma$ phase diagram, consistent with the neutron-diffraction ground state and the observed $T^3$ specific heat. A field along the $a$-axis produces transitions at 1.67 T and 3.8 T; DMRG magnetization curves and static spin structure factors identify the lower one as a quantum critical point with enhanced fluctuations and the upper one as entry into the spin-polarized state.","pith_inferences":["If the order-by-disorder picture is right, then isostructural Ce chalcogenides with slightly different Ce-Se bond angles should sit at different points on the $K/J$ axis, and a compound with a weaker Kitaev term could fall into a spin-liquid region; a systematic search across that family would be a direct test.","The 1.67 T critical field should leave a signature in field-dependent inelastic neutron scattering: a continuum of excitations appearing near the critical field would support the spinon scenario the paper floats, whereas sharp magnon modes persisting across the transition would argue against it.","The four-parameter model was fitted to powder data, so a single-crystal inelastic neutron scattering experiment resolving the full dispersion would either confirm the fitted couplings or reveal additional exchange paths; such a measurement would also give an independent check on the quoted $K$ value."],"forward_implications":["KCeSe2 joins KCeS2 and CsCeSe2 as Ce-based triangular-lattice magnets where stripe-yz order coexists with strong Kitaev exchange, so the extracted ratio $K/J \\approx -1.36$ makes this a test case for bond-dependent anisotropy on a frustrated lattice.","The order-by-disorder mechanism implies that in zero field the stripe-yz pattern is selected from a degenerate manifold by quantum fluctuations, so tuning the Kitaev coupling, for instance by chemical substitution, could in principle restore a spin-disordered or spin-liquid state.","A magnetic field along the Kitaev bond direction drives two transitions at 1.67 T and 3.8 T: the lower one is identified as a quantum critical point with multi-fluctuation behavior, and the upper one as the crossover into a spin-polarized phase.","Near 1.67 T the DMRG calculations show disordered moment arrangements and strong quantum fluctuations, conditions the paper suggests could host deconfined spinon excitations, though it notes that dynamical measurements are still needed to confirm this.","The extracted $J$, $K$, and $\\Gamma$ place KCeSe2 inside the stripe-yz region of the triangular-lattice $J$-$K$-$\\Gamma$ phase diagram, directly linking the empirical parameters to the theoretical global phase diagram."],"supporting_citations":[{"why":"Supplies the anisotropic spin-1/2 Hamiltonian of Eq. 1 used to model the low-energy magnetism of KCeSe2.","marker":"[47]"},{"why":"Provides the transformation of Eq. 2 between the four exchange parameters and the extended Kitaev-Heisenberg parameters, and the phase diagram in which stripe-yz order is identified.","marker":"[48]"},{"why":"Gives the triangular-lattice spin-liquid and magnetic-order phase diagram used to interpret the fitted parameters.","marker":"[52]"},{"why":"Supplies the exact-diagonalization and thermal-pure-quantum solver used for specific-heat simulations and low-energy spectral calculations.","marker":"[53]"},{"why":"Provides the updated version of the same quantum lattice solver used for those finite-temperature simulations.","marker":"[54]"},{"why":"Provides the linear spin-wave theory code used as a comparison for the inelastic neutron scattering spectra.","marker":"[55]"},{"why":"Supplies the tensor-network library used for the DMRG phase diagram, magnetization curves, and spin structure factors.","marker":"[56]"},{"why":"Provides the corresponding code release of the tensor-network library used in the DMRG calculations.","marker":"[57]"},{"why":"Provides the prior spin-wave analysis of the sister compound KCeS2 whose inelastic neutron scattering lineshape KCeSe2 resembles, anchoring the excitation assignment.","marker":"[40]"},{"why":"Documents the analogous field-induced quantum criticality in CsCeSe2, providing the comparison case for the 1.67 T transition in KCeSe2.","marker":"[42]"}],"fun_headline_variants":["Triangular lattice KCeSe2: Kitaev exchange dominates","KCeSe2 reveals dominant Kitaev coupling on triangular lattice","Field-induced quantum transitions in Kitaev-triangular KCeSe2","Dominant Kitaev term and two phase transitions in KCeSe2"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes the four nearest-neighbor exchange couplings in Eq. 1, with the fitted values, capture all relevant low-energy magnetism of KCeSe2, so no significant further-neighbor couplings, multi-ion terms, or neglected processes shift the derived Kitaev and Heisenberg values, and that the finite clusters used in exact diagonalization and DMRG represent the thermodynamic limit.","fun_headline_variants_meta":{"raw":{"variants":["Triangular lattice KCeSe2: Kitaev exchange dominates","KCeSe2 reveals dominant Kitaev coupling on triangular lattice","Field-induced quantum transitions in Kitaev-triangular KCeSe2","Dominant Kitaev term and two phase transitions in KCeSe2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000278,"raw_usage":{"total_tokens":1696,"prompt_tokens":1033,"completion_tokens":663,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":649,"completion_tokens_details":{"reasoning_tokens":589}},"tokens_in":649,"tokens_out":663,"duration_ms":6753,"temperature":1.0,"reasoning_tokens":589,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:44:44.915302+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single-crystal inelastic neutron scattering measurement that resolves the full spin-wave dispersion would falsify the central claim if its mode energies and intensities cannot be reproduced by Eq. 1 with the fitted parameters, or if an equally good fit requires adding a second-neighbor exchange that changes $K$ by more than its statistical uncertainty.","supporting_citations":[{"cited_title":"Topography of Spin Liquids on a Triangular Lat- tice,","cited_arxiv_id":null,"evidence_quote":"Gives the triangular-lattice spin-liquid and magnetic-order phase diagram used to interpret the fitted parameters."},{"cited_title":"Quan- tum lattice model solver HΦ,","cited_arxiv_id":null,"evidence_quote":"Supplies the exact-diagonalization and thermal-pure-quantum solver used for specific-heat simulations and low-energy spectral calculations."},{"cited_title":"Update of HΦ: Newly added functions and methods in versions 2 and 3,","cited_arxiv_id":null,"evidence_quote":"Provides the updated version of the same quantum lattice solver used for those finite-temperature simulations."},{"cited_title":"Linear spin wave theory for single-q in- commensurate magnetic structures,","cited_arxiv_id":null,"evidence_quote":"Provides the linear spin-wave theory code used as a comparison for the inelastic neutron scattering spectra."},{"cited_title":"Spin-wave dynamics in the KCeS 2 de- lafossite: A theoretical description of powder inelastic neutron- scattering data,","cited_arxiv_id":null,"evidence_quote":"Provides the prior spin-wave analysis of the sister compound KCeS2 whose inelastic neutron scattering lineshape KCeSe2 resembles, anchoring the excitation assignment."},{"cited_title":"Quantum Spin Dynamics Due to Strong Kitaev Interactions in the Triangular-Lattice An- tiferromagnet CsCeSe2,","cited_arxiv_id":null,"evidence_quote":"Documents the analogous field-induced quantum criticality in CsCeSe2, providing the comparison case for the 1.67 T transition in KCeSe2."}],"review_version":1}