Spin Supersolid States

• Spin supersolid states are defined by the simultaneous presence of staggered longitudinal spin modulation and coherent transverse order, breaking both lattice and continuous symmetries.
• They emerge from a competition between interactions that favor density-wave order and the delocalizing effects of transverse exchanges, with evidence from one-dimensional to three-dimensional systems.
• Diagnostic techniques such as structure-factor analysis and spin stiffness measurements confirm the coexistence of diagonal and off-diagonal long-range order in these phases.

Spin supersolid states are magnetic phases in which two orders coexist microscopically: a longitudinal spin modulation that breaks lattice translational symmetry, and a transverse coherent spin component that breaks a continuous spin-rotation symmetry. In the standard Matsubara–Matsuda mapping, this is the magnetic analogue of a bosonic supersolid, with $S_i^z$ playing the role of density and $S_i^+$ the role of boson creation; the spin supersolid is therefore the coexistence of diagonal and off-diagonal order rather than a mere interpolation between a magnetization plateau and a canted antiferromagnet (Huang et al., 5 Jan 2026). The subject now spans one-dimensional benchmark models, frustrated triangular-lattice quantum magnets, metallic rare-earth systems, and spin-orbit-coupled condensates, with the strongest material evidence concentrated in easy-axis triangular antiferromagnets and their metallic extensions (Rossini et al., 2011, Chen et al., 2024, Xu et al., 25 Mar 2026).

1. Concept and symmetry structure

The defining criterion is simultaneous breaking of a discrete lattice symmetry and a continuous internal symmetry. In the conventional spin-$\tfrac12$ triangular-lattice setting, the solid component is a three-sublattice modulation in $S^z$ at the $K$ points of the Brillouin zone, while the superfluid component is transverse order in the $xy$ plane associated with spontaneous breaking of spin $U(1)$ symmetry about the field axis (Huang et al., 5 Jan 2026). In one-dimensional realizations, the same logic applies with a different ordering wavevector: the spin-supersolid phase in the anisotropic spin-1 Heisenberg chain is defined by finite staggered longitudinal order at $q=\pi$ together with finite spin stiffness (Rossini et al., 2011).

This distinction separates spin supersolids from nearby phases. A UUD plateau, or more generally a commensurate collinear phase, has the solid component but lacks transverse phase coherence. A conventional spin superfluid or canted easy-plane state has broken $U(1)$ symmetry but no density-wave-like modulation of $S^z$. The supersolid requires both. The review literature formulates this directly in terms of structure-factor order parameters,

$S_i^+$0

with both finite in the supersolid regime (Huang et al., 5 Jan 2026).

The notion also extends beyond dipolar transverse order. In spin-1 triangular systems with strong single-ion anisotropy, the off-diagonal component can be quadrupolar rather than dipolar, producing a nematic supersolid characterized by simultaneous longitudinal order and ferroquadrupolar $S_i^+$1-breaking order (Huang et al., 5 Jan 2026). A further generalization appears in the triangular-lattice spin-1 $S_i^+$2 model, where the supersolid breaks lattice translation and an $S_i^+$3 spin symmetry, yielding an “$S_i^+$4-supersolid” with two Goldstone modes rather than one (Wang et al., 25 Sep 2025).

2. Microscopic routes to spin supersolidity

A recurring microscopic mechanism is competition between interactions that favor longitudinal density-wave order and terms that preserve transverse phase coherence. In the spin-1 anisotropic Heisenberg chain,

$S_i^+$5

with $S_i^+$6 and $S_i^+$7, the XY exchange delocalizes spin excitations while the Ising and single-ion anisotropies favor staggered longitudinal order; in the large-$S_i^+$8 regime the model maps to an effective spin-$S_i^+$9 XX chain in a field, making the supersolid interpretable as mobile quasiparticles propagating on an ordered background (Rossini et al., 2011).

In two-dimensional frustrated magnets the canonical minimal model is the easy-axis triangular-lattice XXZ antiferromagnet,

$\tfrac12$0

or equivalent notations with $\tfrac12$1 and $\tfrac12$2 (Xu et al., 2024, Chi et al., 2024). Frustration stabilizes three-sublattice longitudinal order, easy-axis anisotropy favors the solid component, and field tuning drives the familiar sequence

$\tfrac12$3

with the supersolids flanking the UUD plateau in the weak-easy-axis case (Chi et al., 2024, Xiang et al., 15 Apr 2025). In the strongly Ising regime of $\tfrac12$4, correlations stabilize a high-field $\tfrac12$5 ($\tfrac12$6-coplanar) state rather than the V state of simple cluster mean-field theory (Xu et al., 2024).

Three-dimensional routes are now also established. In the realistic model for $\tfrac12$7-AgNiO$\tfrac12$8, a low-field supersolid emerges continuously from a collinear stripe state by softening of a transverse magnon mode at $\tfrac12$9, i.e. by magnon Bose condensation on top of a crystalline background (Seabra et al., 2010). In EuCo$S^z$0Al$S^z$1, by contrast, the effective supersolid is metallic and three-dimensional: interlayer ferromagnetic and intralayer antiferromagnetic RKKY exchange, together with sizable interlayer dipolar anisotropy, generate a stacked-triangular easy-axis model whose zero-field ground state is already a Y supersolid and whose finite-field sequence is Y $S^z$2 UUD $S^z$3 V $S^z$4 polarized (Xi et al., 25 Mar 2026).

The same conceptual structure appears in mixed-spin models and cold-atom analogues. Coupled alternating spin-1 and spin-$S^z$5 chains support a stripe supersolid in which stripe order in longitudinal magnetization coexists with transverse coherence; in the anisotropic limits analyzed in detail, repulsive intrachain interactions are necessary for stripe supersolidity with uniform hoppings, while in the anisotropic-hopping case coupled Ising chains already suffice as a minimal model (Heydarinasab et al., 2017). In spin-orbit-coupled Bose gases, stripe supersolids and supersolid-like crystalline states arise by interference among multiple momentum minima, although those works often emphasize finite-current branches, superstripes, or self-trapped supersolid-like solitons rather than bulk magnetic supersolids (Lyu et al., 2023, Adhikari, 2021, Wang et al., 2019, Kaur et al., 2022).

3. Static diagnostics and order parameters

The most robust identification strategy uses one observable for the solid component and another for the superfluid component. In the spin-1 chain, the longitudinal structure factor

$S^z$6

defines the solid order parameter

$S^z$7

while the off-diagonal sector is quantified by the spin stiffness

$S^z$8

obtained from the curvature of the twisted-boundary ground-state energy (Rossini et al., 2011). The phase classification is then exact in the operational sense: solid if $S^z$9 and $K$0, superfluid if $K$1 and $K$2, and supersolid if both are nonzero (Rossini et al., 2011).

On the triangular lattice, the same structure is used in momentum space at $K$3. In DMRG studies with impurities, the diagonal and transverse channels are measured as

$K$4

so that the supersolid is the regime where both remain finite (Huang et al., 3 Sep 2025). Tensor-network work on $K$5 instead uses explicit three-sublattice order parameters,

$K$6

$K$7

with supersolidity identified by $K$8 and $K$9 (Xu et al., 2024).

In $xy$0, the solid order is encoded in a complex three-sublattice order parameter,

$xy$1

while the superfluid component is the coherent transverse order $xy$2 (Gao et al., 2024). In the spin-1 $xy$3 model, the solid component is carried by three-sublattice order in $xy$4, whereas the superfluid component is an ordered $xy$5 sector; the supersolid is therefore diagnosed by simultaneous $xy$6-point structure in both channels (Wang et al., 25 Sep 2025).

These diagnostics also clarify a methodological point. Longitudinal order can often be extracted reliably with open boundaries or finite cylinders, but stiffness cannot. In the one-dimensional chain, periodic boundary conditions are indispensable because a twist can be gauged away under open boundaries, making $xy$7 trivially; this is precisely why variational MPS with periodic boundary conditions were necessary to directly characterize the supersolid rather than infer it from magnetization profiles (Rossini et al., 2011).

4. Collective modes and dynamical signatures

The most stable dynamical hallmark of a spin supersolid is a Goldstone mode tied to the broken continuous symmetry. In the easy-axis triangular-lattice supersolids, the transverse dynamical structure factor develops a gapless branch at the $xy$8 points; this is seen in tensor-network spectra for $xy$9, in DMRG/TDVP spectra for the Y and V supersolids of the easy-axis triangular antiferromagnet, and in neutron and tensor-network studies of $U(1)$0 (Xu et al., 2024, Huang et al., 3 Sep 2025, Gao et al., 2024). The presence of this mode distinguishes supersolids sharply from the UUD plateau, whose transverse excitations are ordinary gapped magnons (Huang et al., 3 Sep 2025).

Beyond the Goldstone sector, the low-energy spectrum is phase dependent. In the Y supersolid of the easy-axis triangular model, tensor-network calculations resolve two low-energy branches: a true Goldstone branch and a second, weakly gapped branch. At zero field there are two nearly degenerate rotons near $U(1)$1, but under field the roton associated with the Goldstone branch disappears quickly, whereas the roton in the gapped branch persists across the Y phase (Chi et al., 2024). This led to a specific conclusion with broader interpretive value: rotons are not universal signatures of spin supersolidity. The V supersolid retains the Goldstone mode but lacks the low-energy roton structure characteristic of the Y phase (Chi et al., 2024).

In $U(1)$2, the supersolid dynamics are even richer. The low-energy spectrum exhibits a “double magnon-roton” structure: one branch contains the true $U(1)$3 Goldstone mode and a roton near $U(1)$4, while the other contains a pseudo-Goldstone mode and another roton. The pseudo-Goldstone gap is attributed to a sixfold anisotropy of the solid order parameter via order-by-quantum-disorder, whereas the true Goldstone mode tracks the phase of the transverse superfluid component (Gao et al., 2024). The same work connects these low-energy roton branches to the large sub-Kelvin entropy and giant magnetocaloric response of the material (Gao et al., 2024).

For $U(1)$5, tensor-network excitation calculations add a different dynamical signature: in the zero-field Y supersolid, the spectrum contains nearly dispersionless branches at integer multiples of the Ising scale,

$U(1)$6

with visible features at $U(1)$7, $U(1)$8, and $U(1)$9. Real-space inspection of the excited tensors shows that the $q=\pi$0 sector has nonlocal two-spin-flip character, supporting its interpretation as a multi-magnon branch rather than a simple single-particle excitation (Xu et al., 2024).

A further refinement comes from disorder. In the easy-axis triangular antiferromagnet with impurities, the gapless $q=\pi$1-point Goldstone mode remains robust in both low- and high-field supersolids at impurity density $q=\pi$2, whereas comparable disorder splits the lowest magnon bands of the UUD state. This contrast has been proposed as a direct spectroscopic signature of dissipationless spin-superfluid dynamics inside the supersolid (Huang et al., 3 Sep 2025).

The $q=\pi$3-supersolid on the triangular lattice furnishes an instructive counterexample to the one-Goldstone expectation. Because the ordered state breaks $q=\pi$4 down to $q=\pi$5, it supports two Goldstone modes and a symmetry-protected double degeneracy of magnon branches throughout the Brillouin zone (Wang et al., 25 Sep 2025). This indicates that the collective mode structure of spin supersolids is controlled by the precise internal symmetry that is broken, not merely by the coexistence criterion.

5. Realizations across dimensions and platforms

The most quantitatively developed experimental realizations are frustrated triangular-lattice quantum magnets. In $q=\pi$6, the effective model is the easy-axis XXZ antiferromagnet with $q=\pi$7, $q=\pi$8, and $q=\pi$9; the field sequence Y $U(1)$0 UUD $U(1)$1 V is supported by thermodynamics, DMRG/XTRG, neutron diffraction, and spectroscopy, and the two supersolid regimes are further reflected in a giant magnetocaloric effect (Xiang et al., 15 Apr 2025). The neutron propagation vector in the supersolid regimes has in-plane $U(1)$2 character, directly resolving the three-sublattice solid component, while interlayer incommensurability is interpreted as indicative of the superfluid component (Xiang et al., 15 Apr 2025).

A complementary regime is realized in the strongly Ising material $U(1)$3, described by $U(1)$4, $U(1)$5, and $U(1)$6 (Xu et al., 2024). The zero-field state is identified as a Y supersolid by the coexistence of $U(1)$7 longitudinal order and a gapless Goldstone mode, while the field-driven high-field supersolid is predicted to be $U(1)$8-type rather than V-type once correlations are treated beyond mean field (Chen et al., 2024, Xu et al., 2024). The exact $U(1)$9 thermodynamic-limit status of the transverse order parameter is, however, unusually delicate: finite-$S^z$0 scaling suggests that the superfluid order at exactly zero field may become extremely small or even vanish, whereas a tiny field $S^z$1 T stabilizes an unambiguous SSY state (Xu et al., 2024).

Spin supersolidity is not confined to two dimensions. The anisotropic spin-1 Heisenberg chain provides a one-dimensional benchmark in which direct coexistence of finite $S^z$2 and finite $S^z$3 can be established numerically, and the supersolid-to-solid transition displays stiffness critical behavior consistent with $S^z$4 (Rossini et al., 2011). Realistic three-dimensional models also support supersolids. In the layered triangular-lattice model for $S^z$5-AgNiO$S^z$6, the low-field supersolid is continuously connected to the collinear stripe phase and emerges by magnon condensation at a soft mode, with a finite-temperature transition consistent with 3D XY criticality (Seabra et al., 2010). In EuCo$S^z$7Al$S^z$8, the metallic spin supersolid extends the phenomenon from Mott-insulating magnets to a highly conductive metal; transport anomalies, Hall response, and quantum oscillations track the phase structure, while a stacked-triangular RKKY-dipolar model identifies a zero-field Y supersolid and a field-induced V supersolid as genuine three-dimensional states (Xu et al., 25 Mar 2026, Xi et al., 25 Mar 2026).

Related phenomena also appear in ultracold gases, though the terminology is often more cautious. In Raman spin-orbit-coupled Bose gases, the stripe supersolid is the $S^z$9 member of a broader family of supercurrent-carrying supersolid states with quasimomentum-dependent current and instability thresholds (Lyu et al., 2023). In spin-1 and spin-2 spinor condensates with Rashba or Dresselhaus coupling, mean-field calculations find supersolid-like superstripe, superlattice, multi-ring, triangular-lattice, and square-lattice states, typically in trapped or self-bound settings where the crystalline modulation is finite-size rather than a bulk thermodynamic order (Adhikari, 2021, Kaur et al., 2022). These works broaden the vocabulary of spin supersolidity but also reinforce the distinction between bulk magnetic supersolids and supersolid-like patterned condensates.

6. Methods, controversies, and current directions

Progress in the field has depended on methods that can access both symmetry sectors simultaneously. For one-dimensional rings, variational MPS with periodic boundary conditions made it possible to compute stiffness directly and thereby sharpen the definition of the supersolid region (Rossini et al., 2011). For frustrated two-dimensional magnets, infinite PEPS and tangent-space excitation ansätze have become central because they operate directly in the thermodynamic limit and yield momentum-resolved spectra (Xu et al., 2024). Finite- and real-time DMRG on cylinders, including TDVP, has been used to extract both static order and dynamical structure factors in clean and disordered supersolids (Huang et al., 3 Sep 2025). XTRG and related thermal tensor-network methods connect these phases to entropy accumulation and magnetocaloric response (Xiang et al., 15 Apr 2025). Classical Monte Carlo, iDMRG, and flavor spin-wave theory extend the reach to three-dimensional metallic systems and non-XXZ symmetry classes (Xi et al., 25 Mar 2026, Wang et al., 25 Sep 2025).

Two interpretive cautions recur across the literature. First, indirect signatures are insufficient. Earlier work on the spin-1 chain and on triangular magnets often inferred supersolidity from magnetization textures or plateau-adjacent anomalies, but later studies emphasized that genuine identification requires simultaneous evidence for both longitudinal order and stiffness or a Goldstone mode (Rossini et al., 2011, Huang et al., 3 Sep 2025). Second, rotons and continua are not decisive by themselves. Roton minima are prominent in some Y supersolids, but absent in the V supersolid of the easy-axis triangular model; conversely, continuum-like low-energy spectra may result from finite experimental broadening rather than intrinsic fractionalization (Chi et al., 2024).

A further subtlety concerns explicit symmetry breaking. In the SOC triangular-lattice model with symmetry-allowed anisotropic exchanges, infinitesimal spin-orbit coupling destroys the strict $S_i^+$00 supersolid by lifting the accidental $S_i^+$01 manifold and opening a pseudo-Goldstone gap. Yet the resulting $S_i^+$02 anisotropy is irrelevant over a finite temperature window for the Y and $S_i^+$03 states, allowing thermally stabilized KT supersolids, whereas the $S_i^+$04 anisotropy of the V state does not (Park et al., 28 Jan 2026). This is a strong reminder that “supersolid” can denote distinct infrared structures at zero and nonzero temperature once microscopic $S_i^+$05 symmetry is broken.

Open issues are now less about existence in principle than about classification and directness of evidence. The high-field phase in Ising-like cobaltates remains under active discussion, with correlation-driven $S_i^+$06 order replacing the simple V picture in one case (Xu et al., 2024). In $S_i^+$07, the case for supersolidity is strong but the transverse ordered component is still not directly imaged in neutron diffraction, and polarized neutron measurements would sharpen that identification (Xiang et al., 15 Apr 2025). Metallic supersolids raise additional questions about the feedback between itinerant carriers and ordered local moments beyond qualitative Weiss-field descriptions (Xu et al., 25 Mar 2026). More broadly, the field is now moving from establishing coexistence toward understanding universality classes, spin transport, impurity robustness, and materials design, especially in systems where supersolidity is entangled with magnetocaloric functionality or dissipationless spin dynamics (Huang et al., 5 Jan 2026).