S=1 Shastry–Sutherland Model

• The S=1 Shastry–Sutherland model is a frustrated quantum-spin system on a decorated square lattice that interpolates between dimer and magnetically ordered states via intermediate phases.
• Its anisotropic and field-dependent variants reveal exotic phases such as quantum paramagnets, dual spin supersolids, and unique m=1/6 plateaus not present in S=1/2 systems.
• Advanced numerical techniques like DMRG and CMFT have mapped its complex phase boundaries, highlighting subtle symmetry effects and nontrivial spin transport behaviors.

The $S=1$ Shastry–Sutherland model denotes a family of frustrated quantum-spin systems defined on the two-dimensional Shastry–Sutherland lattice (SSL), a square lattice with an orthogonal covering of diagonal bonds. In its simplest isotropic Heisenberg form, it interpolates between a product-dimer regime and magnetically ordered states through intermediate phases that are absent in the classical limit. In anisotropic and field-dependent variants, the enlarged on-site Hilbert space of spin one further stabilizes phases not available for $S=1/2$, including a quantum paramagnet, two distinct spin supersolids, an “offspring” $m=1/6$ plateau, and chirality-split magnons in a collinear Néel phase. Closely related generalized constructions also admit exact trimer-singlet product ground states for integer spin (Wu et al., 23 May 2026, Siam et al., 11 Jul 2025, Su et al., 2014, Richter et al., 2012).

1. Lattice, conventions, and model families

The SSL is obtained by placing alternating diagonal bonds on a square lattice so that orthogonal dimers live on two sets of diagonals that do not touch each other. The square plaquettes without a diagonal bond are the “empty plaquettes,” which form two checkerboard sublattices. In the isotropic $S=1$ Heisenberg study, the Hamiltonian is

$$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$

with $J$ the nearest-neighbor inter-dimer coupling on square edges and $J'$ the diagonal intra-dimer coupling; energy units are set by $J'=1$ (Wu et al., 23 May 2026).

For $S=1$, the literature uses several closely related but not identical conventions. In the easy-axis magnon study, the model is a generalized SSL Heisenberg antiferromagnet with single-ion anisotropy,

$$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$

where $S=1/2$0 labels intra-dimer bonds, $S=1/2$1 and $S=1/2$2 label two inequivalent inter-dimer directions, and $S=1/2$3 is easy-axis. In the strongly anisotropic field-driven study, the Hamiltonian is an $S=1/2$4 $S=1/2$5 model with single-ion anisotropy and longitudinal field,

$S=1/2$6

These differing notations are a persistent source of confusion; the diagonal-bond coupling is denoted by $S=1/2$7 in the isotropic Heisenberg work, but by $S=1/2$8 in the easy-axis magnon work (Siam et al., 11 Jul 2025, Su et al., 2014).

Model variant	Hamiltonian ingredients	Distinctive $S=1/2$9 phenomena
Isotropic Heisenberg SSL	$m=1/6$0, $m=1/6$1, $m=1/6$2	dimer, MPD, spiral, NAF
Easy-axis generalized SSL	$m=1/6$3, $m=1/6$4, $m=1/6$5, $m=1/6$6	Néel (I), chirality-split magnons, SSE/SNE
Ising-like $m=1/6$7 SSL in field	$m=1/6$8, $m=1/6$9, $S=1$0, $S=1$1	QP, SS1, SS2, $S=1$2, $S=1$3, $S=1$4 plateaus
Modified trimer SSL	$S=1$5, $S=1$6 on decorated lattice	exact trimer-singlet product phase

At the isotropic point $S=1$7, each $S=1$8 bond hosts an exact $S=1$9 singlet, and the bond expectation value is

$$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$0

The corresponding product state is the natural $$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$1 analogue of the dimer limit (Wu et al., 23 May 2026).

2. Zero-field isotropic Heisenberg phase diagram

For the isotropic $$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$2 SSL Heisenberg model, density matrix renormalization group (DMRG) and cluster mean-field theory (CMFT) identify the zero-field phase sequence

$$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$3

as $$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$4 increases (Wu et al., 23 May 2026).

The transition from the dimer phase to the MPD phase occurs at $$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$5 and is first order. It is located from sharp features in $$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$6 and jumps in the bipartite entanglement entropy

$$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$7

The MPD-to-spiral transition is likewise first order at $$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$8, again diagnosed by discontinuities in $$H = J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j + J'\sum_{\langle\!\langle i,j\rangle\!\rangle'}\mathbf{S}_i\cdot\mathbf{S}_j, \qquad g=\frac{J}{J'},$$9 and $J$0. The spiral-to-NAF transition at $J$1 is continuous and corresponds to commensurate locking of an incommensurate ordering wavevector to $J$2; the energy derivative shows no clear singularity within finite-size resolution, while the magnetic Bragg peak evolves continuously and $J$3 shows a sharp drop followed by a plateau (Wu et al., 23 May 2026).

The numerical evidence comes from SU(2)-symmetric DMRG on finite cylinders of widths $J$4 and lengths up to $J$5, retaining up to 8000 states with truncation errors below $J$6. Open clusters are also used to expose symmetry-broken bond patterns that are partially hidden on cylinders. CMFT calculations employ $J$7 clusters with $J$8 or $J$9 and a DMRG impurity solver implemented via ITensor, with truncation errors below $J'$0 (Wu et al., 23 May 2026).

Placed in a broader $J'$1-$J'$2 context, the $J'$3 sequence differs qualitatively from both the $J'$4 and classical limits. The global diagram assembled from $J'$5, $J'$6, $J'$7, and $J'$8 shows that the spiral regime expands with increasing $J'$9, while the dimer and MPD windows shrink and ultimately vanish in the classical limit. This supports the interpretation that quantum-disordered phases are progressively suppressed as $J'=1$0 increases (Wu et al., 23 May 2026).

3. Microscopic characterization of the MPD, spiral, and Néel regimes

The primary local diagnostic is the bond energy

$J'=1$1

In the dimer phase, $J'=1$2 is strongly negative on the diagonal $J'=1$3 bonds and much weaker on the inter-dimer $J'=1$4 bonds. In the MPD phase, DMRG finds strong intradimer correlations together with weak tetramerization on the empty plaquettes. This tetramerization appears as a small but finite alternation of the four edge-bond energies around an empty plaquette, with two symmetry-related patterns residing on the two checkerboard plaquette sublattices (Wu et al., 23 May 2026).

A convenient MPD order parameter is

$J'=1$5

on the two empty-plaquette sublattices. On open clusters, $J'=1$6 is small but finite and its sign selects one of the two degenerate plaquette patterns. On cylinders, the superposition of degenerate patterns strongly suppresses the apparent order parameter. This clarifies a common misconception: the MPD phase is not a pure plaquette state. It combines robust dimers with only weak empty-plaquette tetramerization, and is therefore distinct from the $J'=1$7 empty-plaquette phase (Wu et al., 23 May 2026).

Real-space spin correlations,

$J'=1$8

further separate the phases. In the MPD regime, correlations decay exponentially and display a sawtooth-like modulation along the $J'=1$9 direction, consistent with tetramerization on empty plaquettes. In the spiral regime, correlations show power-law decay with incommensurate oscillations. In the NAF phase, they exhibit power-law decay without oscillations, consistent with commensurate $S=1$0 order on finite cylinders (Wu et al., 23 May 2026).

The static structure factor,

$S=1$1

is broad in the dimer and MPD phases, develops a sharp peak at $S=1$2 in the spiral regime, and locks to $S=1$3 in the NAF phase. The spiral ordering wavevector

$S=1$4

moves continuously toward $S=1$5 as $S=1$6 increases. The pitch angle is $S=1$7 and the real-space period along $S=1$8 is $S=1$9. At $$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$0, CMFT finds an $$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$1 magnetic unit cell, consistent with a commensurate approximation to an incommensurate spiral with $$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$2 (Wu et al., 23 May 2026).

4. Easy-axis Néel (I) order, altermagnetic magnons, and spin transport

A distinct $$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$3 SSL regime arises when easy-axis single-ion anisotropy stabilizes a collinear Néel (I) state. In the generalized model with antiferromagnetic $$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$4, $$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$5, $$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$6 and $$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$7, the representative parameters

$$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$8

place the system deep inside a collinear $$H=\sum_{(i,j)}J_{ij}\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2,$$9-axis antiferromagnet with a four-site magnetic unit cell. Sites 1 and 3 are spin-up and sites 2 and 4 are spin-down, and the ordered state belongs to the plane group $S=1/2$00 (Siam et al., 11 Jul 2025).

Linear spin-wave theory around this state uses a four-sublattice Holstein–Primakoff expansion and produces an $S=1/2$01 bosonic Bogoliubov–de Gennes Hamiltonian. Because there are four spins per unit cell, there are four positive-energy magnon bands. For the representative parameters, the spectrum separates into a lower pair and an upper pair, with each pair consisting of two bands of opposite chirality. Magnon chirality is labeled by

$S=1/2$02

so that the bands are denoted $S=1/2$03 (Siam et al., 11 Jul 2025).

The nontrivial result is that opposite-chirality magnon bands are nondegenerate throughout the Brillouin zone without spin–orbit coupling, Dzyaloshinskii–Moriya interactions, or an external field. The mechanism is symmetry-based. In the Néel (I) SSL, opposite-spin sublattices are related by $S=1/2$04 rather than by inversion centered between opposite-spin sites, and inversion centers lie at mid-edges and map each sublattice to itself. The chirality splitting alternates sign along orthogonal momentum directions, a magnonic realization of altermagnetic alternation (Siam et al., 11 Jul 2025).

Transport is treated within a relaxation-time Kubo formalism for the magnon spin current driven by a thermal gradient. The resulting longitudinal spin Seebeck and transverse spin Nernst responses are finite because the two circularly polarized magnons have different thermal populations and group velocities when $S=1/2$05. Using $S=1/2$06 and $S=1/2$07, the study finds that $S=1/2$08 is maximal near $S=1/2$09 and minimal near $S=1/2$10, while $S=1/2$11 vanishes when the thermal gradient is aligned with crystal axes and peaks near $S=1/2$12. When $S=1/2$13, the restored $S=1/2$14 symmetry enforces $S=1/2$15 and the spin Nernst coefficient vanishes for all $S=1/2$16. By contrast, the intrinsic thermal Hall conductivity remains negligible or zero because Berry curvature is symmetry-forbidden or cancels pairwise in this setting (Siam et al., 11 Jul 2025).

5. Uniaxial anisotropy, magnetic field, and plateau physics

In the strongly Ising-like $S=1/2$17 $S=1/2$18 SSL with single-ion anisotropy, stochastic series expansion quantum Monte Carlo on $S=1/2$19 tori with $S=1/2$20–$S=1/2$21 reveals a much richer field-driven structure than in the isotropic Heisenberg problem. The study focuses on $S=1/2$22, with ground-state behavior accessed by low temperatures using $S=1/2$23 and, in some scans, $S=1/2$24 (Su et al., 2014).

At zero field, four phases appear as functions of $S=1/2$25 and $S=1/2$26: an Ising antiferromagnet with ordering vector $S=1/2$27, a transverse XY-ordered phase labeled “SF,” a large-$S=1/2$28 triplet-dimer phase, and a quantum paramagnet stabilized by sufficiently large $S=1/2$29. The SF–QP transition is continuous in the $S=1/2$30 universality class, with representative critical values $S=1/2$31 at $S=1/2$32 and $S=1/2$33 at $S=1/2$34 (Su et al., 2014).

Under a longitudinal field, the magnetization process exhibits several plateau and supersolid phases. For small $S=1/2$35, the sequence is

$S=1/2$36

With increasing $S=1/2$37, SS1 and PL1 disappear and a $S=1/2$38 plateau PL2 emerges. For $S=1/2$39, a $S=1/2$40 plateau PL3 appears below PL2, and a second supersolid SS2 is stabilized between them. For very large $S=1/2$41, where the zero-field state is triplet-dimer-like, a distinct high-field $S=1/2$42 plateau PL4 appears; unlike PL1, PL4 is a quantum-disordered dimer plateau without Bragg peaks in either longitudinal or transverse structure factors (Su et al., 2014).

The $S=1/2$43 and $S=1/2$44 plateaus are particularly characteristic of the $S=1/2$45 case. In the Ising limit, PL2 is a period-3 stripe state with ordering wavevector

$S=1/2$46

PL3 is obtained from the same parent stripe pattern by replacing each $S=1/2$47 stripe by an $S=1/2$48 stripe, preserving the period-3 order. This “offspring” mechanism is impossible for $S=1/2$49 because the $S=1/2$50 state does not exist. SS2 inherits the same period-3 longitudinal order while developing finite spin stiffness $S=1/2$51, thereby combining diagonal and off-diagonal long-range order (Su et al., 2014).

This variant also resolves another common ambiguity in the literature: the large-$S=1/2$52 “dimer” phase need not be a singlet-dimer phase. In the present sign convention, the transverse exchange is ferromagnetic, so the large-$S=1/2$53 state is described as a triplet-dimer analogue of the canonical $S=1/2$54 singlet-dimer phase (Su et al., 2014).

6. Exact generalized constructions, comparisons, and open directions

A related but distinct integer-spin extension replaces each SSL diagonal by an equilateral triangle, producing a decorated square lattice of coupled spin trimers. For integer spin, and in particular for $S=1/2$55, the local trimer has a unique singlet ground state with energy

$S=1/2$56

Because the inter-trimer couplings satisfy a balance condition, the product

$S=1/2$57

of local trimer singlets is an exact eigenstate of the full lattice Hamiltonian and, within a finite parameter window, the exact ground state. For $S=1/2$58, exact diagonalization with periodic boundary conditions gives $S=1/2$59 and $S=1/2$60 for $S=1/2$61 on $S=1/2$62 sites, while a rigorous inner interval is

$S=1/2$63

Within this trimer-singlet product phase, the ground-state energy per site is

$S=1/2$64

independent of $S=1/2$65 (Richter et al., 2012).

Although this decorated-trimer model is not the standard SSL Heisenberg model, it is part of the broader $S=1/2$66 Shastry–Sutherland landscape because it preserves the underlying orthogonal-dimer logic while changing the local building block from dimers to trimers. It also exhibits robust magnetization plateaus at $S=1/2$67 and $S=1/2$68, and a ferrimagnetic phase for sufficiently negative $S=1/2$69 (Richter et al., 2012).

Across the literature, several comparisons are now clear. First, $S=1/2$70 physics is not obtained by merely weakening quantum fluctuations relative to $S=1/2$71. The isotropic model supports an MPD phase and an incommensurate spiral between the dimer and Néel limits; the anisotropic field-driven model supports a quantum paramagnet, two supersolids, and an offspring $S=1/2$72 plateau; and the easy-axis ordered regime supports nonrelativistic chirality-split magnon bands without SOC or field (Wu et al., 23 May 2026, Su et al., 2014, Siam et al., 11 Jul 2025). Second, the role of anisotropy is decisive: easy-axis single-ion anisotropy stabilizes Néel (I) order and altermagnetic magnon transport, whereas easy-plane single-ion anisotropy stabilizes the quantum paramagnet and reshapes the field-induced plateau hierarchy (Siam et al., 11 Jul 2025, Su et al., 2014).

The experimental context is correspondingly broad. Materials mentioned in this literature include rare-earth-based $S=1/2$73 ($S=1/2$74–$S=1/2$75), $S=1/2$76, $S=1/2$77, and Pr-based melilites such as $S=1/2$78 and $S=1/2$79. Proposed signatures include short-range scattering without sharp Bragg peaks in the MPD regime, a sharp incommensurate peak at $S=1/2$80 in the spiral regime, conventional $S=1/2$81 order in the NAF regime, angle-dependent spin Seebeck and spin Nernst responses in the Néel (I) regime, and plateau structures at $S=1/2$82, $S=1/2$83, and $S=1/2$84 in the anisotropic field-driven problem (Wu et al., 23 May 2026, Siam et al., 11 Jul 2025, Su et al., 2014).

Open questions identified in these works include the determination of precise exchange ratios in larger-$S=1/2$85 SSL compounds, the role of exchange anisotropy and interlayer couplings, finite-temperature crossovers in MPD and spiral regimes, possible field-tuned versions of the MPD and spiral phases, interaction and disorder effects beyond linear spin-wave theory, and a quantitative mapping of the full $S=1/2$86 phase diagram in the altermagnetic Néel (I) regime (Wu et al., 23 May 2026, Siam et al., 11 Jul 2025).