pith. machine review for the scientific record. sign in

arxiv: 2604.24687 · v1 · submitted 2026-04-27 · ❄️ cond-mat.str-el

Recognition: unknown

Field-induced jammed polyhex spin liquid in the honeycomb Ising antiferromagnet

Nicholas Franklin , Jacob Richards , Harry Lane
Authors on Pith no claims yet

Pith reviewed 2026-05-08 01:56 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords Ising antiferromagnethoneycomb latticepolyhex tilingjammed spin liquidground state degeneracymagnetic fieldquantum fluctuations
0
0 comments X

The pith

A magnetic field induces extensive ground-state degeneracy in the honeycomb Ising antiferromagnet that maps exactly to polyhex tilings and forms a jammed spin liquid with no local zero modes.

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

The authors examine the ground states of the honeycomb lattice Ising antiferromagnet when a magnetic field is applied. They find that over a range of field values the ground states are extensively degenerate and can be mapped one-to-one onto all possible ways to tile the plane with polyhexes. Because no local spin flip connects one such ground state to another, the manifold is jammed rather than liquid in the usual sense. Monte Carlo runs show that the spins can still rearrange, but the diffusion rate drops as the field is increased. Adding a small transverse exchange term lifts part of the degeneracy and selects a subspace of aperiodic 12-hex tilings.

Core claim

At finite longitudinal magnetic field the classical ground-state manifold of the honeycomb Ising antiferromagnet is in exact correspondence with the set of all polyhex tilings of the plane. These states form a jammed spin liquid: the degeneracy is extensive yet there exist no local zero-energy moves that link distinct configurations. Monte Carlo sampling of the manifold demonstrates that spin diffusion remains possible but is suppressed by increasing field strength. A small transverse coupling partially lifts the degeneracy, selecting an extensive set of non-periodic tilings composed of 12-hex units.

What carries the argument

The exact mapping between Ising spin configurations and polyhex tilings of the honeycomb lattice, which guarantees the absence of local zero modes.

If this is right

  • The ground-state degeneracy exists over a finite interval of field strengths.
  • Spin diffusion constant is continuously tunable by the applied field.
  • Quantum fluctuations from transverse couplings select non-periodic 12-hex tilings.
  • The phase may be realized in the FePX3 family of compounds.

Where Pith is reading between the lines

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

  • This jammed character may lead to unusual slow relaxation dynamics at low temperatures not explored in the work.
  • The field-tunable diffusion suggests possible applications in controlling magnetic response in 2D materials.
  • Similar mappings might appear in other frustrated lattices with field-induced degeneracies.
  • Experimental detection could involve measuring the field dependence of the specific heat or susceptibility in candidate materials.

Load-bearing premise

The correspondence between the Ising ground states and the full set of polyhex tilings is exact, with no additional constraints imposed by the lattice or the field that would reduce the degeneracy.

What would settle it

A numerical count of ground states on a large finite honeycomb cluster that differs from the known number of polyhex tilings for the same area would falsify the mapping.

Figures

Figures reproduced from arXiv: 2604.24687 by Harry Lane, Jacob Richards, Nicholas Franklin.

Figure 1
Figure 1. Figure 1: FIG. 1. Phase diagram of the Ising Honeycomb model, based view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Magnetization as a function of temperature show view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. (a) view at source ↗
read the original abstract

We analyze the ground state properties of the honeycomb Ising antiferromagnet in an external magnetic field. We demonstrate the existence of extensive ground state degeneracy at finite field that maps to a polyhex tiling problem. This state is shown to be a jammed spin liquid, with no local zero modes connecting ground states. Through Monte Carlo simulations, we explore the properties of these states and show that the spin diffusion can be controlled by the magnetic field strength. By considering quantum fluctuations, we demonstrate that transverse coupling partially lifts the ground state degeneracy, selecting an extensive subspace of non-periodic tilings of 12-hexes. We suggest that the jammed polyhex spin liquid phase exists in an experimentally realizable region of parameter space and may be present in the FePX$_3$ compounds.

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 / 3 minor

Summary. The manuscript analyzes the ground state properties of the honeycomb Ising antiferromagnet in an external magnetic field. It demonstrates extensive ground state degeneracy at finite field that maps to a polyhex tiling problem. This state is shown to be a jammed spin liquid with no local zero modes connecting ground states. Monte Carlo simulations explore the properties, showing field-controlled spin diffusion. Quantum fluctuations via transverse coupling partially lift the degeneracy, selecting an extensive subspace of non-periodic 12-hex tilings. The phase is suggested to exist in an experimentally realizable region and may be present in FePX3 compounds.

Significance. If the exact mapping to polyhex tilings and the confirmation of no local zero modes hold, this provides a concrete example of a field-induced jammed classical spin liquid with combinatorial degeneracy. The link to tiling problems and the Monte Carlo demonstration of field-tunable diffusion offer testable predictions. The quantum perturbation analysis adds a mechanism for partial degeneracy lifting. The work is internally consistent on the classical side with no evident overcounting in the energy minimization.

major comments (2)
  1. [§3] §3 (mapping derivation): The one-to-one correspondence between field-tuned Ising ground states and polyhex tilings is asserted via energy minimization, but the explicit enumeration or bijection proof (including handling of boundary conditions on finite lattices) is not detailed enough to confirm no overcounting or missed constraints; this is load-bearing for the extensive degeneracy claim.
  2. [§4] §4 (Monte Carlo): The validation of absent local zero modes and field-dependent diffusion lacks reported error bars, system-size scaling, and the precise move set (e.g., single flips vs. loop updates) used to probe connectivity; without these, the 'jammed' characterization and diffusion control remain qualitative.
minor comments (3)
  1. [Abstract] Abstract and introduction: The term 'jammed polyhex spin liquid' is used before any definition or contrast to conventional spin liquids; a one-sentence clarification would improve accessibility.
  2. [Figures] Figure captions (e.g., diffusion plots): Axes lack explicit units or labels for field strength values; legends should specify the Monte Carlo parameters (temperature, sweeps) for reproducibility.
  3. [References] References: Prior works on honeycomb Ising antiferromagnets in field and on polyhex enumeration algorithms are not cited; adding 2-3 key references would contextualize the novelty.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the positive assessment of its internal consistency and potential significance. We address each major comment below.

read point-by-point responses

  1. Referee: [§3] §3 (mapping derivation): The one-to-one correspondence between field-tuned Ising ground states and polyhex tilings is asserted via energy minimization, but the explicit enumeration or bijection proof (including handling of boundary conditions on finite lattices) is not detailed enough to confirm no overcounting or missed constraints; this is load-bearing for the extensive degeneracy claim.

    Authors: We agree that a more explicit demonstration of the bijection is warranted given its central role in establishing extensive degeneracy. The manuscript derives the mapping by minimizing the field-dependent energy, which restricts each hexagon to configurations with a fixed number of down spins that exactly match the edge placements in a polyhex tiling. To strengthen this, the revised manuscript will add an appendix containing (i) explicit enumeration of all ground states on small periodic clusters (e.g., 2×2 and 3×3 supercells) showing one-to-one correspondence with valid polyhex tilings, and (ii) a concise argument that, under periodic boundary conditions in the thermodynamic limit, every valid tiling corresponds to a unique ground state with no overcounting or additional constraints. These additions will be made without changing the reported results. revision: yes

  2. Referee: [§4] §4 (Monte Carlo): The validation of absent local zero modes and field-dependent diffusion lacks reported error bars, system-size scaling, and the precise move set (e.g., single flips vs. loop updates) used to probe connectivity; without these, the 'jammed' characterization and diffusion control remain qualitative.

    Authors: We concur that additional technical details are needed to make the Monte Carlo evidence fully quantitative. The simulations employ single-spin-flip Metropolis dynamics to establish the absence of local zero modes (distinct ground states remain disconnected under single flips) and to measure field-dependent diffusion via the mean-squared displacement of spin flips. In the revision we will (i) report statistical error bars obtained from independent runs, (ii) present finite-size scaling of the connectivity measure across lattices from 20×20 to 100×100 to confirm that the jammed property persists, and (iii) explicitly state the move set, including the use of supplementary loop updates solely for equilibration diagnostics. These changes will render the field-controlled diffusion quantitative while leaving the qualitative conclusions unchanged. revision: yes

Circularity Check

0 steps flagged

Minor self-citation present but derivation chain remains independent

full rationale

The paper derives the field-induced degeneracy by direct minimization of the classical Ising Hamiltonian on the honeycomb lattice, mapping configurations to polyhex tilings via explicit energy rules without redefining inputs in terms of outputs. Monte Carlo sampling of diffusion and zero-mode absence is performed on the resulting ensemble rather than fitted to it. Quantum lifting via transverse terms is treated as a separate perturbative analysis. Any self-citations are peripheral and do not carry the central mapping or degeneracy count, satisfying the criteria for a low circularity score.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The claim rests on the standard Ising antiferromagnet Hamiltonian plus the assumption that finite-field ground states map exactly onto polyhex tilings; the jammed property and quantum lifting are additional postulates without independent evidence supplied in the abstract.

axioms (2)
  • domain assumption Spins are classical Ising variables on the honeycomb lattice with nearest-neighbor antiferromagnetic coupling.
    Standard starting point for the model; invoked implicitly throughout the abstract.
  • ad hoc to paper The external field induces extensive degeneracy that is in one-to-one correspondence with polyhex tilings.
    Central mapping asserted without derivation details in the provided abstract.
invented entities (1)
  • jammed polyhex spin liquid no independent evidence
    purpose: Name for the field-induced degenerate ground-state manifold with no local zero modes.
    New descriptive term introduced to characterize the phase.

pith-pipeline@v0.9.0 · 5423 in / 1269 out tokens · 39635 ms · 2026-05-08T01:56:33.676996+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

49 extracted references · 2 canonical work pages

  1. [1]

    H. Yan, O. Benton, A. H. Nevidomskyy, and R. Moess- ner, Classification of classical spin liquids: Detailed for- malism and suite of examples, Phys. Rev. B109, 174421 (2024)

    2024

  2. [2]

    H. Yan, O. Benton, R. Moessner, and A. H. Nevidom- skyy, Classification of classical spin liquids: Typology and resulting landscape, Phys. Rev. B110, L020402 (2024)

    2024

  3. [3]

    Lozano-G´ omez, O

    D. Lozano-G´ omez, O. Benton, M. J. P. Gingras, and H. Yan, An atlas of classical pyrochlore spin liquids (2025), arXiv:2411.03547 [cond-mat.str-el]

    work page arXiv 2025

  4. [4]

    Mulder, R

    A. Mulder, R. Ganesh, L. Capriotti, and A. Paramekanti, Spiral order by disorder and lattice nematic order in a frustrated heisenberg antiferromagnet on the honeycomb lattice, Phys. Rev. B81, 214419 (2010)

    2010

  5. [5]

    Reuther, R

    J. Reuther, R. Thomale, and S. Rachel, Spiral order in the honeycomb iridate Li2IrO3, Phys. Rev. B90, 100405 (2014)

    2014

  6. [6]

    Niggemann, M

    N. Niggemann, M. Hering, and J. Reuther, Classical spi- ral spin liquids as a possible route to quantum spin liq- uids, J. Phys. Condens. Matter32, 024001 (2019)

    2019

  7. [7]

    Yan and J

    H. Yan and J. Reuther, Low-energy structure of spiral spin liquids, Phys. Rev. Res.4, 023175 (2022)

    2022

  8. [8]

    and Bak, P., Two-dimensional ising model with competing interactions : floating phase, walls and dislo- cations, J

    Villain, J. and Bak, P., Two-dimensional ising model with competing interactions : floating phase, walls and dislo- cations, J. Phys. France42, 657 (1981)

    1981

  9. [9]

    Moessner, S

    R. Moessner, S. L. Sondhi, and P. Chandra, Two- dimensional periodic frustrated ising models in a trans- verse field, Phys. Rev. Lett.84, 4457 (2000)

    2000

  10. [10]

    Selke, The ANNNI model — theoretical analysis and experimental application, Phys

    W. Selke, The ANNNI model — theoretical analysis and experimental application, Phys. Rep.170, 213 (1988)

    1988

  11. [11]

    M. E. Fisher and W. Selke, Infinitely Many Commensu- rate Phases in a Simple Ising Model, Phys. Rev. Lett.44, 1502 (1980)

    1980

  12. [12]

    G. H. Wannier, Antiferromagnetism. the triangular ising net, Phys. Rev.79, 357 (1950)

    1950

  13. [13]

    S. Jin, A. Sen, and A. W. Sandvik, Ashkin-teller critical- ity and pseudo-first-order behavior in a frustrated ising model on the square lattice, Phys. Rev. Lett.108, 045702 (2012)

    2012

  14. [14]

    S. T. Bramwell and M. J. P. Gingras, Spin ice state in frustrated magnetic pyrochlore materials, Science294, 1495 (2001)

    2001

  15. [15]

    Binder and A

    K. Binder and A. P. Young, Spin glasses: Experimen- tal facts, theoretical concepts, and open questions, Rev. Mod. Phys.58, 801 (1986)

    1986

  16. [16]

    Klich, S.-H

    I. Klich, S.-H. Lee, and K. Iida, Glassiness and exotic entropy scaling induced by quantum fluctuations in a disorder-free frustrated magnet, Nat. Commun.5, 3497 (2014)

    2014

  17. [17]

    Bilitewski, M

    T. Bilitewski, M. E. Zhitomirsky, and R. Moessner, Jammed spin liquid in the bond-disordered kagome anti- ferromagnet, Phys. Rev. Lett.119, 247201 (2017)

    2017

  18. [18]

    Bilitewski, M

    T. Bilitewski, M. E. Zhitomirsky, and R. Moessner, Dy- namics and energy landscape of the jammed spin liquid, Phys. Rev. B99, 054416 (2019)

    2019

  19. [19]

    Parisi and F

    G. Parisi and F. Zamponi, Mean-field theory of hard sphere glasses and jamming, Rev. Mod. Phys.82, 789 (2010). 6

    2010

  20. [20]

    Yan and R

    H. Yan and R. Pohle, ClassicalZ 2 spin liquid on the generalized four-color kitaev model, Phys. Rev. Res.7, L012052 (2025)

    2025

  21. [21]

    Kondev and C

    J. Kondev and C. L. Henley, Four-coloring model on the square lattice: A critical ground state, Phys. Rev. B52, 6628 (1995)

    1995

  22. [22]

    Lozano-G´ omez, Y

    D. Lozano-G´ omez, Y. Iqbal, and M. Vojta, A classi- cal chiral spin liquid from chiral interactions on the py- rochlore lattice, Nat. Commun.15, 10162 (2024)

    2024

  23. [23]

    D. S. Rokhsar and S. A. Kivelson, Superconductivity and the quantum hard-core dimer gas, Phys. Rev. Lett.61, 2376 (1988)

    1988

  24. [24]

    Moessner and S

    R. Moessner and S. L. Sondhi, Ising models of quantum frustration, Phys. Rev. B63, 224401 (2001)

    2001

  25. [25]

    Sachdev, Spin-peierls ground states of the quantum dimer model: A finite-size study, Phys

    S. Sachdev, Spin-peierls ground states of the quantum dimer model: A finite-size study, Phys. Rev. B40, 5204 (1989)

    1989

  26. [26]

    Moessner, S

    R. Moessner, S. L. Sondhi, and P. Chandra, Phase di- agram of the hexagonal lattice quantum dimer model, Phys. Rev. B64, 144416 (2001)

    2001

  27. [27]

    P. W. Leung, K. C. Chiu, and K. J. Runge, Colum- nar dimer and plaquette resonating-valence-bond orders in the quantum dimer model, Phys. Rev. B54, 12938 (1996)

    1996

  28. [28]

    D. A. Ivanov, Vortexlike elementary excitations in the rokhsar-kivelson dimer model on the triangular lattice, Phys. Rev. B70, 094430 (2004)

    2004

  29. [29]

    Senthil and M

    T. Senthil and M. P. A. Fisher,Z 2 gauge theory of electron fractionalization in strongly correlated systems, Phys. Rev. B62, 7850 (2000)

    2000

  30. [30]

    Gardner, Mathematical games, Sci

    M. Gardner, Mathematical games, Sci. Am.216, 124 (1967)

    1967

  31. [31]

    J. Yang, A. Samarakoon, S. Dissanayake, H. Ueda, I. Klich, K. Iida, D. Pajerowski, N. P. Butch, Q. Huang, J. R. D. Copley, and S.-H. Lee, Spin jam induced by quantum fluctuations in a frustrated magnet, Proc. Natl. Acad. Sci. USA112, 11519 (2015)

    2015

  32. [32]

    Samarakoon, T

    A. Samarakoon, T. J. Sato, T. Chen, G.-W. Chern, J. Yang, I. Klich, R. Sinclair, H. Zhou, and S.-H. Lee, Aging, memory, and nonhierarchical energy landscape of spin jam, Proc. Natl. Acad. Sci. USA113, 11806 (2016)

    2016

  33. [33]

    Piyakulworawat, A

    C. Piyakulworawat, A. Thennakoon, J. Yang, H. Yoshizawa, D. Ueta, T. J. Sato, K. Sheng, W.- T. Chen, W.-W. Pai, K. Matan, and S.-H. Lee, Zero-point entropies of spin-jam and spin-glass states in a frustrated magnet, Phys. Rev. B109, 104420 (2024)

    2024

  34. [34]

    A. R. Wildes, K. C. Rule, R. I. Bewley, M. Enderle, and T. J. Hicks, The magnon dynamics and spin exchange pa- rameters of FePS3, J. Phys. Condens. Matter24, 416004 (2012)

    2012

  35. [35]

    A. R. Wildes, D. Lan¸ con, M. K. Chan, F. Weickert, N. Harrison, V. Simonet, M. E. Zhitomirsky, M. V. Gvozdikova, T. Ziman, and H. M. Rønnow, High field magnetization of FePS 3, Phys. Rev. B101, 024415 (2020)

    2020

  36. [36]

    A. R. Wildes, B. F˚ ak, U. B. Hansen, M. Enderle, J. R. Stewart, L. Testa, H. M. Rønnow, C. Kim, and J.-G. Park, Spin wave spectra of single crystal CoPS 3, Phys. Rev. B107, 054438 (2023)

    2023

  37. [37]

    L. Chen, X. Teng, D. Hu, F. Ye, G. E. Granroth, M. Yi, J.-H. Chung, R. J. Birgeneau, and P. Dai, Thermal evo- lution of spin excitations in honeycomb ising antiferro- magnetic FePSe3, Npj Quantum Mater.9, 40 (2024)

    2024

  38. [38]

    Le Mardel´ e, A

    F. Le Mardel´ e, A. El Mendili, M. E. Zhitomirsky, I. Mo- helsky, D. Jana, I. Plutnarova, Z. Sofer, C. Faugeras, M. Potemski, and M. Orlita, Transverse and longitudi- nal magnons in the strongly anisotropic antiferromagnet FePSe3, Phys. Rev. B109, 134410 (2024)

    2024

  39. [39]

    Kud¯ o and S

    T. Kud¯ o and S. Katsura, A method of determining the orderings of the ising model with several neighbor in- teractions under the magnetic field and applications to hexagonal lattices, Prog. Theor. Phys.56, 435 (1976)

    1976

  40. [40]

    S. W. Golomb and L. D. Baumert, Backtrack program- ming, J. ACM12, 516–524 (1965)

    1965

  41. [41]

    Scott,Programming a Combinatorial Puzzle, Tech

    D. Scott,Programming a Combinatorial Puzzle, Tech. Rep. 1 (Princeton University, 1958)

    1958

  42. [42]

    D. E. Knuth, Dancing links (2000), arXiv:cs/0011047 [cs.DS]

    work page Pith review arXiv 2000

  43. [43]

    This leads to an under-count since there is no require- ment that notches and protrusions on opposite side of the supercell tessellate in the limit that the supercell vol- ume tends to infinity

  44. [44]

    Inclusion of four symmetry elements reducesAby a factor of four but does not affect the scaling

    Our procedure treats states that differ by the symmetry elements of the supercell as distinct. Inclusion of four symmetry elements reducesAby a factor of four but does not affect the scaling

  45. [45]

    This scaling fails in the limit of smallN s where the poly- hex size is comparable with the supercell size

  46. [46]

    R. H. Swendsen and J.-S. Wang, Replica monte carlo sim- ulation of spin-glasses, Phys. Rev. Lett.57, 2607 (1986)

    1986

  47. [47]

    Dahlbom, H

    D. Dahlbom, H. Zhang, C. Miles, S. Quinn, A. Niraula, B. Thipe, M. Wilson, S. Matin, H. Mankad, S. Hahn, D. Pajerowski, S. Johnston, Z. Wang, H. Lane, Y. W. Li, X. Bai, M. Mourigal, C. D. Batista, and K. Barros, Sunny.jl: A julia package for spin dynamics, J. Open Source Softw.10, 8138 (2025)

    2025

  48. [48]

    G. C. Rhoads, Planar tilings by polyominoes, polyhexes, and polyiamonds, J. Comput. Appl. Math.174, 329 (2005)

    2005

  49. [49]

    Hermele, M

    M. Hermele, M. P. A. Fisher, and L. Balents, Py- rochlore photons: Theu(1) spin liquid in as= 1 2 three- dimensional frustrated magnet, Phys. Rev. B69, 064404 (2004)

    2004