pith. sign in

arxiv: 2605.18982 · v1 · pith:SY2NYQB3new · submitted 2026-05-18 · ❄️ cond-mat.str-el

Magnetic phases in the J₁-J₂ antiferromagnetic XY model on the honeycomb lattice

I. V. Lukin , M. O. Luhanko , Yu. V. Slyusarenko , A. G. Sotnikov This is my paper

Pith reviewed 2026-05-20 08:10 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords J1-J2 XY modelhoneycomb latticeantiferromagnetphase diagramNéel orderspiral phasecorner transfer matrixiPEPS
0
0 comments X

The pith

The J1-J2 XY antiferromagnet on the honeycomb lattice supports Néel, Ising, collinear, and incommensurate spiral phases with transitions set by the frustration ratio J2/J1.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper maps the ground-state phase diagram of the antiferromagnetic XY model with nearest-neighbor J1 and next-nearest-neighbor J2 couplings on the honeycomb lattice. It employs numerical tensor-network methods to track how magnetic order changes with increasing J2/J1 and to compare energies of different candidate states. The work finds that a collinear phase is stable over dimerized alternatives in the frustrated regime and that this collinear phase ends in a continuous transition to an incommensurate spiral.

Core claim

The authors establish that the model realizes Néel order at small J2/J1, an Ising phase, then a collinear phase whose energy lies below that of dimerized states; the collinear phase undergoes a second-order transition to an incommensurate spiral phase at larger J2/J1.

What carries the argument

Corner transfer matrix renormalization group algorithm with a two-site unit cell together with the infinite spiral projected entangled pair states ansatz used to compute variational energies and order parameters.

If this is right

  • The located phase boundaries supply reference points for comparing with other frustrated two-dimensional magnets.
  • The continuous nature of the collinear-to-spiral transition implies that the spiral wave vector evolves without an abrupt jump.
  • Dimerization is energetically disfavored relative to collinear order throughout the regime of competing interactions.
  • The same computational framework can be reused to study related models with altered spin anisotropy or lattice geometry.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The phase sequence may help interpret magnetic ordering in layered van der Waals materials whose spins approximate the XY limit on honeycomb lattices.
  • Increasing the bond dimension in the same ansatz would provide a direct numerical test of whether the transition remains continuous.
  • The suppression of dimer order relative to Heisenberg versions of the model points to the role of XY anisotropy in stabilizing collinear over valence-bond states.

Load-bearing premise

The numerical renormalization procedure and the chosen ansatz converge to the true ground state without significant truncation or unit-cell bias across the full range of J2/J1.

What would settle it

Exact diagonalization or density-matrix renormalization group on cylinders showing either a dimerized state lower in energy than the collinear state or a discontinuous jump in the ordering wave vector at the putative second-order transition.

Figures

Figures reproduced from arXiv: 2605.18982 by A. G. Sotnikov, I. V. Lukin, M. O. Luhanko, Yu. V. Slyusarenko.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Honeycomb lattice with the two-site unit cell in [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) Honeycomb PEPS state with the two-site unit cell [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Definitions of different boundary tensors appearing [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Illustration of one step of the CTMRG update: (a) [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Computation of correlation function with CTM ten [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Ground-state energy (a) and staggered magnetiza [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. The behavior of the energy (a) and position [PITH_FULL_IMAGE:figures/full_fig_p007_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Phase diagram of the [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗
read the original abstract

We study ground-state properties and phase diagram of the $J_{1}$-$J_{2}$ antiferromagnetic XY model on the honeycomb lattice by means of the developed corner transfer matrix renormalization group algorithm with the two-site unit cell and the infinite spiral projected entangled pair states ansatz. We identify the main phases: N\'{e}el, Ising, collinear, and incommensurate spiral phases, as well as the transitions between them, as functions of the ratio $J_{2}/J_{1}$. In the regime of competing types of ordering, we show that the energies of the dimerized states are systematically higher than the energies in the collinear phase. This collinear phase transforms to the incommensurate spiral phase through the second-order phase transition upon a further increase of $J_2/J_1$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The manuscript investigates the ground-state phase diagram of the J₁-J₂ antiferromagnetic XY model on the honeycomb lattice. Employing a corner transfer matrix renormalization group (CTMRG) algorithm with a two-site unit cell together with an infinite spiral projected entangled pair states (iPEPS) ansatz, the authors identify Néel, Ising, collinear, and incommensurate spiral phases, map the transitions as functions of J₂/J₁, demonstrate that dimerized states lie higher in energy than the collinear phase in competing-order regimes, and characterize the collinear-to-incommensurate-spiral transition as second-order.

Significance. If the numerical results prove robust, the work supplies a detailed phase diagram for a frustrated XY antiferromagnet on the honeycomb lattice, clarifying the stability of collinear order relative to dimerized states and the order of the transition into spiral order. The direct minimization via tensor-network methods on infinite systems, with no free parameters or fitted quantities, constitutes a methodological strength that supports the reported phase boundaries and transition characters.

major comments (1)
  1. The central claims—phase boundaries, the systematic energy advantage of the collinear phase over dimerized states, and the second-order character of the collinear-to-spiral transition—rest on the convergence of the two-site CTMRG and infinite spiral iPEPS ansatz to the true ground state for all J₂/J₁. The manuscript provides no bond-dimension extrapolations, truncation-error bounds, or comparisons against larger unit cells or exact diagonalization benchmarks, leaving open the possibility of ansatz bias or truncation artifacts in the frustrated regime where Néel, collinear, spiral, and dimerized tendencies compete.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive feedback. We address the single major comment below and have prepared revisions to strengthen the presentation of numerical convergence.

read point-by-point responses

  1. Referee: The central claims—phase boundaries, the systematic energy advantage of the collinear phase over dimerized states, and the second-order character of the collinear-to-spiral transition—rest on the convergence of the two-site CTMRG and infinite spiral iPEPS ansatz to the true ground state for all J₂/J₁. The manuscript provides no bond-dimension extrapolations, truncation-error bounds, or comparisons against larger unit cells or exact diagonalization benchmarks, leaving open the possibility of ansatz bias or truncation artifacts in the frustrated regime where Néel, collinear, spiral, and dimerized tendencies compete.

    Authors: We agree that explicit documentation of convergence is essential for claims in the frustrated regime. In the revised manuscript we will add (i) bond-dimension extrapolations of the energy and order parameters for representative values of J₂/J₁ in each phase, (ii) truncation-error estimates obtained from the CTMRG and iPEPS algorithms, and (iii) a direct comparison of our infinite-system results with exact-diagonalization data on finite clusters up to 24 sites. These additions will quantify the residual truncation error and demonstrate that the reported phase boundaries and transition order remain stable within the stated precision. We note that the two-site unit cell is the natural choice for the honeycomb lattice and that the spiral iPEPS ansatz is constructed to accommodate incommensurate wave vectors; we will nevertheless clarify why larger supercells are not required for the phases studied. revision: yes

Circularity Check

0 steps flagged

Numerical tensor-network minimization of Hamiltonian yields self-contained phase diagram

full rationale

The work computes ground-state energies and order parameters directly via CTMRG (two-site cell) and infinite spiral iPEPS variational optimization of the J1-J2 XY Hamiltonian for varying J2/J1. Phase boundaries, dimer vs. collinear energy comparisons, and the collinear-to-spiral transition order are outputs of this minimization; no parameter is fitted to a target observable and then re-used as a prediction, no self-citation supplies a uniqueness theorem that forces the ansatz, and no quantity is defined in terms of itself. The derivation chain therefore remains independent of its own results.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The claim rests on the domain assumption that the honeycomb lattice with nearest-neighbor J1 and next-nearest-neighbor J2 antiferromagnetic XY couplings is the correct microscopic model, plus the technical assumption that the chosen tensor-network ansatzes faithfully represent the ground state. No free parameters or new entities are introduced.

axioms (1)
  • domain assumption The honeycomb lattice with nearest-neighbor J1 and next-nearest-neighbor J2 antiferromagnetic XY interactions defines the Hamiltonian under study.
    This is the standard microscopic model stated in the abstract.

pith-pipeline@v0.9.0 · 5699 in / 1370 out tokens · 60347 ms · 2026-05-20T08:10:05.214780+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages · 1 internal anchor

  1. [1]

    Lacroix, P

    C. Lacroix, P. Mendels, and F. Mila,Introduction to frustrated magnetism: materials, experiments, theory (Springer Science & Business Media, 2011)

    work page 2011

  2. [2]

    C. N. Varney, K. Sun, V. Galitski, and M. Rigol, Kalei- doscope of exotic quantum phases in a frustratedXY model, Phys. Rev. Lett.107, 077201 (2011)

    work page 2011

  3. [3]

    Z. Zhu, D. A. Huse, and S. R. White, Unexpectedz- direction Ising antiferromagnetic order in a frustrated spin-1/2J 1 −J 2 XYmodel on the honeycomb lattice, Phys. Rev. Lett.111, 257201 (2013)

    work page 2013

  4. [4]

    Di Ciolo, J

    A. Di Ciolo, J. Carrasquilla, F. Becca, M. Rigol, and V. Galitski, Spiral antiferromagnets beyond the spin- wave approximation: FrustratedXYand Heisenberg models on the honeycomb lattice, Phys. Rev. B89, 094413 (2014)

    work page 2014

  5. [5]

    Zhu and S

    Z. Zhu and S. R. White, Quantum phases of the frus- tratedXYmodels on the honeycomb lattice, Modern Physics Letters B28, 1430016 (2014)

    work page 2014

  6. [6]

    R. F. Bishop, P. H. Y. Li, and C. E. Campbell, Frustrated spin- 1 2 J1–J2 isotropicXYmodel on the honeycomb lat- tice, Phys. Rev. B89, 214413 (2014)

    work page 2014

  7. [7]

    Huang, X.-Y

    Y. Huang, X.-Y. Dong, D. N. Sheng, and C. S. Ting, Quantum phase diagram and chiral spin liquid in the ex- tended spin- 1 2 honeycomb XY model, Phys. Rev. B103, L041108 (2021)

    work page 2021

  8. [8]

    Oitmaa and R

    J. Oitmaa and R. R. P. Singh, Phase diagram of the frustrated quantum-XYmodel on the honeycomb lattice studied by series expansions: Evidence for proximity to a bicritical point, Phys. Rev. B89, 104423 (2014)

    work page 2014

  9. [9]

    Carrasquilla, A

    J. Carrasquilla, A. D. Ciolo, F. Becca, V. Galitski, and M. Rigol, Nature of the phases in the frustrated XYmodel on the honeycomb lattice, Phys. Rev. B88, 241109(R) (2013)

    work page 2013

  10. [10]

    Satoori, S

    S. Satoori, S. Mahdavifar, and J. Vahedi, Quantum cor- relations in the frustrated XY model on the honeycomb lattice, Sci. Rep.13, 16034 (2023)

    work page 2023

  11. [11]

    A. L. Chernyshev, Quantum stabilization of unexpected ordered phases on the honeycomb lattice, Phys. Rev. B 112, 174413 (2025)

    work page 2025

  12. [12]

    Renormalization algorithms for Quantum-Many Body Systems in two and higher dimensions

    F. Verstraete and J. I. Cirac, Renormalization algorithms for Quantum-Many Body Systems in two and higher dimensions, arXiv preprint cond-mat/0407066 (2004), arXiv:cond-mat/0407066

    work page internal anchor Pith review Pith/arXiv arXiv 2004

  13. [13]

    Verstraete, V

    F. Verstraete, V. Murg, and J. Cirac, Matrix product states, projected entangled pair states, and variational renormalization group methods for quantum spin sys- tems, Advances in Physics57, 143 (2008)

    work page 2008

  14. [14]

    Or´ us, A practical introduction to tensor networks: Matrix product states and projected entangled pair states, Annals of Physics349, 117 (2014)

    R. Or´ us, A practical introduction to tensor networks: Matrix product states and projected entangled pair states, Annals of Physics349, 117 (2014)

    work page 2014

  15. [15]

    J. I. Cirac, D. P´ erez-Garc´ ıa, N. Schuch, and F. Ver- straete, Matrix product states and projected entangled pair states: Concepts, symmetries, theorems, Rev. Mod. Phys.93, 045003 (2021)

    work page 2021

  16. [16]

    Or´ us, Tensor networks for complex quantum systems, Nature Reviews Physics1, 538 (2019)

    R. Or´ us, Tensor networks for complex quantum systems, Nature Reviews Physics1, 538 (2019)

    work page 2019

  17. [17]

    M. C. Ba˜ nuls, Tensor network algorithms: a route map, Annual Review of Condensed Matter Physics14, 173 (2023)

    work page 2023

  18. [18]

    Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics 326, 96 (2011)

    U. Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics 326, 96 (2011)

    work page 2011

  19. [19]

    S. R. White, Density-matrix algorithms for quantum renormalization groups, Phys. Rev. B48, 10345 (1993)

    work page 1993

  20. [20]

    Eisert, M

    J. Eisert, M. Cramer, and M. B. Plenio, Colloquium: Area laws for the entanglement entropy, Rev. Mod. Phys. 82, 277 (2010)

    work page 2010

  21. [21]

    Hasik and P

    J. Hasik and P. Corboz, Incommensurate order with translationally invariant projected entangled-pair states: Spiral states and quantum spin liquid on the anisotropic triangular lattice, Phys. Rev. Lett.133, 176502 (2024)

    work page 2024

  22. [22]

    Schmoll, J

    P. Schmoll, J. Naumann, E. L. Weerda, J. Eisert, and Y. Iqbal, Bathing in a sea of candidate quantum spin liquids: From the gapless ruby to the gapped maple-leaf lattice (2025), arXiv:2407.07145 [cond-mat.str-el]

    work page arXiv 2025

  23. [23]

    Corboz, Variational optimization with infinite pro- jected entangled-pair states, Phys

    P. Corboz, Variational optimization with infinite pro- jected entangled-pair states, Phys. Rev. B94, 035133 (2016)

    work page 2016

  24. [24]

    Vanderstraeten, J

    L. Vanderstraeten, J. Haegeman, P. Corboz, and F. Ver- straete, Gradient methods for variational optimization of projected entangled-pair states, Phys. Rev. B94, 155123 (2016)

    work page 2016

  25. [25]

    Liao, J.-G

    H.-J. Liao, J.-G. Liu, L. Wang, and T. Xiang, Differ- entiable programming tensor networks, Phys. Rev. X9, 031041 (2019)

    work page 2019

  26. [26]

    Naumann, E

    J. Naumann, E. L. Weerda, M. Rizzi, J. Eisert, and P. Schmoll, An introduction to infinite projected entangled-pair state methods for variational ground state simulations using automatic differentiation, SciPost Phys. Lect. Notes86, 1 (2024). 9

    work page 2024

  27. [27]

    Nishino and K

    T. Nishino and K. Okunishi, Corner transfer matrix renormalization group method, J. Phys. Soc. Jpn.65, 891 (1996)

    work page 1996

  28. [28]

    Nishino and K

    T. Nishino and K. Okunishi, Corner transfer matrix al- gorithm for classical renormalization group, J. Phys. Soc. Jpn.66, 3040 (1997)

    work page 1997

  29. [29]

    Or´ us and G

    R. Or´ us and G. Vidal, Simulation of two-dimensional quantum systems on an infinite lattice revisited: Cor- ner transfer matrix for tensor contraction, Phys. Rev. B 80, 094403 (2009)

    work page 2009

  30. [30]

    I. V. Lukin and A. G. Sotnikov, Variational optimiza- tion of tensor-network states with the honeycomb-lattice corner transfer matrix, Phys. Rev. B107, 054424 (2023)

    work page 2023

  31. [31]

    Nyckees, A

    S. Nyckees, A. Rufino, F. Mila, and J. Colbois, Critical line of the triangular Ising antiferromagnet in a field from aC 3-symmetric corner transfer matrix algorithm, Phys. Rev. E108, 064132 (2023)

    work page 2023

  32. [32]

    Gendiar, R

    A. Gendiar, R. Krcmar, S. Andergassen, M. Daniˇ ska, and T. Nishino, Weak correlation effects in the Ising model on triangular-tiled hyperbolic lattices, Phys. Rev. E86, 021105 (2012)

    work page 2012

  33. [33]

    Naumann, J

    J. Naumann, J. Eisert, and P. Schmoll, Variational opti- mization of projected entangled-pair states on the trian- gular lattice, Phys. Rev. B113, 045117 (2026)

    work page 2026

  34. [34]

    Yang and P

    Q. Yang and P. Corboz, Efficient iPEPS simulation on the honeycomb lattice via QR-based CTMRG, Phys. Rev. B113, 085109 (2026)

    work page 2026

  35. [35]

    Nyckees, P

    S. Nyckees, P. Ghosh, and F. Mila, Tensor-network study of the ground state of the maple-leaf Heisen- berg antiferromagnet, Zeitschrift f¨ ur Naturforschung A doi:10.1515/zna-2025-0409 (2026)

    work page doi:10.1515/zna-2025-0409 2025

  36. [36]

    Ghosh and F

    P. Ghosh and F. Mila, Simplex crystal ground state and magnetization plateaus in the spin-1/2 Heisenberg model on the ruby lattice (2025), arXiv:2512.14173 [cond- mat.str-el]

    work page arXiv 2025

  37. [37]

    Ghosh, J

    P. Ghosh, J. A. Koziol, S. Nyckees, K. P. Schmidt, and F. Mila, Symmetry breaking and competing valence bond states in the star lattice Heisenberg antiferromagnet, Phys. Rev. B112, 144423 (2025)

    work page 2025

  38. [38]

    Zhang, Q

    X.-Y. Zhang, Q. Yang, P. Corboz, J. Haegeman, and Y. He, Topological and trivial valence-bond orders in higher-spin Kitaev models (2025), arXiv:2511.07415 [cond-mat.str-el]

    work page arXiv 2025

  39. [39]

    Levin and C

    M. Levin and C. P. Nave, Tensor renormalization group approach to two-dimensional classical lattice models, Phys. Rev. Lett.99, 120601 (2007)

    work page 2007

  40. [40]

    H. C. Jiang, Z. Y. Weng, and T. Xiang, Accurate de- termination of tensor network state of quantum lattice models in two dimensions, Phys. Rev. Lett.101, 090603 (2008)

    work page 2008

  41. [41]

    Jordan, R

    J. Jordan, R. Or´ us, G. Vidal, F. Verstraete, and J. I. Cirac, Classical simulation of infinite-size quantum lat- tice systems in two spatial dimensions, Phys. Rev. Lett. 101, 250602 (2008)

    work page 2008

  42. [42]

    I. V. Lukin and A. G. Sotnikov, Corner transfer ma- trix renormalization group approach in the zoo of Archimedean lattices, Phys. Rev. E109, 045305 (2024)

    work page 2024

  43. [43]

    Corboz, T

    P. Corboz, T. M. Rice, and M. Troyer, Competing states in thet-Jmodel: Uniformd-wave state versus stripe state, Phys. Rev. Lett.113, 046402 (2014)

    work page 2014

  44. [44]

    Paszke, S

    A. Paszke, S. Gross, F. Massa, A. Lerer, J. Bradbury, G. Chanan, T. Killeen, Z. Lin, N. Gimelshein, L. Antiga, A. Desmaison, A. K¨ opf, E. Yang, Z. DeVito, M. Raison, A. Tejani, S. Chilamkurthy, B. Steiner, L. Fang, J. Bai, and S. Chintala, PyTorch: An imperative style, high- performance deep learning library, inAdvances in Neural Information Processing S...

    work page 2019

  45. [45]

    Hasik, D

    J. Hasik, D. Poilblanc, and F. Becca, Investigation of the N´ eel phase of the frustrated Heisenberg antiferromag- net by differentiable symmetric tensor networks, SciPost Phys.10, 012 (2021)

    work page 2021

  46. [46]

    Hasik and G

    J. Hasik and G. B. Mbeng, peps-torch: A li- brary for simulation of strongly correlated systems using iPEPS tensor networks,https://github.com/ jurajHasik/peps-torch(2021)

    work page 2021

  47. [47]

    Jiang, S

    S. Jiang, S. R. White, and A. L. Chernyshev, Quan- tum phases in the honeycomb-latticeJ 1–J3 ferro- antiferromagnetic model, Phys. Rev. B108, L180406 (2023)

    work page 2023

  48. [48]

    J. B. Fouet, P. Sindzingre, and C. Lhuillier, An investiga- tion of the quantumJ 1–J2–J3 model on the honeycomb lattice, Eur. Phys. J. B20, 241 (2001)

    work page 2001

  49. [49]

    Oitmaa and R

    J. Oitmaa and R. R. P. Singh, Phase diagram of theJ 1 − J2 −J 3 Heisenberg model on the honeycomb lattice: A series expansion study, Phys. Rev. B84, 094424 (2011)

    work page 2011