pith. sign in

arxiv: 2605.19831 · v1 · pith:5YPF4JAPnew · submitted 2026-05-19 · ❄️ cond-mat.stat-mech

Ground-state Entropy of the Ising model on a Frustrated lattice

B Sriram Shastry , Bill Sutherland , Fr\'ed\'eric Mila , Afonso Rufino This is my paper

Pith reviewed 2026-05-20 02:11 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords Ising modelground-state entropyShastry-Sutherland latticefrustrated latticeresidual entropyzero-temperature configurationsstatistical mechanics
0
0 comments X

The pith

The 2D Ising model on the Shastry-Sutherland lattice has a non-zero ground-state entropy due to frustration.

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

This paper calculates the ground-state entropy for a two-dimensional Ising model on the Shastry-Sutherland lattice. The lattice geometry creates frustration because competing interactions cannot all be satisfied at once. This leaves a large number of configurations with the same lowest energy, producing residual entropy at zero temperature. The authors also examine a continuous generalization that relaxes one constraint on the allowed ground-state configurations. Such calculations quantify degeneracy in frustrated magnets and clarify their low-temperature thermodynamics.

Core claim

We report the ground-state entropy of a 2-d Ising model on the Shastry-Sutherland lattice. We also study a generalization of this model, where a constraint on the zero temperature allowed configurations is removed continuously.

What carries the argument

The Shastry-Sutherland lattice with Ising spins subject to a specific pattern of frustrated interactions, whose degenerate ground-state manifold is counted via enumeration or transfer-matrix methods.

If this is right

  • The ground states remain degenerate even at absolute zero, producing a finite residual entropy per spin.
  • Relaxing the constraint continuously produces a family of models whose entropy varies smoothly with the relaxation parameter.
  • The degeneracy is a direct geometric consequence of the lattice and interaction pattern rather than an artifact of finite size.
  • Similar counting methods apply to other lattices that share the same local frustration motif.

Where Pith is reading between the lines

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

  • If the entropy is extensive, the number of ground states grows exponentially with system size, a signature of strong frustration.
  • Adding weak quantum fluctuations would likely lift the degeneracy and select a unique ground state or produce a quantum spin liquid.
  • The continuous generalization offers a tunable knob that could be realized in artificial spin-ice arrays or optical-lattice simulators.
  • Comparison with series expansions or tensor-network methods on the same lattice would provide an independent check on the entropy value.

Load-bearing premise

The Shastry-Sutherland lattice geometry and the Ising interaction pattern are taken as given, and the ground-state manifold is assumed to be correctly identified by the enumeration or transfer-matrix method used in the calculation.

What would settle it

An exact count of ground states on a larger finite patch of the lattice, or a Monte Carlo sampling at extremely low temperature, that yields an entropy per site different from the reported value would falsify the result.

Figures

Figures reproduced from arXiv: 2605.19831 by Afonso Rufino, Bill Sutherland, B Sriram Shastry, Fr\'ed\'eric Mila.

Figure 1
Figure 1. Figure 1: The Ising model on the SS lattice. The weights in the Ising energy Eq. (1) are 2α for diagonal bonds, and all others have weight 1. Stretching the bonds leads to other picturizations of this lattice as in Fig. (2) parallel spins residing on the diagonal bond, which are missing at α > 1. We refer to the fraction of these new ferromagnetic diagonal type configurations as nfmd, nfmd = 1 N X <l,m> ⟨(1 + σlσm)⟩… view at source ↗
Figure 2
Figure 2. Figure 2: Variants of the SS lattice. Left Archimedean limit where equilateral triangle and squares with equal bond length tile the 2-d plane. This is one of the 11 Archimedean lattices [11].Right The orthogonal dimer version where the diagonal bonds of Fig. (1) are shrunk relative to the sides of the squares. This picture is relevant to the measured lattice structure of SrCu2(BO3)2 [3, 4]. Since the diagonal bonds … view at source ↗
Figure 3
Figure 3. Figure 3: Two rows of types (A) and (B) with diagonals slanting distinctly. We represent the allowed matrix elements in the notation W (A) (σ ′ i , σ′ i+1|σi , σi+1) with odd i, and W(B) for even i. 2 Transfer matrix for calculating ground-state entropy We now set up a transfer matrix for computing the entropy of the Ising model on the SS-lattice at T=0. It is also useful to consider a generalization where configura… view at source ↗
Figure 4
Figure 4. Figure 4: The configurations at each (A) type square W (A) (σ ′ i , σ′ i+1|σi , σi+1), where σ are the spins on the lower side and σ ′ are the upper side. The configuration ○3 is missing in the entropy calculation for α > 1, while all other configurations are common to the two problems. A magnitude r is assigned to the diagram ○3 , and the rest are assigned 1. Varying r from 0 to 1 gives a generalized model, where o… view at source ↗
Figure 5
Figure 5. Figure 5: Illustration of the CTMRG implementation of the SS lattice Ising model. (a) The partition function of the system is written exactly as the contraction of a square network of fourth-order tensors ai1,i2,i3,i4 . (b) The contraction of all tensors surrounding a site is approximated by the product of finite corner (Ci) and edge (Ti) transfer matrices. formulation of the partition function is obtained by lettin… view at source ↗
Figure 6
Figure 6. Figure 6: Illustration of the exact calculation of the matrix element ⟨ΨA|TB|ΨA⟩ in Tensor Network diagrammatic notation [17]. It is also possible to calculate ⟨ΨA|TB|ΨA⟩ exactly for infinite L, by ex￾ploiting the tensor product structure of |ΨA⟩. Given the formula for TA as a tensor product of two-site operators (12), the leading eigenvector can also be written as the tensor product |ΨA⟩ = ⊗ L/2 i=1|ϕ 2i−1,2i A ⟩, … view at source ↗
Figure 7
Figure 7. Figure 7: Entropy (Σ), correlation length of spin-spin correlations (ξ) and fraction of ferromagnetic bonds (nfmd) versus the parameter r introduced in Eq. (7, 8), calculated with CTMRG. At r = 0 we suppress configurations with parallel spins on the diagonals (see Fig. (4)). This constraint gives all configurations for α>1, and by an elementary calculation gives Σ= 1 2 log 2. This is also the magnitude of the interc… view at source ↗
read the original abstract

We report the ground-state entropy of a 2-d Ising model on the Shastry-Sutherland lattice. We also study a generalization of this model, where a constraint on the zero temperature allowed configurations is removed continuously.

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

Summary. The manuscript reports the ground-state entropy of the 2D Ising model on the Shastry-Sutherland lattice, obtained as (1/N) log W where W is the number of minimizing configurations in the ground-state manifold. It also examines a continuous generalization in which a constraint on the zero-temperature allowed configurations is removed.

Significance. If the reported entropy value is accurate and the manifold is correctly enumerated, the result supplies a concrete benchmark for residual entropy in a frustrated 2D Ising system on the Shastry-Sutherland lattice, which is realized in materials such as SrCu2(BO3)2. The continuous generalization provides a controlled way to study the lifting of degeneracy and the evolution of the entropy, which could connect to broader questions of constrained systems and their thermodynamics.

major comments (2)
  1. [Enumeration / transfer-matrix section] The identification of the ground-state manifold (central to the entropy claim) rests on an enumeration or transfer-matrix procedure whose completeness is not cross-checked by independent brute-force enumeration on small clusters. Any systematic under- or over-counting of W propagates directly into the reported residual entropy per site.
  2. [Results] No error analysis, finite-size scaling, or convergence checks for the entropy value are supplied, which is required to substantiate the numerical result given that the central claim is the specific value of (1/N) log W.
minor comments (1)
  1. [Abstract] The abstract states the result without indicating the numerical value, method, or any supporting detail, which is atypical and reduces immediate accessibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address the two major comments point by point below and describe the changes we will implement.

read point-by-point responses

  1. Referee: [Enumeration / transfer-matrix section] The identification of the ground-state manifold (central to the entropy claim) rests on an enumeration or transfer-matrix procedure whose completeness is not cross-checked by independent brute-force enumeration on small clusters. Any systematic under- or over-counting of W propagates directly into the reported residual entropy per site.

    Authors: We agree that an independent cross-check on small clusters would increase confidence in the completeness of the ground-state manifold. In the revised manuscript we will add a dedicated subsection presenting brute-force enumeration results for small finite clusters (compatible with periodic boundary conditions on the Shastry-Sutherland lattice) and direct comparison of the counted W with the transfer-matrix output for the same sizes. This verification will be shown explicitly for several small system sizes. revision: yes

  2. Referee: [Results] No error analysis, finite-size scaling, or convergence checks for the entropy value are supplied, which is required to substantiate the numerical result given that the central claim is the specific value of (1/N) log W.

    Authors: We acknowledge the absence of these analyses in the original submission. The revised manuscript will include a new figure and accompanying text that reports the entropy per site for a sequence of increasing system sizes, together with a finite-size scaling extrapolation to the thermodynamic limit and a discussion of numerical uncertainties arising from the enumeration procedure. revision: yes

Circularity Check

0 steps flagged

No circularity detected; derivation is a direct count of minimizing configurations

full rationale

The paper reports the ground-state entropy via enumeration or transfer-matrix counting of configurations that minimize all local bond energies on the given Shastry-Sutherland lattice geometry with the stated antiferromagnetic couplings. The abstract and description contain no equations, no fitted parameters renamed as predictions, and no self-citation chains that bear the central claim. The residual entropy is obtained directly as (1/N) log W where W is the counted degeneracy; this is independent of the inputs once the lattice and Hamiltonian are fixed. The method is a standard, externally verifiable computational procedure for frustrated Ising models and does not reduce to its own assumptions by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review supplies no explicit free parameters, axioms, or invented entities; the standard Ising Hamiltonian on the named lattice is presupposed.

pith-pipeline@v0.9.0 · 5561 in / 1017 out tokens · 35765 ms · 2026-05-20T02:11:25.000752+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

19 extracted references · 19 canonical work pages

  1. [1]

    B. S. Shastry and B. Sutherland, Exact Ground State of a Quantum- Mechanical Antiferromagnet, Physica108B, 1069 (1981)

    work page 1981

  2. [2]

    Albrecht and F

    M. Albrecht and F. Mila, First-order transition between magnetic order and valence bond order in a 2D frustrated Heisenberg model, Europhys. Lett. 34 , 145 (1996)

    work page 1996

  3. [3]

    R. W. Smith and D. A. Keszler, Synthesis, structure, and properties of the orthoborateSrCu 2(BO3)2, J. Solid State Chem.93, 430 (1991)

    work page 1991

  4. [4]

    Kageyama, K

    H. Kageyama, K. Yoshimura, N. V. Mushnikov, K. Onizuka, M. Kato, K. Kosuge, C. P. Slichter, T. Goto, and Y. Ueda, Exact Dimer Ground State and Quantized Magnetization Plateaus in the Two-Dimensional Spin SystemSrCu 2(BO3)2, Phys. Rev. Letts.82, 3168 (1999)

    work page 1999

  5. [5]

    Miyahara and K

    S. Miyahara and K. Ueda, Exact Dimer Ground State of the Two Di- mensional Heisenberg Spin SystemSrCu 2(BO3)2 , Phys. Rev. Letts. 82, 3701 (1999)

    work page 1999

  6. [6]

    Koga and N

    A. Koga and N. Kawakami,Quantum Phase Transitions in the Shastry- Sutherland Model forSrCu 2(BO3)2 , Phys. Rev. Lett. 84, 4461, (2000)

    work page 2000

  7. [7]

    Chung, J

    C. Chung, J. Marston and S. Sachdev, Quantum phases of the Shastry- Sutherland antiferromagnet: Application toSrCu 2(BO3)2, Phys. Rev. 64, 134407 (2001). REFERENCES17

    work page 2001

  8. [8]

    Streˇ cka, Exact solution of the spin-1/2 Ising model on the Shastry- Sutherland (orthogonal-dimer) lattice, Phys

    J. Streˇ cka, Exact solution of the spin-1/2 Ising model on the Shastry- Sutherland (orthogonal-dimer) lattice, Phys. Letts.349, 505 (2006)

    work page 2006

  9. [9]

    Rousochatzakis, A

    I. Rousochatzakis, A. M. L´ auchli, and R. Moessner, Quantum mag- netism on the Cairo pentagonal lattice, Phys. Rev. B85, 104415 (2012)

    work page 2012

  10. [10]

    Johannes Kepler,Harmonice Mundi, Linz (Austria): Johann Planck, (1619)

    work page

  11. [11]

    Branko Gr¨ unbaum and G. C. Shephard,Tilings and Patterns, Dover Publications, 2016. (ISBN-13978-0486469812)

    work page 2016

  12. [12]

    Barma and B

    M. Barma and B. S. Shastry, Classical equivalents of one-dimensional quantum-mechanical systems, Phys. Rev. B18, 3351 (1978), see Fig. 1

    work page 1978

  13. [13]

    J., Corner transfer matrices

    Baxter, R. J., Corner transfer matrices. Physica A Statistical Mechanics and its Applications 106, 18–27 (1981)

    work page 1981

  14. [14]

    and Okunishi, K., Corner Transfer Matrix Renormalization Group Method

    Nishino, T. and Okunishi, K., Corner Transfer Matrix Renormalization Group Method. J. Phys. Soc. Jpn. 65, 891–894 (1996)

    work page 1996

  15. [15]

    and Mila, F

    Vanhecke, B., Colbois, J., Vanderstraeten, L., Verstraete, F. and Mila, F. Solving frustrated Ising models using tensor networks. Phys. Rev. Research 3, 013041 (2021)

    work page 2021

  16. [16]

    In a finite T study of the Ising model, the factorrin the weight factor Wcan be (somewhat artificially) generated by adding energy−logr×T for configurations with parallel spins on a diagonal

    work page

  17. [17]

    A line connecting two tensors stands for a contracted index

    In the Tensor Network diagrammatic notation, tensors are represented as squares or circles, while their indices are represented as lines. A line connecting two tensors stands for a contracted index

    work page

  18. [18]

    Wright and B

    John G. Wright and B. S. Shastry, ”DiracQ: A Quantum Many- Body Physics Package”, arXiv:1301.4494[cond-mat.str-el] (2013), Journal of Open Research Software, 3(1) e13 (2015); DOI: https://doi.org/10.5334/jors.cb. REFERENCES18

    work page doi:10.5334/jors.cb 2013

  19. [19]

    Sutherland and B

    B. Sutherland and B. S. Shastry, ‘Exact Solution of a Large Class of Quantum Systems Exhibiting Ground State Singularities, J. Stat. Phys. 33, 477 (1983)

    work page 1983