---
title: Spin Supersolid States
url: https://www.emergentmind.com/topics/spin-supersolid-states
type: topic
---

# Spin Supersolid States

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 [2601.01890]. 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 [1102.3545][2402.15869][2603.24445].

## 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 [2601.01890]. 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 [1102.3545].

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,
\[
\langle m_z^2\rangle = \frac{S^z(\mathbf K)}{N'},\qquad
\langle m_\perp^2\rangle = \frac{S^\perp(\mathbf K)}{N'},
\]
with both finite in the supersolid regime [2601.01890].

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 \(U(1)\)-breaking order [2601.01890]. A further generalization appears in the triangular-lattice spin-1 \(U(2)\) model, where the supersolid breaks lattice translation and an \(SU(2)\) spin symmetry, yielding an “\(SU(2)\)-supersolid” with two Goldstone modes rather than one [2509.20772].

## 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,
\[
\mathcal{H} = \sum_j \left[ \frac{J}{2}\left(S_j^+ S_{j+1}^- + \mathrm{h.c.}\right) + \Delta S_j^z S_{j+1}^z \right] + D\sum_j (S_j^z)^2 - B\sum_j S_j^z ,
\]
with \(J=1\) and \(D=\Delta/2\), the XY exchange delocalizes spin excitations while the Ising and single-ion anisotropies favor staggered longitudinal order; in the large-\(\Delta\) regime the model maps to an effective spin-\(\tfrac12\) XX chain in a field, making the supersolid interpretable as mobile quasiparticles propagating on an ordered background [1102.3545].

In two-dimensional frustrated magnets the canonical minimal model is the easy-axis triangular-lattice XXZ antiferromagnet,
\[
\hat{H}= \sum_{\langle ij\rangle} \big [J_{zz} S^z_i S^z_j + J_{xy}\big(S^x_i S^x_j + S^y_i S^y_j \big)\big] - g\mu_B B_z \sum_{i} S_i^z,
\]
or equivalent notations with \(J_z\) and \(J_{xy}\) [2405.05151][2404.14163]. Frustration stabilizes three-sublattice longitudinal order, easy-axis anisotropy favors the solid component, and field tuning drives the familiar sequence
\[
\text{Y supersolid} \to \text{UUD} \to \text{V or }\Psi \text{ state} \to \text{polarized},
\]
with the supersolids flanking the UUD plateau in the weak-easy-axis case [2404.14163][2504.11298]. In the strongly Ising regime of \(\mathrm{K}_2\mathrm{Co}(\mathrm{SeO}_3)_2\), correlations stabilize a high-field \(\Psi\) (\(\pi\)-coplanar) state rather than the V state of simple cluster mean-field theory [2405.05151].

Three-dimensional routes are now also established. In the realistic model for \(2H\)-AgNiO\(_2\), a low-field supersolid emerges continuously from a collinear stripe state by softening of a transverse magnon mode at \(M'\), i.e. by magnon Bose condensation on top of a crystalline background [1003.3430]. In EuCo\(_2\)Al\(_9\), 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 \(\to\) UUD \(\to\) V \(\to\) polarized [2603.24446].

The same conceptual structure appears in mixed-spin models and cold-atom analogues. Coupled alternating spin-1 and spin-\(\tfrac12\) 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 [1712.06270]. 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 [2311.13786][2104.02197][1909.11871][2201.07500].

## 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^{zz}(q) = \frac{1}{N}\sum_{j,\ell} e^{-iq(j-\ell)} \langle S_j^z S_\ell^z\rangle
\]
defines the solid order parameter
\[
\mathcal{O}_{SDW} = \lim_{N\to\infty} \frac{S^{zz}(\pi)}{N},
\]
while the off-diagonal sector is quantified by the spin stiffness
\[
\rho_s = N \frac{\partial^2 E_0(\phi)}{\partial \phi^2}\Big|_{\phi=0},
\]
obtained from the curvature of the twisted-boundary ground-state energy [1102.3545]. The phase classification is then exact in the operational sense: solid if \(\mathcal{O}_{SDW}\neq0\) and \(\rho_s=0\), superfluid if \(\mathcal{O}_{SDW}=0\) and \(\rho_s\neq0\), and supersolid if both are nonzero [1102.3545].

On the triangular lattice, the same structure is used in momentum space at \(K\). In DMRG studies with impurities, the diagonal and transverse channels are measured as
\[
\langle m_z^2\rangle = \frac{S^z(K)}{L_y^2},\qquad
\langle m_\perp^2\rangle = \frac{S^{xy}(K)}{L_y^2},
\]
so that the supersolid is the regime where both remain finite [2509.03489]. Tensor-network work on \(\mathrm{K}_2\mathrm{Co}(\mathrm{SeO}_3)_2\) instead uses explicit three-sublattice order parameters,
\[
m^{st}_{xy}=\left|\frac{1}{3}\sum_{n\in \mathrm{u.c.}}e^{2in\pi/3}\big(\langle S_n^x\rangle \hat{x}+\langle S_n^{y}\rangle\hat{y}\big)\right|,
\]
\[
m_z^{st}=\left|\frac{1}{3}\sum_{n\in \mathrm{u.c.}}e^{2in\pi/3}\langle S_n^z\rangle\right|,
\]
with supersolidity identified by \(m^{st}_{xy}\neq0\) and \(m_z^{st}\neq0\) [2405.05151].

In \(\mathrm{Na}_2\mathrm{BaCo}(\mathrm{PO}_4)_2\), the solid order is encoded in a complex three-sublattice order parameter,
\[
\Psi \equiv \frac{1}{N}\left( \sum_{i\in A}\langle S_i^z\rangle + \sum_{j\in B}\langle S_j^z\rangle e^{i2\pi/3} + \sum_{k\in C}\langle S_k^z\rangle e^{i4\pi/3} \right),
\]
while the superfluid component is the coherent transverse order \(\Phi \sim \frac{1}{N}\sum_i\langle S_i^+\rangle\) [2404.15997]. In the spin-1 \(U(2)\) model, the solid component is carried by three-sublattice order in \(\lambda^8\), whereas the superfluid component is an ordered \(\boldsymbol{\lambda}^{\mathrm{su2}}=(\lambda^4,\lambda^5,\lambda^3)\) sector; the supersolid is therefore diagnosed by simultaneous \(K\)-point structure in both channels [2509.20772].

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 \(\rho_s=0\) 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 [1102.3545].

## 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 \(K\) points; this is seen in tensor-network spectra for \(\mathrm{K}_2\mathrm{Co}(\mathrm{SeO}_3)_2\), in DMRG/TDVP spectra for the Y and V supersolids of the easy-axis triangular antiferromagnet, and in neutron and tensor-network studies of \(\mathrm{Na}_2\mathrm{BaCo}(\mathrm{PO}_4)_2\) [2405.05151][2509.03489][2404.15997]. The presence of this mode distinguishes supersolids sharply from the UUD plateau, whose transverse excitations are ordinary gapped magnons [2509.03489].

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 \(M\), but under field the roton associated with the Goldstone branch disappears quickly, whereas the roton in the gapped branch persists across the Y phase [2404.14163]. 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 [2404.14163].

In \(\mathrm{Na}_2\mathrm{BaCo}(\mathrm{PO}_4)_2\), the supersolid dynamics are even richer. The low-energy spectrum exhibits a “double magnon-roton” structure: one branch contains the true \(U(1)\) Goldstone mode and a roton near \(M\), 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 [2404.15997]. The same work connects these low-energy roton branches to the large sub-Kelvin entropy and giant magnetocaloric response of the material [2404.15997].

For \(\mathrm{K}_2\mathrm{Co}(\mathrm{SeO}_3)_2\), 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,
\[
\omega_n = nJ_{zz},
\]
with visible features at \(\omega_1\simeq 3~\mathrm{meV}\), \(\omega_2\simeq 6~\mathrm{meV}\), and \(\omega_3\simeq 9~\mathrm{meV}\). Real-space inspection of the excited tensors shows that the \(\omega_2\) sector has nonlocal two-spin-flip character, supporting its interpretation as a multi-magnon branch rather than a simple single-particle excitation [2405.05151].

A further refinement comes from disorder. In the easy-axis triangular antiferromagnet with impurities, the gapless \(K\)-point Goldstone mode remains robust in both low- and high-field supersolids at impurity density \(\sim 1.85\%\), 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 [2509.03489].

The \(SU(2)\)-supersolid on the triangular lattice furnishes an instructive counterexample to the one-Goldstone expectation. Because the ordered state breaks \(SU(2)\) down to \(U(1)\), it supports two Goldstone modes and a symmetry-protected double degeneracy of magnon branches throughout the Brillouin zone [2509.20772]. 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 \(\mathrm{Na}_2\mathrm{BaCo}(\mathrm{PO}_4)_2\), the effective model is the easy-axis XXZ antiferromagnet with \(J_{xy}=0.88~\mathrm{K}\), \(J_z=1.48~\mathrm{K}\), and \(g_c\simeq4.89\); the field sequence Y \(\to\) UUD \(\to\) V is supported by thermodynamics, DMRG/XTRG, neutron diffraction, and spectroscopy, and the two supersolid regimes are further reflected in a giant magnetocaloric effect [2504.11298]. The neutron propagation vector in the supersolid regimes has in-plane \((1/3,1/3)\) character, directly resolving the three-sublattice solid component, while interlayer incommensurability is interpreted as indicative of the superfluid component [2504.11298].

A complementary regime is realized in the strongly Ising material \(\mathrm{K}_2\mathrm{Co}(\mathrm{SeO}_3)_2\), described by \(J_{zz}=2.98~\mathrm{meV}\), \(J_{xy}=0.21~\mathrm{meV}\), and \(g=7.8\) [2405.05151]. The zero-field state is identified as a Y supersolid by the coexistence of \(\sqrt{3}\times\sqrt{3}\) longitudinal order and a gapless Goldstone mode, while the field-driven high-field supersolid is predicted to be \(\Psi\)-type rather than V-type once correlations are treated beyond mean field [2402.15869][2405.05151]. The exact \(B_z=0\) thermodynamic-limit status of the transverse order parameter is, however, unusually delicate: finite-\(D\) scaling suggests that the superfluid order at exactly zero field may become extremely small or even vanish, whereas a tiny field \(B_z=0.5\) T stabilizes an unambiguous SSY state [2405.05151].

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 \(\mathcal O_{SDW}\) and finite \(\rho_s\) can be established numerically, and the supersolid-to-solid transition displays stiffness critical behavior consistent with \(\beta_s=\tfrac12\) [1102.3545]. Realistic three-dimensional models also support supersolids. In the layered triangular-lattice model for \(2H\)-AgNiO\(_2\), 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 [1003.3430]. In EuCo\(_2\)Al\(_9\), 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 [2603.24445][2603.24446].

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 \(k=0\) member of a broader family of supercurrent-carrying supersolid states with quasimomentum-dependent current and instability thresholds [2311.13786]. 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 [2104.02197][2201.07500]. 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 [1102.3545]. 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 [2405.05151]. 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 [2509.03489]. XTRG and related thermal tensor-network methods connect these phases to entropy accumulation and magnetocaloric response [2504.11298]. Classical Monte Carlo, iDMRG, and flavor spin-wave theory extend the reach to three-dimensional metallic systems and non-XXZ symmetry classes [2603.24446][2509.20772].

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 [1102.3545][2509.03489]. 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 [2404.14163].

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 \(T=0\) supersolid by lifting the accidental \(U(1)\) manifold and opening a pseudo-Goldstone gap. Yet the resulting \(Z_6\) anisotropy is irrelevant over a finite temperature window for the Y and \(\Psi\) states, allowing thermally stabilized KT supersolids, whereas the \(Z_3\) anisotropy of the V state does not [2601.20963]. This is a strong reminder that “supersolid” can denote distinct infrared structures at zero and nonzero temperature once microscopic \(U(1)\) 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 \(\Psi\) order replacing the simple V picture in one case [2405.05151]. In \(\mathrm{Na}_2\mathrm{BaCo}(\mathrm{PO}_4)_2\), 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 [2504.11298]. Metallic supersolids raise additional questions about the feedback between itinerant carriers and ordered local moments beyond qualitative Weiss-field descriptions [2603.24445]. 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 [2601.01890].

Source: https://www.emergentmind.com/topics/spin-supersolid-states