pith. sign in

arxiv: 2510.02042 · v2 · submitted 2025-10-02 · ❄️ cond-mat.str-el

The Finite-Temperature Behavior of a Triangular Heisenberg Antiferromagnet

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

classification ❄️ cond-mat.str-el
keywords triangular latticeHeisenberg antiferromagnetspiral spin liquidfinite temperaturephase diagramspecific heat humpNematic Bond Theorystructure factor
0
0 comments X

The pith

Along a specific line of couplings the triangular Heisenberg model develops a ring of degenerate minima that produces spiral spin liquid behavior at intermediate temperatures.

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

The paper uses Nematic Bond Theory to compute the free energy and static structure factor of the classical Heisenberg antiferromagnet on the triangular lattice with interactions up to third neighbors. It constructs the full finite-temperature phase diagram and isolates the special line where the Fourier transform of the exchange has a continuous ring of minima. On that line the system forms a spiral spin liquid whose structure factor evolves with the second-neighbor coupling strength; the low-temperature ordered state that emerges is the single-q spiral carrying the largest spin-wave entropy on the ring. The same calculations identify the broad specific-heat hump inside the 120-degree phase as the temperature at which the correlation length begins to grow exponentially.

Core claim

Along the line J3 = J2/2 the Fourier-transformed exchange exhibits a degenerate ring-like minimum, giving rise to spiral spin liquid behavior at intermediate temperatures. The low-temperature order coincides with the single-q spiral states of maximum spin-wave entropy along the degenerate ring.

What carries the argument

Nematic Bond Theory, which supplies a variational free-energy functional whose minimization yields both the phase boundaries and the momentum-space structure factor at any temperature.

If this is right

  • The broad specific-heat hump inside the 120-degree phase marks the onset of an exponentially growing correlation length.
  • Symmetry-breaking transition temperatures can be tracked continuously across the entire J2-J3 plane.
  • The structure factor of the spiral spin liquid changes systematically with J2 while the ring minimum is preserved.

Where Pith is reading between the lines

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

  • The entropy-maximization rule for selecting the low-temperature spiral may apply to other classical frustrated magnets whose ground-state manifold is a continuous manifold in q-space.
  • Neutron-scattering experiments tuned to the J3 = J2/2 line could directly image the ring-shaped diffuse scattering that defines the spiral spin liquid regime.
  • The same Nematic Bond Theory functional could be applied to the quantum S = 1/2 version to test whether the spiral spin liquid survives zero-point fluctuations.

Load-bearing premise

Minimizing the Nematic Bond Theory free-energy functional produces the correct finite-temperature phases and structure factors for the classical model.

What would settle it

A Monte Carlo simulation or exact enumeration that finds the low-temperature ordered state on the J3 = J2/2 line to be a different single-q spiral than the one with maximum spin-wave entropy would falsify the selection rule.

Figures

Figures reproduced from arXiv: 2510.02042 by Cecilie Glittum, Olav F. Sylju{\aa}sen.

Figure 1
Figure 1. Figure 1: FIG. 1. Triangular lattice with first- ( [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Zero-temperature phase diagram with [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Contour plot of the static structure factor for (a) [PITH_FULL_IMAGE:figures/full_fig_p005_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Phase transition temperature [PITH_FULL_IMAGE:figures/full_fig_p006_8.png] view at source ↗
read the original abstract

We investigate the classical antiferromagnetic Heisenberg model on the triangular lattice with up to third-nearest neighbor exchange couplings using the Nematic Bond Theory. This approach allows us to compute the free energy and the neutron scattering static structure factor at finite temperatures. We map out the phase diagram with a particular emphasis on finite-temperature phase transitions that break lattice-rotational symmetries, spiral spin liquids and the broad specific heat hump that is ubiquitous in the antiferromagnetic 120 degree phase. We identify this specific heat hump as signaling the onset of an exponentially increasing correlation length. Further, we map out the temperature of the specific heat hump and the transition temperatures of the symmetry-breaking transitions throughout the exchange-coupling space. Along the line $J_3 = J_2/2$, the Fourier-transformed exchange coupling exhibits a degenerate ring-like minimum, giving rise to spiral spin liquid behavior at intermediate temperatures. We investigate the structure factor of the spiral spin liquid as function of $J_2$ and identify the corresponding low-temperature order, which coincides with the single-$\vec{q}$ spiral states of maximum spin-wave entropy along the degenerate ring.

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 paper investigates the classical antiferromagnetic Heisenberg model on the triangular lattice with up to third-nearest neighbor exchange couplings using the Nematic Bond Theory. This approach allows computation of the free energy and neutron scattering static structure factor at finite temperatures. The authors map the phase diagram emphasizing lattice-rotational symmetry-breaking transitions, spiral spin liquids, and the specific-heat hump in the 120-degree phase, which they identify as the onset of exponentially growing correlation lengths. Along the J3 = J2/2 line, a degenerate ring-like minimum in the Fourier-transformed exchange gives rise to intermediate-temperature spiral spin liquid behavior, with low-temperature order identified as the single-q spirals of maximum spin-wave entropy.

Significance. If the Nematic Bond Theory free-energy minimization is quantitatively reliable, the manuscript delivers a detailed finite-temperature phase diagram for a classically frustrated model, including predictions for structure factors and specific-heat features relevant to triangular-lattice materials. The explicit connection between the ring degeneracy, spiral spin liquid regime, and entropy-selected single-q ordering constitutes a useful theoretical contribution to the study of degenerate manifolds in spin systems.

major comments (2)
  1. [Abstract and J3 = J2/2 results] The phase boundaries, structure-factor identification of the spiral spin liquid, and low-T single-q selection along J3 = J2/2 all rest on minimization of the NBT free-energy functional. The manuscript provides no benchmarks against Monte Carlo simulations or exact diagonalization for the classical model (see the abstract and the section describing the J3 = J2/2 results). Without such validation, it remains unclear whether the bond-nematic truncation correctly captures soft-mode entropy contributions along the continuous ring degeneracy, which could alter the predicted ordering wave-vector.
  2. [Specific-heat hump discussion] The claim that the specific-heat hump marks the onset of an exponentially increasing correlation length is obtained directly from the same NBT free-energy minimization. An explicit extraction of the correlation length (e.g., from the width of the structure-factor peaks or an auxiliary calculation) as a function of temperature would be required to substantiate this identification and distinguish it from other crossovers.
minor comments (2)
  1. [Methods] Clarify the precise definition of the nematic bond order parameter and its relation to the spin-wave entropy argument used for low-T selection.
  2. [Phase diagram] Add a brief comparison table or figure overlaying NBT transition temperatures with any available literature values for limiting cases (e.g., pure J1 or J1-J2 models).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We respond to each major comment below and indicate the revisions we will make.

read point-by-point responses
  1. Referee: [Abstract and J3 = J2/2 results] The phase boundaries, structure-factor identification of the spiral spin liquid, and low-T single-q selection along J3 = J2/2 all rest on minimization of the NBT free-energy functional. The manuscript provides no benchmarks against Monte Carlo simulations or exact diagonalization for the classical model (see the abstract and the section describing the J3 = J2/2 results). Without such validation, it remains unclear whether the bond-nematic truncation correctly captures soft-mode entropy contributions along the continuous ring degeneracy, which could alter the predicted ordering wave-vector.

    Authors: We agree that explicit benchmarks against Monte Carlo would strengthen the presentation. The Nematic Bond Theory is a variational approximation whose accuracy has been established in earlier applications to related frustrated models; the ring degeneracy itself is an exact feature of the exchange Hamiltonian, and the low-temperature single-q selection follows from the standard spin-wave entropy calculation, which is non-perturbative in that limit. Nevertheless, to address the concern about intermediate-temperature soft-mode contributions, we will add a brief comparison with existing Monte Carlo data for the J2-only and J3-only limits together with a short discussion of the expected accuracy of the bond-nematic truncation along the degenerate line. revision: partial

  2. Referee: [Specific-heat hump discussion] The claim that the specific-heat hump marks the onset of an exponentially increasing correlation length is obtained directly from the same NBT free-energy minimization. An explicit extraction of the correlation length (e.g., from the width of the structure-factor peaks or an auxiliary calculation) as a function of temperature would be required to substantiate this identification and distinguish it from other crossovers.

    Authors: The referee is correct that the identification currently rests on the temperature dependence of the minimized free energy and the associated structure factor. To make the claim more direct, we will add an explicit extraction of the correlation length from the inverse width of the structure-factor peaks (or an equivalent auxiliary calculation) as a function of temperature in the 120-degree phase, thereby confirming the exponential growth and clarifying the nature of the crossover. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation applies NBT as an independent approximation tool

full rationale

The paper applies Nematic Bond Theory to minimize an approximate free-energy functional and extract the structure factor, phase boundaries, specific-heat features, and spiral spin liquid signatures along the J3 = J2/2 line. These quantities are direct computational outputs of the functional rather than inputs renamed as predictions. The low-temperature single-q selection is presented as coinciding with an independent spin-wave entropy calculation along the degeneracy ring, which is not derived from or forced by the finite-T NBT minimization. No self-definitional loops, fitted parameters called predictions, or load-bearing self-citations that reduce the central claims to unverified inputs appear in the derivation chain. The method is used with explicit approximations whose validity is external to the present results.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The calculation assumes the classical limit (S→∞), the validity of the Nematic Bond Theory approximation for the free energy, and that the Fourier transform of the exchange couplings fully determines the degenerate minima. No new particles or forces are introduced.

free parameters (1)
  • J2/J1 and J3/J1
    The two independent ratios of further-neighbor to nearest-neighbor exchange are scanned to produce the phase diagram; they are external parameters of the model rather than fitted quantities.
axioms (1)
  • domain assumption Nematic Bond Theory yields a sufficiently accurate variational free energy for the classical Heisenberg model on the triangular lattice.
    The entire finite-temperature analysis rests on this approximation being reliable enough to locate phase boundaries and identify the specific-heat feature.

pith-pipeline@v0.9.0 · 5730 in / 1416 out tokens · 21594 ms · 2026-05-18T10:29:34.429944+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

27 extracted references · 27 canonical work pages

  1. [1]

    A. O. Scheie, Y. Kamiya, H. Zhang, S. Lee, A. J. Woods, M. O. Ajeesh, M. G. Gonzalez, B. Bernu, J. W. Villanova, J. Xing, Q. Huang, Q. Zhang, J. Ma, E. S. Choi, D. M. Pajerowski, H. Zhou, A. S. Sefat, S. Okamoto, T. Berlijn, L. Messio, R. Movshovich, C. D. Batista, and D. A. Ten- nant, Nonlinear magnons and exchange Hamiltonians of the delafossite proxima...

  2. [2]

    A. O. Scheie, M. Lee, K. Wang, P. Laurell, E. S. Choi, D. Pajerowski, Q. Zhang, J. Ma, H. D. Zhou, S. Lee, S. M. Thomas, M. O. Ajeesh, P. F. S. Rosa, A. Chen, V. S. Zapf, M. Heyl, C. D. Batista, E. Dagotto, J. E. Moore, and D. A. Tennant, Spectrum and low-energy gap in triangular quantum spin liquid NaYbSe 2, arXiv preprint arXiv:2406.17773 10.48550/arXiv...

  3. [3]

    T. Xie, A. A. Eberharter, J. Xing, S. Nishimoto, M. Brando, P. Khanenko, J. Sichelschmidt, A. A. Tur- rini, D. G. Mazzone, P. G. Naumov, L. D. Sanjeewa, N. Harrison, A. S. Sefat, B. Normand, A. M. L¨ auchli, A. Podlesnyak, and S. E. Nikitin, Complete field-induced spectral response of the spin-1/2 triangular-lattice an- tiferromagnet csybse2, npj Quantum ...

  4. [4]

    S. Yano, J. Yang, K. Iida, C.-W. Wang, A. G. Man- ning, D. Ueta, and S. Itoh, Spin reorientation and in- terplanar interactions of the two-dimensional triangular- lattice Heisenberg antiferromagnetsh−(Lu,Y)MnO 3 and h−(Lu,Sc)FeO 3, Phys. Rev. B110, 134444 (2024)

  5. [5]

    Rastelli, A

    E. Rastelli, A. Tassi, and L. Reatto, Noncollinear mag- netic order and spin wave spectrum in presence of com- peting exchange interactions, Journal of Magnetism and Magnetic Materials15-18, 357 (1980)

  6. [6]

    Jolicoeur, E

    T. Jolicoeur, E. Dagotto, E. Gagliano, and S. Bacci, Ground-state properties of the s=1/2 heisenberg antifer- romagnet on a triangular lattice, Phys. Rev. B42, 4800 (1990)

  7. [7]

    S. E. Korshunov, Chiral phase of the Heisenberg anti- ferromagnet with a triangular lattice, Phys. Rev. B47, 6165 (1993)

  8. [8]

    Glittum and O

    C. Glittum and O. F. Sylju˚ asen, Arc-shaped structure factor in theJ 1−J2−J3 classical Heisenberg model on the triangular lattice, Phys. Rev. B104, 184427 (2021)

  9. [9]

    Sylju˚ asen and Jens Paaske, Nematic Bond Theory of Heisenberg Helimagnets, Phys- ical Review Letters119, 157202 (2017)

    Michael Schecter and Olav F. Sylju˚ asen and Jens Paaske, Nematic Bond Theory of Heisenberg Helimagnets, Phys- ical Review Letters119, 157202 (2017)

  10. [10]

    O. F. Sylju˚ asen, J. Paaske, and M. Schecter, Interplay between magnetic and vestigial nematic orders in the lay- eredJ 1 −J 2 classical Heisenberg model, Phys. Rev. B99, 174404 (2019)

  11. [11]

    N. D. Mermin and H. Wagner, Absence of Ferro- magnetism or Antiferromagnetism in One- or Two- Dimensional Isotropic Heisenberg Models, Phys. Rev. Lett.17, 1133 (1966)

  12. [12]

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

  13. [13]

    C. L. Henley, Ordering due to disorder in a frustrated vector antiferromagnet, Phys. Rev. Lett.62, 2056 (1989)

  14. [14]

    Chandra, P

    P. Chandra, P. Coleman, and A. I. Larkin, Ising tran- sition in frustrated heisenberg models, Phys. Rev. Lett. 64, 88 (1990)

  15. [15]

    Kawamura, A

    H. Kawamura, A. Yamamoto, and T. Okubo, Z 2-Vortex Ordering of the Triangular-Lattice Heisenberg Antifer- romagnet, Journal of the Physical Society of Japan79, 023701 (2010)

  16. [16]

    Kawamura and S

    H. Kawamura and S. Miyashita, Phase transition of the two-dimensional heisenberg antiferromagnet on the tri- angular lattice, Journal of the Physical Society of Japan 53, 4138 (1984), https://doi.org/10.1143/JPSJ.53.4138

  17. [17]

    Aoyama and H

    K. Aoyama and H. Kawamura, Spin current as a probe of the Z 2-vortex topological transition in the classical Heisenberg antiferromagnet on the triangular lattice, Physical Review Letters124, 047202 (2020)

  18. [18]

    N. D. Andriushin, S. E. Nikitin, Ø. S. Fjellv˚ ag, J. S. White, A. Podlesnyak, D. S. Inosov, M. C. Rahn, M. Schmidt, M. Baenitz, and A. S. Sukhanov, Obser- vation of the spiral spin liquid in a triangular-lattice ma- terial, Nature Communications16, 2619 (2025)

  19. [19]

    Bergman, J

    D. Bergman, J. Alicea, E. Gull, S. Trebst, and L. Ba- lents, Order-by-disorder and spiral spin-liquid in frus- trated diamond-lattice antiferromagnets, Nature Physics 3, 487 (2007)

  20. [20]

    Okumura, H

    S. Okumura, H. Kawamura, T. Okubo, and Y. Mo- tome, Novel Spin-Liquid States in the Frustrated Heisen- berg Antiferromagnet on the Honeycomb Lattice, Jour- nal of the Physical Society of Japan79, 114705 (2010), https://doi.org/10.1143/JPSJ.79.114705

  21. [21]

    Yan and J

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

  22. [22]

    M. G. Gonzalez, A. Fancelli, H. Yan, and J. Reuther, Magnetic properties of the spiral spin liquid and sur- rounding phases in the square lattice XY model, Phys. Rev. B110, 085106 (2024)

  23. [23]

    L. E. Svistov, A. I. Smirnov, L. A. Prozorova, L. N. Demi- anets, and A. Y. Shapiro, Quasi-two-dimensional antifer- romagnet on a triangular lattice RbFe(MoO4)2, Physical Review B67, 094434 (2003)

  24. [24]

    Ishii, S

    R. Ishii, S. Tanaka, K. Onuma, Y. Nambu, M. Toku- naga, T. Sakakibara, N. Kawashima, Y. Maeno, C. Bro- holm, D. P. Gautreaux, J. Y. Chan, and S. Nakat- suji, Successive phase transitions and phase diagrams for 9 the quasi-two-dimensional easy-axis triangular antiferro- magnet Rb4Mn(MoO4)3, Europhysics Letters94, 17001 (2011)

  25. [25]

    Kadowaki, H

    H. Kadowaki, H. Takei, and K. Motoya, Double-Q 120 degrees structure in the Heisenberg antiferromagnet on rhombohedrally stacked triangular lattice LiCrO2, Jour- nal of Physics: Condensed Matter7, 6869 (1995)

  26. [26]

    A. V. Ushakov, D. A. Kukusta, A. N. Yaresko, and D. I. Khomskii, Magnetism of layered chromium sulfides MCrS 2 (M= Li, Na, K, Ag, and Au): A first-principles study, Phys. Rev. B87, 014418 (2013)

  27. [27]

    Mohanty, K

    S. Mohanty, K. M. Ranjith, C. S. Saramgi, Y. Skourski, B. B¨ uchner, H. J. Grafe, and R. Nath, Ground state properties of a spin- 5 2 frustrated triangular lattice an- tiferromagnet NH 4Fe(PO3F)2 (2025), arXiv:2502.06422 [cond-mat.mtrl-sci]