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
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
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 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.
The abstract, sentence by sentence
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
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| 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 was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
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.