pith. sign in

arxiv: 2604.14767 · v1 · submitted 2026-04-16 · ❄️ cond-mat.str-el

Nontrivial three-sublattice magnetization in the easy-axis spin-1/2 XXZ antiferromagnet on the triangular lattice

Masahiro Kadosawa , Masaaki Nakamura , Yukinori Ohta , Satoshi Nishimoto This is my paper

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

classification ❄️ cond-mat.str-el
keywords XXZ antiferromagnettriangular latticeeasy-axis anisotropythree-sublattice orderquantum fluctuationsDMRGY stateIsing manifold
0
0 comments X

The pith

Quantum fluctuations select a nontrivial three-sublattice state in the easy-axis XXZ antiferromagnet on the triangular lattice that persists even at infinite anisotropy.

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

The paper applies density-matrix renormalization group methods to finite clusters of the spin-1/2 XXZ antiferromagnet on the triangular lattice, mapped to one-dimensional chains via spiral boundary conditions. It tracks the evolution of three-sublattice magnetizations as the easy-axis anisotropy parameter Delta increases from the isotropic Heisenberg point. At large Delta the local moments stay close to a (2m, -m, -m) pattern whose extrapolated values are 0.41873 and 0.20832, well below classical saturation. The authors also show that this Y-state configuration has lower energy than the up-down-down state and matches unbiased thermodynamic-limit energies. This establishes that the ground state remains a quantum-selected ordered structure rather than becoming a trivial saturated collinear Ising state.

Core claim

In the easy-axis regime, the ordered moments remain close to a zero-magnetization three-sublattice structure of the form (2m,-m,-m) over a broad range of Delta. Extrapolation in 1/Delta shows that the positive sublattice moment stays well below the classical saturation value 1/2, approaching 0.41873 as Delta to infinity, while the magnitude of the negative sublattice moment approaches 0.20832. The Y state is favored at zero field, and independent thermodynamic-limit energy calculations yield an energy consistent with the Y-state solution. These results show that the easy-axis ground state does not simply cross over to a trivially saturated collinear Ising state, but instead remains a trivial

What carries the argument

Three-sublattice moments extracted from symmetry-broken local magnetization profiles on spiral-boundary DMRG clusters, extrapolated in 1/Delta to the infinite-anisotropy limit.

Where Pith is reading between the lines

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

  • This selection from the degenerate manifold may operate in other anisotropic frustrated spin systems.
  • Neutron scattering or NMR in candidate materials could detect the specific (2m, -m, -m) moment imbalance.
  • Finite-temperature or field-dependent extensions might show how the selected order evolves or melts.

Load-bearing premise

The extrapolation of sublattice moments in 1/Delta to the infinite-anisotropy limit accurately captures the thermodynamic-limit behavior without significant finite-size or boundary-condition artifacts.

What would settle it

A direct calculation or measurement of the sublattice moments reaching exactly the classical saturation values of 0.5 and 0 in the thermodynamic limit at very large anisotropy would falsify the nontrivial order claim.

Figures

Figures reproduced from arXiv: 2604.14767 by Masaaki Nakamura, Masahiro Kadosawa, Satoshi Nishimoto, Yukinori Ohta.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic illustrations of the magnetic structures relevant to [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) Schematic illustration of the [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Anisotropy dependence of the positive and negative sublat [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Local spin moment [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: shows that both 2m1 and m2 vary monotonically with 1/∆. Their dependence is approximately linear over the intermediate range shown, whereas the slope gradually de￾creases toward the Ising limit 1/∆ → 0, indicating a tendency toward saturation. Extrapolating the data to 1/∆ → 0, we ob￾tain finite values for both quantities. In particular, the extrap￾olated value of the positive sublattice moment is estimate… view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Scaled energies of the Y state and UDD state as functions of [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Finite-size scaling of the net magnetization per site [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Comparison of the sublattice magnetization [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
read the original abstract

We investigate the ground-state magnetic structure of the spin-$1/2$ XXZ antiferromagnet on the triangular lattice in the easy-axis regime using the density-matrix renormalization group. By applying spiral boundary conditions, we exactly map finite $L\times L$ clusters onto one-dimensional chains while avoiding the spatial anisotropy inherent in cylindrical geometries. From symmetry-broken local magnetization profiles, we extract the three-sublattice moments and track their evolution with anisotropy. At the isotropic point, we obtain a positive sublattice moment of $0.21671$, consistent with previous numerical estimates. In the easy-axis regime, the ordered moments remain close to a zero-magnetization three-sublattice structure of the form $(2m,-m,-m)$ over a broad range of $\Delta$. Extrapolation in $1/\Delta$ shows that the positive sublattice moment stays well below the classical saturation value $1/2$, approaching $0.41873$ as $\Delta\to\infty$, while the magnitude of the negative sublattice moment approaches $0.20832$. We further compare the energies of the Y state and the up-down-down state and find that the Y state is favored at zero field. Independent thermodynamic-limit energy calculations, performed without assuming any particular ordered pattern, yield an energy consistent with the Y-state solution. These results show that the easy-axis ground state does not simply cross over to a trivially saturated collinear Ising state, but instead remains a nontrivial three-sublattice ordered state selected from the macroscopically degenerate Ising manifold by quantum fluctuations.

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

2 major / 2 minor

Summary. The manuscript investigates the ground state of the spin-1/2 XXZ antiferromagnet on the triangular lattice in the easy-axis regime (Δ>1) via DMRG. Finite L×L clusters are mapped to 1D chains using spiral boundary conditions. From symmetry-broken local magnetization profiles the authors extract three-sublattice moments, extrapolate them linearly in 1/Δ, and report limiting values 0.41873 and 0.20832 (below classical saturation 0.5 and 0.25). They conclude that the ground state remains a nontrivial zero-magnetization Y-state selected by quantum fluctuations rather than crossing over to a saturated collinear Ising configuration. Supporting evidence includes energy comparisons between the Y-state and up-down-down state plus independent thermodynamic-limit energy calculations performed without assuming a particular ordered pattern.

Significance. If the reported extrapolations accurately represent the thermodynamic limit, the work clarifies the persistence of quantum order-by-disorder selection within the macroscopically degenerate Ising manifold on the triangular lattice. The spiral-boundary-condition mapping that enables controlled DMRG on 2D clusters and the independent energy calculations that do not presuppose order are methodological strengths that strengthen the numerical evidence.

major comments (2)
  1. [Abstract and extrapolation section] Abstract and the section describing the 1/Δ extrapolation: the sublattice moments are extrapolated only in 1/Δ on finite-L clusters; no prior L→∞ scaling at fixed large Δ is reported. Because the central claim that the moments remain below saturation (0.41873 and 0.20832) rests on these intercepts being the true thermodynamic-limit values, the absence of explicit finite-size analysis leaves open the possibility that L-dependent corrections alter the intercepts and weaken the evidence against a saturated collinear state.
  2. [Energy comparison and thermodynamic-limit calculations] Section on energy comparisons and independent thermodynamic-limit calculations: while the Y-state is shown to be lower in energy than the up-down-down state and independent calculations are stated to be consistent with the Y-state, the manuscript does not specify the range of system sizes, boundary conditions, or extrapolation procedure used for those thermodynamic-limit energies. This information is required to assess convergence and to confirm that the energy comparison is not itself affected by the same finite-size issues raised for the moments.
minor comments (2)
  1. [Abstract] The abstract provides no information on the DMRG bond dimensions, truncation errors, or convergence checks employed for the mapped 1D chains. These technical details should be added (or referenced to a methods section) so that readers can evaluate the numerical accuracy of the reported moments.
  2. Figures displaying local magnetization profiles should explicitly label the three sublattices and indicate the spiral boundary conditions used; this would improve clarity when readers compare the extracted moments to the classical (2m,−m,−m) structure.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We are pleased that the referee recognizes the significance of our results and the strengths of our approach. We respond to each major comment below, indicating the revisions we will make to address the concerns.

read point-by-point responses

  1. Referee: [Abstract and extrapolation section] Abstract and the section describing the 1/Δ extrapolation: the sublattice moments are extrapolated only in 1/Δ on finite-L clusters; no prior L→∞ scaling at fixed large Δ is reported. Because the central claim that the moments remain below saturation (0.41873 and 0.20832) rests on these intercepts being the true thermodynamic-limit values, the absence of explicit finite-size analysis leaves open the possibility that L-dependent corrections alter the intercepts and weaken the evidence against a saturated collinear state.

    Authors: We acknowledge the validity of this observation. Our primary extrapolation was performed in 1/Δ because the spiral boundary conditions provide excellent convergence for the local magnetizations even on moderate L. To strengthen our claim, we will add in the revised manuscript an analysis of the L-dependence at several fixed large values of Δ, including extrapolations to L→∞. This will demonstrate that finite-size corrections do not alter the conclusion that the moments remain below saturation in the thermodynamic limit. We will update the abstract and the relevant section to reflect this additional evidence. revision: yes

  2. Referee: [Energy comparison and thermodynamic-limit calculations] Section on energy comparisons and independent thermodynamic-limit calculations: while the Y-state is shown to be lower in energy than the up-down-down state and independent calculations are stated to be consistent with the Y-state, the manuscript does not specify the range of system sizes, boundary conditions, or extrapolation procedure used for those thermodynamic-limit energies. This information is required to assess convergence and to confirm that the energy comparison is not itself affected by the same finite-size issues raised for the moments.

    Authors: We agree that more details are needed for reproducibility and assessment. The thermodynamic-limit energy calculations were carried out using DMRG on cylindrical clusters with periodic boundary conditions along the short direction and open boundaries along the long direction, for widths L=4 to L=12 and lengths up to several times the width, with finite-size extrapolation in 1/L. We will include a detailed description of the system sizes, boundary conditions, and extrapolation procedure in the revised manuscript, along with the specific energy values obtained. revision: yes

Circularity Check

0 steps flagged

No significant circularity in numerical derivation chain

full rationale

The paper obtains its central results via direct DMRG simulations on finite L×L clusters mapped to 1D chains with spiral boundary conditions, extracts symmetry-broken local magnetization profiles to determine three-sublattice moments, and performs a linear extrapolation in 1/Δ to the infinite-anisotropy limit. These extrapolated values (0.41873 and 0.20832) are compared against classical saturation (0.5 and 0.25) as an independent benchmark, and an additional energy calculation without assuming any ordered pattern is shown to be consistent with the Y state. No step in the chain reduces by construction to its inputs, no parameters are fitted and then relabeled as predictions, and no load-bearing self-citations or imported uniqueness theorems are invoked. The procedure is a standard, externally falsifiable numerical workflow whose outputs are independent of the target claim.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The claim rests on the validity of DMRG as an unbiased solver for the mapped chains and on the assumption that finite-size extrapolation in 1/Δ converges to the correct thermodynamic-limit order parameters.

axioms (2)
  • domain assumption DMRG truncation error is negligible for the system sizes employed
    Standard assumption in DMRG studies; no explicit error bars or bond-dimension scaling shown in abstract
  • domain assumption Spiral boundary conditions introduce no spurious anisotropy or topological effects that alter the bulk order
    The paper states this mapping avoids cylindrical anisotropy, but the assumption is not independently verified in the provided abstract

pith-pipeline@v0.9.0 · 5603 in / 1383 out tokens · 53637 ms · 2026-05-10T10:28:37.546704+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

42 extracted references · 42 canonical work pages

  1. [1]

    Three-sublattice magnetization In contrast to bipartite lattices, where the magnetic order can be characterized by the staggered magnetization, as in our previous analyses of the square- and honeycomb-lattice cases [24, 25], the triangular-lattice XXZ model is expected to exhibit a three-sublattice magnetic structure. In the present analysis, we assume th...

    work page

  2. [2]

    III B, the DMRG calculations with pinning fields are performed in sectors of fixed total magne- tizationS z tot

    Energy As discussed in Sec. III B, the DMRG calculations with pinning fields are performed in sectors of fixed total magne- tizationS z tot. Because the spin configuration near the system edges is directly constrained by the pinning fields, the total energy of a finite system is not always the most suitable quan- tity for characterizing the bulk magnetic ...

    work page

  3. [3]

    A. P. Ramirez, Strongly geometrically frustrated magnets, An- nual Review of Materials Research24, 453 (1994)

    work page 1994

  4. [4]

    Balents, Spin liquids in frustrated magnets, Nature464, 199 (2010)

    L. Balents, Spin liquids in frustrated magnets, Nature464, 199 (2010)

    work page 2010

  5. [5]

    Lacroix, P

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

    work page 2011

  6. [6]

    Miyashita and H

    S. Miyashita and H. Kawamura, Phase transitions of anisotropic heisenberg antiferromagnets on the triangular lattice, Journal of the Physical Society of Japan54, 3385 (1985)

    work page 1985

  7. [7]

    Sellmann, X.-F

    D. Sellmann, X.-F. Zhang, and S. Eggert, Phase diagram of the antiferromagnetic xxz model on the triangular lattice, Physical Review B91, 081104 (2015)

    work page 2015

  8. [8]

    Wessel and M

    S. Wessel and M. Troyer, Supersolid hard-core bosons on the triangular lattice, Physical Review Letters95, 127205 (2005)

    work page 2005

  9. [9]

    Heidarian and K

    D. Heidarian and K. Damle, Persistent supersolid phase of hard- core bosons on the triangular lattice, Physical Review Letters 95, 127206 (2005)

    work page 2005

  10. [10]

    Boninsegni and N

    M. Boninsegni and N. Prokof’ev, Supersolid phase of hard- core bosons on a triangular lattice, Phys. Rev. Lett.95, 237204 (2005)

    work page 2005

  11. [11]

    R. G. Melko, A. Paramekanti, A. A. Burkov, A. Vishwanath, D. N. Sheng, and L. Balents, Supersolid order from disorder: Hard-core bosons on the triangular lattice, Phys. Rev. Lett.95, 127207 (2005)

    work page 2005

  12. [12]

    D. A. Huse and V . Elser, Simple variational wave functions for two-dimensional heisenberg spin-1/2 antiferromagnets, Phys. Rev. Lett.60, 2531 (1988)

    work page 1988

  13. [13]

    Bernu, C

    B. Bernu, C. Lhuillier, and L. Pierre, Signature of n ´eel order in exact spectra of quantum antiferromagnets on finite lattices, Phys. Rev. Lett.69, 2590 (1992)

    work page 1992

  14. [14]

    Bernu, P

    B. Bernu, P. Lecheminant, C. Lhuillier, and L. Pierre, Exact spectra, spin susceptibilities, and order parameter of the quan- tum heisenberg antiferromagnet on the triangular lattice, Phys. Rev. B50, 10048 (1994)

    work page 1994

  15. [15]

    Capriotti, A

    L. Capriotti, A. E. Trumper, and S. Sorella, Long-range n ´eel order in the triangular heisenberg model, Phys. Rev. Lett.82, 3899 (1999)

    work page 1999

  16. [16]

    Zheng, J

    W. Zheng, J. O. Fjæ restad, R. R. P. Singh, R. H. McKenzie, and R. Coldea, Excitation spectra of the spin-1/2 triangular-lattice heisenberg antiferromagnet, Phys. Rev. B74, 224420 (2006)

    work page 2006

  17. [17]

    S. R. White and A. L. Chernyshev, N´eel order in square and tri- angular lattice heisenberg models, Phys. Rev. Lett.99, 127004 (2007)

    work page 2007

  18. [18]

    Q. Li, H. Li, J. Zhao, H.-G. Luo, and Z. Y . Xie, Magnetization of the spin-1/2 heisenberg antiferromagnet on the triangular lat- tice, Physical Review B105, 184418 (2022)

    work page 2022

  19. [19]

    Huang, X

    J. Huang, X. Qian, and M. Qin, On the magnetization of the 120◦ order of the spin-1/2triangular lattice heisenberg model: a dmrg revisit, J. Phys.: Condens. Matter36, 185602 (2024)

    work page 2024

  20. [20]

    Zheng, J

    W. Zheng, J. O. Fjærestad, R. R. P. Singh, R. H. McKenzie, and R. Coldea, Excitation spectra of the spin- 1 2 triangular-lattice heisenberg antiferromagnet, Phys. Rev. B74, 224420 (2006)

    work page 2006

  21. [21]

    Gao, Y .-C

    Y . Gao, Y .-C. Fan, H. Li, F. Yang, X.-T. Zeng, X.-L. Sheng, R. Zhong, Y . Qi, Y . Wan, and W. Li, Spin supersolidity in nearly ideal easy-axis triangular quantum antiferromagnet na2baco (po4) 2, npj Quantum Materials7, 89 (2022)

    work page 2022

  22. [22]

    Sheng, L

    J. Sheng, L. Wang, A. Candini, W. Jiang, L. Huang, B. Xi, J. Zhao, H. Ge, N. Zhao, Y . Fu, J. Ren, J. Yang, P. Miao, X. Tong, D. Yu, S. Wang, Q. Liu, M. Kofu, R. Mole, G. Bi- asiol, D. Yu, I. A. Zaliznyak, J.-W. Mei, and L. Wu, Two- dimensional quantum universality in the spin-1/2 triangular- lattice quantum antiferromagnet na2baco(po4)2, Proceedings of t...

    work page doi:10.1073/pnas.2211193119 2022

  23. [23]

    M. Zhu, L. M. Chinellato, V . Romerio, N. Murai, S. Ohira- Kawamura, C. Balz, Z. Yan, S. Gvasaliya, Y . Kato, C. D. Batista, and A. Zheludev, Wannier states and spin supersolid physics in the triangular antiferromagnet k2co(seo3)2, npj Quantum Materials10, 74 (2025)

    work page 2025

  24. [24]

    P. W. Anderson, An approximate quantum theory of the antifer- romagnetic ground state, Phys. Rev.86, 694 (1952)

    work page 1952

  25. [25]

    Sachdev,Quantum phases of matter(Cambridge University Press, 2023)

    S. Sachdev,Quantum phases of matter(Cambridge University Press, 2023)

    work page 2023

  26. [26]

    Kadosawa, M

    M. Kadosawa, M. Nakamura, Y . Ohta, and S. Nishimoto, Study of staggered magnetization in the spin-s square-lattice heisen- berg model using spiral boundary conditions, Journal of the Physical Society of Japan92, 023701 (2023)

    work page 2023

  27. [27]

    Kadosawa, M

    M. Kadosawa, M. Nakamura, Y . Ohta, and S. Nishimoto, Comparing quantum fluctuations in the spin- 1 2 and spin-1 xxz heisenberg models on the square and honeycomb lattices, Phys. Rev. B110, 134418 (2024)

    work page 2024

  28. [28]

    G. H. Wannier, Antiferromagnetism. the triangular ising net, Physical Review79, 357 (1950)

    work page 1950

  29. [29]

    Villain, J., Bidaux, R., Carton, J.-P., and Conte, R., Order as an effect of disorder, J. Phys. France41, 1263 (1980). 11

    work page 1980

  30. [30]

    Shender, Antiferromagnetic garnets with fluctuationally in- teracting sublattices, Soviet Journal of Experimental and Theo- retical Physics56, 178 (1982)

    E. Shender, Antiferromagnetic garnets with fluctuationally in- teracting sublattices, Soviet Journal of Experimental and Theo- retical Physics56, 178 (1982)

    work page 1982

  31. [31]

    A. Sen, P. Dutt, K. Damle, and R. Moessner, Variational wave- function study of the triangular lattice supersolid, Phys. Rev. Lett.100, 147204 (2008)

    work page 2008

  32. [32]

    Nakamura, S

    M. Nakamura, S. Masuda, and S. Nishimoto, Characterization of topological insulators based on the electronic polarization with spiral boundary conditions, Phys. Rev. B104, L121114 (2021)

    work page 2021

  33. [33]

    Kadosawa, M

    M. Kadosawa, M. Nakamura, Y . Ohta, and S. Nishimoto, One- dimensional projection of two-dimensional systems using spiral boundary conditions, Phys. Rev. B107, L081104 (2023)

    work page 2023

  34. [34]

    S. R. White, Density matrix formulation for quantum renormal- ization groups, Phys. Rev. Lett.69, 2863 (1992)

    work page 1992

  35. [35]

    Xiang, J

    T. Xiang, J. Lou, and Z. Su, Two-dimensional algorithm of the density-matrix renormalization group, Phys. Rev. B64, 104414 (2001)

    work page 2001

  36. [36]

    Iqbal, W.-J

    Y . Iqbal, W.-J. Hu, R. Thomale, D. Poilblanc, and F. Becca, Spin liquid nature in the heisenbergJ 1 −J 2 triangular antifer- romagnet, Phys. Rev. B93, 144411 (2016)

    work page 2016

  37. [37]

    Ulaga, J

    M. Ulaga, J. Kokalj, T. Tohyama, and P. Prelov ˇsek, Easy-axis heisenberg model on the triangular lattice: From a supersolid to a gapped solid, Phys. Rev. B111, 174442 (2025)

    work page 2025

  38. [38]

    A. V . Chubukov and D. I. Golosov, Quantum theory of an anti- ferromagnet on a triangular lattice in a magnetic field, Journal of Physics: Condensed Matter3, 69 (1991)

    work page 1991

  39. [39]

    Yamamoto, G

    D. Yamamoto, G. Marmorini, and I. Danshita, Quantum phase diagram of the triangular-latticexxzmodel in a magnetic field, Phys. Rev. Lett.112, 127203 (2014)

    work page 2014

  40. [40]

    Miyashita, Magnetic properties of ising-like heisenberg an- tiferromagnets on the triangular lattice, Journal of the Physical Society of Japan55, 3605 (1986)

    S. Miyashita, Magnetic properties of ising-like heisenberg an- tiferromagnets on the triangular lattice, Journal of the Physical Society of Japan55, 3605 (1986)

    work page 1986

  41. [41]

    Kleine, E

    B. Kleine, E. M ¨uller-Hartmann, K. Frahm, and P. Fazekas, Spin-wave analysis of easy-axis quantum antiferromagnets on the triangular lattice, Zeitschrift f¨ur Physik B Condensed Matter 87, 103 (1992)

    work page 1992

  42. [42]

    Sheng and C

    Q. Sheng and C. L. Henley, Ordering due to disorder in a tri- angular heisenberg antiferromagnet with exchange anisotropy, Journal of Physics: Condensed Matter4, 2937 (1992)

    work page 1992