{"id":"d6fb418c-aca9-4d14-8b3e-19e0edf26240","arxiv_id":"2510.12667","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"DMRG and exact-diagonalization results on up to 72 sites indicate a gapped zero-field ground state for the easy-axis triangular-lattice Heisenberg model with α≲0.3–0.5, and a crossover/transition to gapless behavior at larger α.","lead":"This paper uses large-scale numerical simulations to argue that a triangular-lattice magnet near the Ising limit has a gapped 'spin solid' ground state at zero field, not a supersolid with spontaneous transverse spin order. It also shows that honeycomb- and square-lattice versions reproduce much of the magnetization curve but stay gapless at the corresponding point.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gapped-solid claim rests on finite-size extrapolations whose asymptotic form is untested; a direct one-magnon-gap extrapolation is needed.","rationale":"The reader's weakest assumption matches mine: the finite-size extrapolation. I find no basis to call the authors fraudulent or the argument circular; the magnetization-curve and stiffness data are plausibly consistent with a small gap. But the claim's load-bearing support is the intercept of a quadratic fit and a 4-point 1/N extrapolation. A direct gap measurement is the natural arbiter, because it avoids both fit-shape and scaling-ansatz assumptions. The discrepancy between α* estimates (0.5 vs 0.3) and the authors' own 'transition or crossover' caveat support keeping the verdict conditional rather than accepting the phase diagram as established. The effective-model comparison (Sec. II) is interesting but does not bear directly on the TL gap, since the mapping is only approximate. I therefore agree with the reader's conditional verdict and suggest no change.","tokens_in":13550,"tokens_out":6211,"duration_ms":53210,"concrete_test":"At α=0.1 on the triangular lattice, directly compute the one-magnon gap Δ1=E(S_z=1)-E(S_z=0) with DMRG for N=72, 96, and 120, and extrapolate Δ1(N)=Δ∞+a/N+b/N^2. Also compute ρ_s at the same sizes for twist angles θ=0.05 and 0.1 to verify finite-θ convergence. If Δ∞≤0 within error, or ρ_s extrapolates clearly positive, then the h*=Δ1>0 claim for α≪1 is not supported. This tests the extrapolation scheme directly instead of relying on fit intercepts.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central new result — a gapped spin solid at h=0 for α≪1 — is inferred from two finite-size extrapolations, neither of which is validated against a direct gap or against larger sizes. Section III A / Appendix A extract h* by quadratic fits to small-m DMRG/ED data for N=30–72. h* is an intercept of a fitted polynomial, not a direct measurement of Δ1=E(S_z=1)-E0. The text notes the data 'collapse quite well', but a good collapse does not fix the intercept; a gapless m∝h response with downward curvature can produce a spurious nonzero h* under a quadratic fit. No error bars are given for h* in Fig. 5. Section III B extrapolates ρ_s/(αJ) with a 1/N ansatz from ED N=18,30,36 plus a single DMRG point at N=48. With only four points and no stated curvature correction or twist-angle convergence (θ=0.1 is fixed), the vanishing extrapolated ρ_s for α<0.3 is fragile. The two probes disagree on the boundary (α*≈0.5 from h* vs ≈0.3 from ρ_s), leaving the phase diagram Fig. 7 internally under-constrained. The authors themselves state it is 'beyond present numerical capabilities to clarify whether we are dealing with a transition or a crossover'. At α=0.1 the conclusion may still be right, but it is not settled by data in this paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the spin-1/2 easy-axis (XXZ) Heisenberg model on the triangular lattice (TL), with exchange anisotropy α = J⊥/Jz < 1 and a longitudinal field h, motivated by experiments on K2Co(SeO3)2 (α ≈ 0.07). It first compares the full TL model with effective models on honeycomb and square lattices obtained by freezing one third of the spins; at the correspondence point (m = 1/2 on the bipartite lattices, which corresponds to h = 0 on the TL), these effective models remain gapless with finite transverse magnetization m⊥. The central and more conjectural claim is that the TL model itself has a gapped ground state at h = 0 for α ≪ 1, with h* = Δ1 = ζαJ, and a transition/crossover to a gapless supersolid at α* ≲ 0.5. Evidence is drawn from (i) the zero-crossing of polynomial fits to magnetization curves m(h) from ED/DMRG on N = 30–72 (Sec. III A, Appendix A), and (ii) the 1/N-extrapolated spin stiffness ρs/(αJ), which vanishes for α ≲ 0.3 (Sec. III B). LSWT is shown to fail at h ∼ 0, with the failure attributed to effective magnon repulsion. The resulting phase diagram (Fig. 7) contains a gapped spin solid for h < h*(α). The authors explicitly concede that a transition vs a crossover at α* cannot be distinguished with present numerics.","tokens_in":13969,"tokens_out":18656,"duration_ms":145434,"significance":"Should the gapped-solid scenario survive, the paper would help settle an active controversy: the easy-axis TL model at α ≪ 1 would not be a supersolid at h = 0, directly affecting the interpretation of KCSO neutron-scattering and thermodynamic experiments, and it identifies a qualitative failure mechanism for LSWT (magnon repulsion / Mott-like gap). The numerical work is substantial: systematic ED/DMRG magnetization curves up to N = 72, a spin-stiffness analysis including a DMRG point at N = 48, and finite-size scaling of m⊥ on the effective honeycomb model. The multi-probe design (m(h), ρs, m⊥) and the explicit α > 0 vs α < 0 comparison are strengths, as is the falsifiable prediction h* = ζαJ with ζ extracted from data. The weaknesses are that both gap diagnostics are extrapolations with untested asymptotic forms and no error bars, and the two probes give inconsistent α* values (≈0.5 vs ≈0.3). The central claim is therefore plausible but not established by the present data; the stress-test concern about the unvalidated finite-size extrapolation lands.","major_comments":[{"comment":"The gapped-solid claim rests on the marginal field h* extracted as the zero-crossing of a polynomial fit to m(h) for N = 30–72 (Fig. 4; Appendix A), not on a direct measurement of the one-magnon gap Δ1 = E(S_z = 1) − E0. For a gapless system with downward curvature in m(h), a polynomial fit can produce a spurious positive zero-crossing; the good data collapse does not fix the intercept. At α = 0.1 the extrapolated h* ≈ ζαJ ≈ 0.01–0.02 J (Fig. 5) is comparable to or smaller than the finite-size gap scale expected for a gapless system on the N = 72 cluster (∝ αJ/L, of order 0.07 J), so h* is not yet separated from finite-size effects. Fig. 5 has no error bars and the fit form/window are not specified. I recommend adding a direct one-magnon-gap analysis: extrapolate Δ1(N) = E(S_z = 1) − E0, already available in the DMRG runs restricted to S_z^tot ≤ 4, and test the sensitivity of h* to the f","section":"Sec. III A / Figs. 4,5; Appendix A"},{"comment":"The spin stiffness is extrapolated to N → ∞ with a bare 1/N ansatz using only four points (N = 18, 30, 36 by ED; N = 48 by DMRG), a fixed twist θ = 0.1 with no convergence check in θ, and no error bars. The conclusion that ρs/(αJ) → 0 for α ≲ 0.3 is fragile; alternative scalings (e.g., 1/N^2, exponential convergence, or an added curvature term) should be tested and shown not to change the zero crossing. In addition, the two probes are quantitatively inconsistent — α* ≈ 0.5 from h*(α) in Sec. III A vs α* ≈ 0.3 from ρs in Sec. III B — which under-determines the boundary drawn in Fig. 7. This discrepancy should be reconciled or presented as an explicit uncertainty range for α*.","section":"Sec. III B / Fig. 6"},{"comment":"The authors state that it is 'beyond present numerical capabilities to clarify whether we are dealing with a transition or a crossover at a particular α*.' This ambiguity is load-bearing: the abstract and the phase diagram assert a gapped GS phase for h < h*(α), but a genuine thermodynamic phase requires h*(∞) > 0. Given the strong size dependence in Figs. 4 and 6 and the mismatch between the two probes, the abstract's wording ('confirm the existence of the gap', 'indicate a transition/crossover') overstates what the data establish. Unless the direct-gap extrapolation requested above is supplied, the central claim should be presented as 'consistent with a gapped solid'.","section":"Sec. V / Fig. 7"}],"minor_comments":[{"comment":"The procedure producing h* is not fully described: specify the polynomial degree, the h/(αJ) fitting window, and whether the fit is applied to the pooled finite-N data or to an extrapolated m(h). Add error estimates (e.g., bootstrap over N or over the fit range).","section":"Fig. 5 (Sec. III A)"},{"comment":"The stiffness formula appears as 'ρs = (1/N)∂²E0/∂²θ', which should read ∂²E0/∂θ² evaluated at θ = 0. Also, the fixed value θ = 0.1 is used without a convergence check; a brief test of θ-dependence would strengthen the DMRG points in Fig. 6.","section":"Eq. (3), Sec. III B"},{"comment":"In HJ,BC the two hopping terms are printed identically (a_i a†_j appears twice); presumably the second should be a†_i a_j. Please correct and verify the subsequent algebra.","section":"Eq. (8), Sec. IV"},{"comment":"The caption says 'DSSP' in two places; the text defines DSSF (dynamical spin structure factor). Unify the acronym.","section":"Fig. 2 caption"},{"comment":"The legend 'N 18, 30, 36, ∞, 48' is confusing: ∞ sits between cluster sizes, and the DMRG N = 48 points are listed last although the text introduces them as crosses. Reorder and explain the symbols in the caption.","section":"Fig. 6 legend"},{"comment":"The abstract says 'several additional numerical studies ... confirm the existence of the gap at α ≪ 1.' The cited literature is divided (Refs. 13–15 and 27 report gapless/supersolid behavior in parts of this regime). 'Confirm' overstates the current state of evidence; consider 'support' or 'are consistent with'.","section":"Abstract / Sec. I"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is in scope for a condensed-matter journal and addresses a timely, genuinely contested question. The authors are honest about the limitations, but the headline claim is pushed further in the abstract and in Fig. 7 than the finite-size evidence allows. The decisive missing piece is a direct one-magnon-gap extrapolation; this is feasible from data already computed (DMRG restricted to S_z^tot ≤ 4). If the authors add it and tighten the abstract, I would support acceptance. If not, the paper should at minimum be reframed as evidence 'consistent with' a gapped solid."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper follows up on the authors' earlier work [25,26] on the easy-axis triangular-lattice Heisenberg model. What's actually new is a larger-DMRG dataset (N=72) for the magnetization curve, a spin-stiffness analysis, and a comparison of the full lattice to effective honeycomb and square-lattice models. The effective-model comparison is the cleanest part: HcL and SqL remain gapless at the correspondence point, while the TL data point to a gapped regime at small α. That contrast is informative and helps delimit where the gap is special to the triangular geometry.\n\nThe paper is honest. It explicitly states it cannot distinguish a transition from a crossover, and it gives the raw data for the magnetization curves in an appendix. The internal consistency of the two probes (m(h) and ρs) is a point in its favor.\n\nThe soft spot is the extrapolation. The gap h* is not measured directly; it's the intercept of a quadratic fit to m(h) data up to N=72 with only S_z^tot ≤ 4. A gapless m ∝ h response with downward curvature could produce a spurious nonzero intercept under that fit. The spin stiffness ρs is extrapolated with a 1/N ansatz from ED N=18,30,36 and one DMRG point at N=48, and that extrapolation gives α*≈0.3 while the h* probe gives α*≈0.5. The discrepancy is not discussed in detail. Without error bars or a direct gap calculation (e.g., from the energy difference E(S_z=1)-E0), the central claim remains plausible but not settled.\n\nThe paper's own text flags the limitation: 'beyond present numerical capabilities to clarify whether we are dealing with a transition or a crossover.' That's candid but also an admission that the phase diagram in Fig. 7 is under-constrained.\n\nWho is this for? Researchers working on KCSO and on the phase diagram of the XXZ model on the triangular lattice. They'll get a useful update with a coherent set of numerical results, even if the headline conclusion is not new.\n\nFor peer review: yes, send it out. The question is contested, the data are substantial, and the effective-model comparison is worth refereeing. The referee should push for direct gap estimates and error-aware extrapolations, but this is not a desk-reject paper.","headline":"A careful but not decisive numerical case for a gapped ground state in the easy-axis triangular lattice; the effective-model comparison is the cleanest part, the extrapolations are the weak link.","tokens_in":14417,"tokens_out":2418,"would_cite":false,"duration_ms":19427,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":null,"created_at":"2026-08-04T09:52:40.482552+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":null,"supporting_citations":[],"review_version":1}