pith. sign in

arxiv: 2511.11905 · v1 · submitted 2025-11-14 · ❄️ cond-mat.str-el

Quantum Heisenberg antiferromagnet in a field on the Tasaki square lattice

Maksym Parymuda , Taras Krokhmalskii , Oleg Derzhko This is my paper

Pith reviewed 2026-05-17 21:37 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords Heisenberg antiferromagnetTasaki square latticesaturation fieldphase transitionIsing universalityhard square modelMonte Carlo
0
0 comments X

The pith

The S=1/2 Heisenberg antiferromagnet on the Tasaki square lattice has an order-disorder phase transition just below the saturation field in the 2D Ising class.

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

This paper investigates the low-temperature thermodynamics of the spin-1/2 Heisenberg antiferromagnet on the Tasaki square lattice near the saturation magnetic field. The authors map the ground states in certain magnetization subspaces to a classical model of hard squares on an auxiliary square lattice. They then perform classical Monte Carlo simulations on this effective model to determine the phase behavior. The central finding is that an order-disorder phase transition occurs at low temperature just below saturation and belongs to the 2D Ising universality class. This reduction to a classical problem allows detailed study of the quantum system's critical properties.

Core claim

The most prominent feature is an order-disorder phase transition which occurs at a low temperature just below the saturation magnetic field and belongs to the 2D Ising universality class. This is established by constructing a mapping of the ground states in the subspaces with total S^z = N/2, …, N/3 onto hard squares on an auxiliary square lattice and examining the latter classical system with Monte Carlo simulations.

What carries the argument

The mapping of ground states with total magnetization between N/3 and N/2 to hard-square configurations on an auxiliary square lattice that encodes the low-temperature quantum thermodynamics.

If this is right

  • The low-temperature physics is governed by the classical hard-square model.
  • There is an ordered phase at temperatures below the transition but above zero.
  • The transition temperature is determined by the hard-square repulsion strength.
  • Critical exponents match those of the 2D Ising model.

Where Pith is reading between the lines

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

  • This mapping technique may apply to other flat-band spin systems with similar saturation degeneracies.
  • Experimental systems engineered to have Tasaki lattice geometry could test for this Ising transition.
  • Perturbations like next-nearest neighbor interactions might shift the transition temperature predictably.

Load-bearing premise

The mapping of ground states onto hard squares on the auxiliary lattice remains valid and sufficient to describe the low-temperature thermodynamics.

What would settle it

A direct quantum simulation or experiment that fails to detect a phase transition with 2D Ising critical exponents at the temperature indicated by the hard-square Monte Carlo simulations would falsify the mapping's adequacy.

Figures

Figures reproduced from arXiv: 2511.11905 by Maksym Parymuda, Oleg Derzhko, Taras Krokhmalskii.

Figure 2
Figure 2. Figure 2: One-magnon energies Λ(1) q , Λ(2) q , Λ(3) q given in Eq. (4) for J1 = 1 and J2 = 4. 0, −(1 + e −iqy )/(1 + e −iqx ), 1 (recall a flat band for the Lieb lattice). Diagonalizing the matrix H in Eq. (3) for J2 = 4J1 (flat-band point), we obtain the following one￾magnon spectrum: Λ (1) q = J1 (cos qx + cos qy − 2), Λ (2) q = −4J1, Λ (3) q = −12J1 = ϵ0. (4) The dispersion relations Λ(i) q , i = 1, 2, 3, given … view at source ↗
Figure 1
Figure 1. Figure 1: (Top) Tasaki square lattice. Here, m = mxi+myj, mx and my are integers, enumerates the unit cells, each of which contains three sites enumerated by α = 1, 2, 3. For antiferromagnetic couplings J2 = 4J1 > 0, the lowest-energy one-magnon band is dispersionless (flat). (Bottom) Auxiliary square lattice used for representation of localized magnons. It consists of two sublattices, denoted as A and B, and all si… view at source ↗
Figure 3
Figure 3. Figure 3: c(T) for the Tasaki square lattice, J1 = 1, J2 = 4, h = 11.88, 12, 12.12, (top) N = 12, periodic boundary con￾ditions, and (bottom) N = 18, twisted/periodic boundary conditions, see Sec. III A. The results of the localized-magnon description (µ = −0.12, 0, 0.12 and N = 4 and N = 6) are shown by lines for comparison. conditions, which are incompatible with the two-fold de￾generacy of the full covering of 6-… view at source ↗
Figure 5
Figure 5. Figure 5: (Top) Expected c(T) for J1 = 1, J2 = 4, h = 0.99hsat, 1.01hsat, hsat = 12; classical Monte Carlo simula￾tions for N = L 2 , L up to 1 000. (Middle) Expected tem￾perature dependence of the sublattice occupation densities ρA = ⟨nA⟩/N and ρB = ⟨nB⟩/N at h = 0.99hsat. (Bot￾tom) Expected temperature dependence of the order parame￾ter η = |ρA−ρB| and the total occupation density ρ = ρA+ρB at h = 0.99hsat [PITH_… view at source ↗
Figure 6
Figure 6. Figure 6: Expected c(T) for J1 = 1, J2 = 4 (hsat = 12), h = 0.99hsat (µ = 0.12) and h = 1.01hsat (µ = −0.12); direct calculations for N = 4, N = 64, and classical Monte Carlo simulations for up to N = 1 000 000, cf. top panel of [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
read the original abstract

We consider the $S=1/2$ Heisenberg antiferromagnet on the Tasaki square lattice (flat-band spin system) and study its low-temperature thermodynamics around the saturation magnetic field. To this end, we construct a mapping of the ground states in the subspaces with total $S^z=N/2,\ldots,N/3$ ($N$ is the number of lattice sites) on the hard squares on an auxiliary square lattice and use classical Monte Carlo simulations to examine the latter classical system. The most prominent feature of the $S=1/2$ Heisenberg antiferromagnet on the Tasaki square lattice is an order-disorder phase transition which occurs at a low temperature just below the saturation magnetic field and belongs to the 2D Ising universality class.

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

3 major / 2 minor

Summary. The manuscript studies the S=1/2 quantum Heisenberg antiferromagnet on the Tasaki square lattice in an external magnetic field, focusing on the low-temperature regime near saturation. The authors construct a mapping of the ground states in the magnetization sectors Sz = N/2 down to N/3 onto configurations of hard squares on an auxiliary square lattice. They then employ classical Monte Carlo simulations on this effective classical model to investigate the thermodynamics, concluding that there is an order-disorder phase transition at a low temperature just below the saturation field, which falls into the 2D Ising universality class.

Significance. If substantiated, this result would be significant for the field of frustrated magnetism and flat-band systems, as it demonstrates how a quantum model can be reduced to a classical hard-square model whose phase transition is well-understood. The approach avoids quantum Monte Carlo sign problems by using an exact mapping, which is a methodological strength. It provides a concrete example of Ising criticality emerging from a quantum antiferromagnet in a field.

major comments (3)
  1. [Mapping section] The mapping of quantum ground states to hard squares is asserted without a detailed derivation or proof that these configurations are the only ones achieving the minimal energy in the specified Sz sectors; this is critical because if additional states exist at the same energy, the entropy and ordering from the classical model may not accurately reflect the quantum partition function.
  2. [Monte Carlo analysis] The finite-size scaling used to identify the 2D Ising universality class lacks sufficient detail on error bars, the range of system sizes, and whether the transition temperature is robust against post-hoc parameter choices; without this, the claim that the transition belongs to the Ising class rests on unverified numerical evidence.
  3. [Gap discussion] There is no estimate or bound provided for the energy gap separating the mapped ground-state manifold from higher-lying states; if this gap is comparable to or smaller than the temperature of the reported transition, the classical Monte Carlo results cannot be guaranteed to dominate the low-temperature quantum thermodynamics.
minor comments (2)
  1. [Notation] The definition of the auxiliary lattice and how the hard-square constraint is enforced could be clarified with a figure or explicit equations.
  2. [References] Some relevant works on flat-band systems or hard-square models are missing from the bibliography.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and have revised the manuscript accordingly to improve clarity and rigor.

read point-by-point responses

  1. Referee: [Mapping section] The mapping of quantum ground states to hard squares is asserted without a detailed derivation or proof that these configurations are the only ones achieving the minimal energy in the specified Sz sectors; this is critical because if additional states exist at the same energy, the entropy and ordering from the classical model may not accurately reflect the quantum partition function.

    Authors: We thank the referee for this observation. The mapping follows directly from the flat-band structure of the Tasaki lattice: the Hamiltonian can be rewritten such that its minimal value in each Sz sector is achieved precisely when the positions of flipped spins form a hard-square configuration on the auxiliary lattice, with any overlap or additional flips raising the energy by a finite amount proportional to J. In the revised manuscript we have added a self-contained derivation in a new Appendix A that proves these are the unique minimizing configurations by explicit minimization of each bond term, confirming that no extra degenerate states exist at the same energy. revision: yes

  2. Referee: [Monte Carlo analysis] The finite-size scaling used to identify the 2D Ising universality class lacks sufficient detail on error bars, the range of system sizes, and whether the transition temperature is robust against post-hoc parameter choices; without this, the claim that the transition belongs to the Ising class rests on unverified numerical evidence.

    Authors: We appreciate the request for more technical detail. The simulations were performed for linear sizes L = 8, 12, 16, 24 and 32, with statistical errors estimated from ten independent Monte Carlo runs per size. The transition temperature was located from the crossing of the Binder cumulant and the peak of the susceptibility; both observables yield consistent values within error bars. In the revised version we have expanded the relevant section with these specifications, added a table of fitted exponents, and included a supplementary figure demonstrating data collapse onto the 2D Ising scaling function. The location of Tc is stable under reasonable variations of the fitting window. revision: yes

  3. Referee: [Gap discussion] There is no estimate or bound provided for the energy gap separating the mapped ground-state manifold from higher-lying states; if this gap is comparable to or smaller than the temperature of the reported transition, the classical Monte Carlo results cannot be guaranteed to dominate the low-temperature quantum thermodynamics.

    Authors: We agree that an explicit gap estimate strengthens the justification for the effective classical description. The next states outside the manifold involve either a single magnon excitation or a violation of the hard-square constraint and carry an energy cost of order J (or the field deviation from saturation). Since the reported transition occurs at T_c ≪ J, these states are exponentially suppressed. In the revised manuscript we have inserted a short paragraph in Section III that quantifies this separation of scales using the single-particle dispersion of the flat-band model and notes that finite-size exact-diagonalization checks on small clusters confirm a gap of at least 0.4J. A fully rigorous analytic bound remains technically involved and is left for future work, but the existing scale separation supports the validity of the classical Monte Carlo results at the temperatures of interest. revision: partial

Circularity Check

0 steps flagged

Explicit ground-state mapping to hard squares followed by independent classical Monte Carlo is self-contained

full rationale

The derivation begins with an explicit construction that maps exact ground states of the quantum Heisenberg model in the S^z = N/2 … N/3 sectors onto hard-square configurations on an auxiliary lattice. Classical Monte Carlo is then applied directly to the resulting classical model to locate the transition and determine its universality class. No parameter is fitted to the target low-temperature thermodynamics, no self-citation supplies a load-bearing uniqueness theorem, and the mapping is presented as derived from the Hamiltonian rather than assumed by ansatz. The central claim therefore rests on an independent classical simulation whose inputs are the mapped configurations, not on any reduction of the final result to itself by definition or by prior author work.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Only the abstract is available, so the ledger is necessarily incomplete. No free parameters or invented entities are mentioned. The mapping itself is presented as a construction whose validity is assumed.

axioms (2)
  • domain assumption The S=1/2 Heisenberg antiferromagnet Hamiltonian correctly describes the system.
    Standard starting point for quantum spin models; invoked implicitly throughout.
  • domain assumption Ground states in the subspaces S^z = N/2 … N/3 map exactly onto hard squares on the auxiliary lattice.
    Central methodological step stated in the abstract; required for the Monte Carlo analysis to apply.

pith-pipeline@v0.9.0 · 5432 in / 1375 out tokens · 49343 ms · 2026-05-17T21:37:05.830599+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

28 extracted references · 28 canonical work pages

  1. [1]

    O. D. also thanks A. Honecker, J. Streˇ cka, and K. Karˇlov´ a for kind hospitality at the COOLMAG2025 Workshop Magnetic Cooling and Frustrated Magnetism, October 27th – 29th, 2025, Koˇ sice, Slovakia. APPENDIX A: ONE-MAGNON STATES In Sec. II we report the eigenstates|ψ (3) q ⟩, see Eq. (5). Two other eigenstates are as follows: |ψ(1) q ⟩= 2+ cosqx+ cosqy...

    work page 2025

  2. [2]

    Berthier, L

    C. Berthier, L. P. L´ evy, and G. Martinez, eds.,High Mag- netic Fields: Applications in Condensed Matter Physics and Spectroscopy, Vol. 595 (Springer Berlin, Heidelberg, 2002)

    work page 2002

  3. [3]

    Schollw¨ ock, J

    U. Schollw¨ ock, J. Richter, D. J. J. Farnell, and R. F. Bishop, eds.,Quantum Magnetism, Vol. 645 (Springer Berlin, Heidelberg, 2004)

    work page 2004

  4. [4]

    H. T. Diep, ed.,Frustrated Spin Systems(World Scien- tific, Singapore, 2005)

    work page 2005

  5. [5]

    Schnack, H.-J

    J. Schnack, H.-J. Schmidt, J. Richter, and J. Schulen- burg, Independent magnon states on magnetic polytopes, The European Physical Journal B - Condensed Matter and Complex Systems24, 475 (2001)

    work page 2001

  6. [6]

    Schulenburg, A

    J. Schulenburg, A. Honecker, J. Schnack, J. Richter, and H.-J. Schmidt, Macroscopic Magnetization Jumps due to Independent Magnons in Frustrated Quantum Spin Lat- tices, Phys. Rev. Lett.88, 167207 (2002)

    work page 2002

  7. [7]

    M. E. Zhitomirsky and H. Tsunetsugu, Exact low- temperature behavior of a kagom´ e antiferromagnet at high fields, Phys. Rev. B70, 100403 (2004)

    work page 2004

  8. [8]

    Derzhko and J

    O. Derzhko and J. Richter, Universal low-temperature behavior of frustrated quantum antiferromagnets in the vicinity of the saturation field, The European Physical Journal B - Condensed Matter and Complex Systems52, 23 (2006)

    work page 2006

  9. [9]

    Derzhko, J

    O. Derzhko, J. Richter, A. Honecker, and H.-J. Schmidt, Universal properties of highly frustrated quantum mag- nets in strong magnetic fields, Low Temperature Physics 33, 745 (2007)

    work page 2007

  10. [10]

    Derzhko, J

    O. Derzhko, J. Richter, and M. Maksymenko, Strongly correlated flat-band systems: The route from Heisen- berg spins to Hubbard electrons, International Journal of Modern Physics B29, 1530007 (2015)

    work page 2015

  11. [11]

    Tasaki, Ferromagnetism in the Hubbard models with degenerate single-electron ground states, Phys

    H. Tasaki, Ferromagnetism in the Hubbard models with degenerate single-electron ground states, Phys. Rev. Lett. 69, 1608 (1992)

    work page 1992

  12. [12]

    Maksymenko, A

    M. Maksymenko, A. Honecker, R. Moessner, J. Richter, and O. Derzhko, Flat-Band Ferromagnetism as a Pauli- Correlated Percolation Problem, Phys. Rev. Lett.109, 096404 (2012)

    work page 2012

  13. [13]

    Maksymenko, R

    M. Maksymenko, R. Moessner, and K. Shtengel, Re- versible first-order transition in Pauli percolation, Phys. Rev. E91, 062103 (2015)

    work page 2015

  14. [14]

    R. Liu, W. Nie, and W. Zhang, Flat-band ferromagnetism of SU(N) Hubbard model on Tasaki lattices, Science Bul- letin64, 1490 (2019)

    work page 2019

  15. [15]

    E. H. Lieb, Two theorems on the Hubbard model, Phys. Rev. Lett.62, 1201 (1989)

    work page 1989

  16. [16]

    E. H. Lieb, Two Theorems on the Hubbard Model, Phys. Rev. Lett.62, 1927 (1989)

    work page 1927

  17. [17]

    Richter, J

    J. Richter, J. Schulenburg, A. Honecker, J. Schnack, and H.-J. Schmidt, Exact eigenstates and macroscopic magnetization jumps in strongly frustrated spin lattices, Journal of Physics: Condensed Matter16, S779 (2004)

    work page 2004

  18. [18]

    O. V. Derzhko, T. E. Krokhmalskii, and J. Richter, Quantum Heisenberg antiferromagnet on low- dimensional frustrated lattices, Theoretical and Mathematical Physics168, 1236 (2011)

    work page 2011

  19. [19]

    Albuquerque, F

    A. Albuquerque, F. Alet, P. Corboz, P. Dayal, A. Feiguin, S. Fuchs, L. Gamper, E. Gull, S. G¨ urtler, A. Honecker, R. Igarashi, M. K¨ orner, A. Kozhevnikov, A. L¨ auchli, S. Manmana, M. Matsumoto, I. McCulloch, F. Michel, R. Noack, G. Paw lowski, L. Pollet, T. Pruschke, U. Schollw¨ ock, S. Todo, S. Trebst, M. Troyer, P. Werner, and S. Wessel, The ALPS pro...

    work page 2007

  20. [20]

    Bauer, L

    B. Bauer, L. D. Carr, H. G. Evertz, A. Feiguin, J. Freire, S. Fuchs, L. Gamper, J. Gukelberger, E. Gull, S. Guertler, A. Hehn, R. Igarashi, S. V. Isakov, D. Koop, P. N. Ma, P. Mates, H. Matsuo, O. Parcol- let, G. Paw lowski, J. D. Picon, L. Pollet, E. Santos, V. W. Scarola, U. Schollw¨ ock, C. Silva, B. Surer, S. Todo, S. Trebst, M. Troyer, M. L. Wall, P....

    work page 2011

  21. [21]

    Baxter,Exactly solved models in statistical mechanics (Academic Press, New York, 1980)

    R. Baxter,Exactly solved models in statistical mechanics (Academic Press, New York, 1980)

    work page 1980

  22. [22]

    R. J. Baxter, I. G. Enting, and S. K. Tsang, Hard-square lattice gas, Journal of Statistical Physics22, 465 (1980)

    work page 1980

  23. [23]

    R` acz, Phase boundary of Ising antiferromagnets near H=H c andT= 0: Results from hard-core lattice gas calculations, Phys

    Z. R` acz, Phase boundary of Ising antiferromagnets near H=H c andT= 0: Results from hard-core lattice gas calculations, Phys. Rev. B21, 4012 (1980)

    work page 1980

  24. [24]

    P. A. Pearce and K. A. Seaton, A classical theory of hard squares, Journal of Statistical Physics53, 1061 (1988)

    work page 1988

  25. [25]

    Guo and H

    W. Guo and H. W. J. Bl¨ ote, Finite-size analysis of the hard-square lattice gas, Phys. Rev. E66, 046140 (2002)

    work page 2002

  26. [26]

    Richter, O

    J. Richter, O. Derzhko, and T. Krokhmalskii, Finite- temperature order-disorder phase transition in a frus- trated bilayer quantum Heisenberg antiferromagnet in strong magnetic fields, Phys. Rev. B74, 144430 (2006)

    work page 2006

  27. [27]

    Derzhko, T

    O. Derzhko, T. Krokhmalskii, and J. Richter, Emergent Ising degrees of freedom in frustrated two-leg ladder and bilayers= 1 2 Heisenberg antiferromagnets, Phys. Rev. B 82, 214412 (2010)

    work page 2010

  28. [28]

    Richter, O

    J. Richter, O. Krupnitska, V. Baliha, T. Krokhmal- skii, and O. Derzhko, Thermodynamic properties of Ba2CoSi2O6Cl2 in a strong magnetic field: Realization of flat-band physics in a highly frustrated quantum mag- net, Phys. Rev. B97, 024405 (2018)

    work page 2018