Pith. Sign in

REVIEW

Unraveling square and kagome lattice frustrated magnetism by scalable tensor networks

Strongly Correlated Electrons · Submitted 6 Nov 2026 · 77 references

Pith's verdict: A strong numerical methods paper with an excellent QMC benchmark, but the kagome spin-liquid stability claim rests on a scaling proxy the authors acknowledge but do not independently test.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on November 9, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2611.09881 (version 1)
Pith IDXSNKLEIE, Pith's permanent record of this paper and its review
CitationJutho Haegeman, Xing-yu Zhang, Yuchi He. "Unraveling square and kagome lattice frustrated magnetism by scalable tensor networks." arXiv:2611.09881 (2026). https://pith.science/paper/XSNKLEIE
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Tensor-network fits give 2Δm = 3.95(11) for the kagome Heisenberg model — above the 3 threshold, so the gapless spin liquid is stable — and 2Δm = 1.020(16) for the square lattice, below it.

The reading

This paper claims that the correlation length imposed on infinite tensor-network states by finite bond dimensions can stand in for system size, and that fitting the decay of the squared local magnetization against that length yields the scaling dimension of the spin operator, a quantitative probe of whether a gapless spin liquid is stable. The fitted exponents are $2\Delta m = 3.95(11)$ for the kagome Heisenberg model and $2\Delta m = 1.020(16)$ for the square-lattice $J_1$-$J_2$ model at $J_2/J_1 = 0.5$. Since the space-time dimension $d+1 = 3$ separates an irrelevant from a relevant spin-spin coupling, the kagome value indicates a stable gapless spin liquid, while the square-lattice value says the coupling is relevant and leaves competing scenarios open for the intermediate regime. The claim matters because both models are long-standing testbeds where the existence and stability of spin liquids have resisted exact methods; the paper reaches $10^{-8}$ energy precision in the unfrustrated benchmark, identifies isotropic plaquette order in the VBS region, and attributes the new reach to parallelization and randomized low-rank contraction techniques.

What carries the argument

The load-bearing mechanism is finite-correlation-length scaling applied to infinite projected entangled-pair states (iPEPS). At finite bond dimension $D$ and environment bond dimension $\chi$, the transfer matrix of the infinite tensor network has a finite leading-eigenvalue gap, so correlations decay over a length $\xi(D,\chi)$ that acts as a soft infrared cutoff; for gapless phases this cutoff plays the role that system size plays in finite-size scaling. The paper uses two scaling laws, $E(D,\chi) = E_\infty + e'/\xi^3$ for the energy density and $m^2(D,\chi) \propto \xi^{-2\Delta m}$ for the squared staggered magnetization, and the stability test is $2\Delta m > 3$ (coupling irrelevant) versus $2\Delta m < 3$ (relevant). The numerical machinery carrying this out is a $\pi$-mirror-symmetric $2\times2$ ansatz built from a single independent tensor, VUMPS environment determination, a randomized QB loop contraction that lowers the next-nearest-neighbor coupling cost from $O(D^6\chi^3)$ to $O(D^4\chi^3)$, and a two-level parallelization whose strong scaling is near-ideal, reaching bond dimensions $D = 8$–$16$ with environments up to $\chi = 2048$.

Extended reading notes

Core claim

The paper's central quantitative claim is that finite-correlation-length scaling of the squared local magnetization, $m^2(D,\chi) \propto [\xi(D,\chi)]^{-2\Delta m}$, provides a proxy for the infrared scaling dimension of the spin-spin coupling, and that the fitted exponents split the two models: $2\Delta m = 3.95(11)$ for kagome, $2\Delta m = 1.020(16)$ for the square lattice at $J_2/J_1 = 0.5$. A coupling whose scaling dimension exceeds the space-time dimension $d+1 = 3$ is irrelevant, so the kagome fit supports a stable gapless U(1) Dirac spin liquid, whereas the square-lattice fit makes the spin-spin coupling relevant and the gapless state infrared-unstable; the authors state their conclusion as a 'possibly stable' gapless spin liquid on kagome and 'competing scenarios' for the square-lattice intermediate regime. In support, the paper reports oscillatory valence-bond susceptibilities coexisting with antiferromagnetic ones on kagome, isotropic plaquette order in the square-lattice VBS region for $J_2/J_1 \ge 0.52$, and an energy extrapolation linear in $1/\xi^3$ at $J_2/J_1 = 0.5$, the signature of linear dispersion, with the unfrustrated benchmark agreeing with quantum Monte Carlo to eight significant digits.

Load-bearing premise

The load-bearing premise is that the slope of $\log m^2$ versus $\log(1/\xi)$ is twice the infrared scaling dimension of the spin operator, so that crossing the space-time dimension 3 decides whether the spin-spin coupling is irrelevant and the gapless phase stable; the authors themselves label this identification a proxy. If that proxy or the finite-correlation-length scaling form fails, the kagome stability conclusion does not follow, and the kagome valence-bond susceptibility additionally rests on only three bond dimensions, $D = 8, 9, 10$.

Editorial extensions

If this is right

  • If $2\Delta m = 3.95(11)$ is correct, the kagome Heisenberg ground state is a gapless U(1) spin liquid whose spin correlations decay almost as fast as free Dirac-fermion bilinears ($2\Delta m$ near 4), with the $4\pi$-flux monopole either irrelevant or suppressed by a small prefactor.
  • The square-lattice $J_1$-$J_2$ model at $J_2/J_1 = 0.5$ is not stable under the same criterion ($2\Delta m = 1.020(16) < 3$), so the observed linear dispersion is compatible with proximity to the Néel–spin-liquid critical point or a direct transition, which is why the intermediate regime stays open.
  • The valence-bond-solid region of the square-lattice model is isotropic plaquette order for $J_2/J_1 \ge 0.52$, narrowing the symmetric region on the VBS side relative to structure-factor-based estimates.
  • The randomized loop contraction reveals an approximate low-rank structure in iPEPS correlation contractions that the paper suggests may be general, making longer-range couplings affordable at essentially the same cost as nearest-neighbor ones.
  • Because temperature and measurement resolution act as infrared cutoffs in experiments, iPEPS data at matched cutoff values can guide spin-liquid identification in scattering and transport measurements.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial extension: sweeping the same $2\Delta m$ versus 3 test across $J_2/J_1$ in the square-lattice model would turn the proxy into a phase-diagram diagnostic — a crossing above 3 would pinpoint a stable gapless regime, while its absence would favor a direct (deconfined) transition.
  • Editorial extension: the kagome exponent's proximity to 4 suggests the infrared spin response is almost free-fermion-like; if so, neutron-scattering continua in kagome candidate materials should show the same algebraic decay, giving an experimental target for the $\xi^{-2\Delta m}$ scaling seen numerically.
  • Editorial extension: the valence-bond susceptibility margin is thin ($\Delta_b = 1.75(25)$ against the $\Delta_b > 1.5$ criterion, within one error bar), and the paper itself calls for larger $D$; the next bond-dimension increment may decide whether the oscillatory valence-bond correlations are a true coexisting susceptibility or a precursor of weak order.
Original abstract

We extend the scalability of variational infinite two-dimensional tensor network methods with application to investigate frustrated magnetism of the square-lattice $J_1$-$J_2$ model and kagome-lattice Heisenberg model. In the sign-problem-free limit of the square-lattice model, we achieve high ground-state energy precision comparable to quantum Monte Carlo methods. With no obstruction beyond this limit, our method obtains direct estimates of order parameters. For the square-lattice $J_1$-$J_2$ model, we determine isotropic plaquette order in the valence-bond solid region. For the model near $J_2/J_1=0.5$, and the kagome lattice Heisenberg model, our data indicate linear-dispersion gapless states. To probe the candidate spin-liquid nature of the highly frustrated cases, we go beyond inferring the absence of symmetry-breaking orders and analyze susceptibilities that characterize gapless spin liquids. By estimating scaling dimensions, our data indicate a stable gapless spin liquid on the kagome lattice, and competing scenarios for the intermediate regime of the square-lattice model. Our developments include parallelization and randomized low-rank methods and largely expand the frontier of high-precision variational numerics for quantum phases of matter.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-11-09: 1 request, 6,362 input tokens.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.69 · human participants 0.01 · controlled comparison 0.15 · randomized assignment 0.02 · states quantitative result 0.31 · compares with prior results 0.90 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

We extend the scalability of variational infinite two-dimensional tensor network methods with application to investigate frustrated magnetism of the square-lattice $J_1$-$J_2$ model and kagome-lattice Heisenberg model. In the sign-problem-free limit of the square-lattice model, we achieve high ground-state energy precision comparable to quantum Monte Carlo methods. With no obstruction beyond this limit, our method obtains direct estimates of order parameters. For the square-lattice $J_1$-$J_2$ model, we determine isotropic plaquette order in the valence-bond solid region. For the model near $J_2/J_1=0.5$, and the kagome lattice Heisenberg model, our data indicate linear-dispersion gapless states. To probe the candidate spin-liquid nature of the highly frustrated cases, we go beyond inferring the absence of symmetry-breaking orders and analyze susceptibilities that characterize gapless spin liquids. By estimating scaling dimensions, our data indicate a stable gapless spin liquid on the kagome lattice, and competing scenarios for the intermediate regime of the square-lattice model. Our developments include parallelization and randomized low-rank methods and largely expand the frontier of high-precision variational numerics for quantum phases of matter.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We extend the scalability of variational infinite two-dimensional tensor network methods with application to investigate frustrated magnetism of the square-lattice $J_1$-$J_2$ model and kagome-lattice Heisenberg model. contribution 0.97 0.93 0.09 asserted 0.99
s1 In the sign-problem-free limit of the square-lattice model, we achieve high ground-state energy precision comparable to quantum Monte Carlo methods. result 0.96 0.95 0.09 asserted 1.00
s2 With no obstruction beyond this limit, our method obtains direct estimates of order parameters. method 0.39 0.85 0.06 asserted 0.99
s3 For the square-lattice $J_1$-$J_2$ model, we determine isotropic plaquette order in the valence-bond solid region. result 1.00 0.96 0.10 asserted 1.00
s4 For the model near $J_2/J_1=0.5$, and the kagome lattice Heisenberg model, our data indicate linear-dispersion gapless states. result 1.00 0.94 0.18 hedged possibility 0.99
s5 To probe the candidate spin-liquid nature of the highly frustrated cases, we go beyond inferring the absence of symmetry-breaking orders and analyze susceptibilities that characterize gapless spin liquids. contribution 0.58 0.94 0.08 asserted 0.77
s6 By estimating scaling dimensions, our data indicate a stable gapless spin liquid on the kagome lattice, and competing scenarios for the intermediate regime of the square-lattice model. result 0.69 0.95 0.41 hedged possibility 0.98
s7 Our developments include parallelization and randomized low-rank methods and largely expand the frontier of high-precision variational numerics for quantum phases of matter. contribution 0.91 0.93 0.12 asserted 1.00

What kind of paper the classifier reports this is: simulation study 0.63 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.