Classification of Classical Spin Liquids: Typology and Resulting Landscape

Classical spin liquids (CSL) lack long-range magnetic order and are characterized by an extensive ground state degeneracy. We propose a classification scheme of CSLs based on the structure of the flat bands of their Hamiltonians. Depending on absence or presence of the gap from the flat band, the CSL are classified as algebraic or fragile topological, respectively. Each category is further classified: the algebraic case by the nature of the emergent Gauss's law at the gap-closing point(s), and the fragile topological case by the homotopy of the eigenvector winding around the Brillouin zone. Previously identified instances of CSLs fit snugly into our scheme, which finds a landscape where algebraic CSLs are located at transitions between \fragile topological ones. It also allows us to present a new, simple family of models illustrating that landscape, which hosts both fragile topological and algebraic CSLs, as well as transitions between them.


Introduction.
Research into magnets without longrange order has a long history, from the study of the effects of disorder in spin glasses [1,2], to the proposal of resonating valence bond states [3,4] that underpins much of modern research in strongly frustrated magnets.Classical spin liquids (CSLs) can arise when the ground state manifold of a continuous spin model is extensively degenerate, representing the extreme limit of the consequences of frustration, when fluctuations between ground states preclude any form of order altogether [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19].Even though CSLs tend to be unstable to perturbations at T = 0, their large entropy at low energies can allow them to dominate the surrounding phase diagram at finite T .This makes them extremely relevant to the finite temperature physics of real frustrated systems.In addition, they can often be usefully thought of as parent states, or intermediate temperature limits, of quantum spin liquids which arise when quantum fluctuations introduce dynamics between the classical ground states [20][21][22][23][24][25][26][27][28][29].
It is therefore important to understand and classify the spin liquids.While classification schemes for QSLs have been successful, notably using the projective symmetry group [30], and the modern perspective of gapped QSLs corresponding to a topological quantum field theory has led to further efforts in classifying QSLs as symmetryenriched topological phases [31,32], no similarly comprehensive classification exists for the CSLs.Previous works have classified frustrated classical spin systems using constraint counting [7], linearization around given spin configurations [33], supersymmetric connections between models [34] or using topological invariants built for specific cases [17].A scheme which generalizes across different models and types of CSL, and which depends on the physics of the CSL as a whole rather than individual spin configurations within it remains to be found.Concretely, the question whether CSLs have hidden topological properties like QSLs is an open one.Further, it is not clear how to place spin liquids with algebraically-decaying and exponentially-decaying correlations in a single scheme.
We present a classification scheme to address these issues, providing a concrete and relatively simple framework to distinguish different CSL states.Our work also offers practical tools to diagnose known CSLs, as well as to construct new ones with desired properties.We demonstrate in particular that certain CSLs can be characterised by a topological invariant that persists as long as the lowest flat bands of the spectrum of the Hamiltonian are separated by a gap from the higher dispersive band(s).We term these "fragile topological" classical spin liquids (FT-CSLs), as adding additional spins to the unit cell can render them topologically trivial without closing the spectral gap.
Algebraic spin liquids can be viewed as inhabiting the arXiv:2305.00155v3[cond-mat.str-el] 1 Jun 2023 boundaries between the FT-CSLs where the spectral gap closes, illustrated schematically in Fig. 1(a).This work presents a largely concept-based, non-technical account of the central narrative underpinning the classification scheme.We will present an extended discussion, including a detailed description alongside the more technical aspects, together with a broader set of new CSLs obtained within this scheme, in a long companion paper [35].
It turns out that an adequate description of such models can often be attained in a self-consistent Gaussian (also known as soft-spin, or large-N ) approximation [8,11], which is inspired by the Luttinger-Tisza method [37][38][39].This amounts to abandoning the hard spin-length constraint |S i | = 1, and replacing it with an average constraint ⟨S 2 ⟩ = 1.This results in a solvable Hamiltonian bilinear in spin variables, which can be diagonalised in momentum space, yielding a spectrum with a band structure on which our classification is based [40].
The extensive ground state degeneracy of the CSL is reflected in the existence of (at least one) flat band at the bottom of the spectrum in momentum space -because all the plane-wave states with a given momentum q in this band are degenerate.The central concept behind our classification scheme is that the unique physics of different CSLs is encoded in their corresponding band structures.The most basic distinction is whether or not the bottom flat bands are separated by a gap from the higher energy ones, illustrated in Fig. 1.Spectra with a band gap have short-range classical spin correlations.The absence of the gap, on the other hand, results in an algebraic CSL, of which the aforementioned U(1) Coulomb phase is one example.We shall formulate the classification of both these types of CSLs, starting with the latter.
Algebraic CSLs.We begin by examining algebraic CSLs.Due to the presence of gap-closing points, these systems exhibit spin correlations that decay algebraically.The well-known examples of Heisenberg antiferromagnet on the pyrochlore lattice [6-10, 23, 36, 41-44], and also various other models [11,12,[14][15][16][17][18][19], belong to this category.While the gap-closing manifold can be more general, such as lines or surfaces [14], here we limit our discussion to gap-closing points only.By examining the eigenvector configuration around these points, we demonstrate how to extract a generalized Gauss's law, which plays a crucial role in describing the low-energy effective field theory of the CSL.The key steps of this procedure are outlined next.We consider a system with N spins per unit cell.The Hamiltonian for CSL systems can generally be expressed as a constrainer Hamiltonian: (1) C(R), which we call a constrainer, involves the sum of spins with different coefficients in a local region around the unit cell located at R (see Eq. ( 11) for an example).The Hamiltonian is the translationally invariant sum of such squared constrainers.The ground states of the system are spin configurations such that all constrainers are zero.The constrainer formalism is quite a general form of writing Heisenberg Hamiltonians that realize CSLs.Unless stated otherwise, the following discussion focuses on two spatial dimensions.
For simplicity, we assume there to be one constrainer per unit cell, but its generalization is not difficult.The spectrum of the Hamiltonian then contains N −1 bottom flat bands and one top dispersive band with a gap-closing point which we set to be at q = q 0 .The eigenvector T(q) of the top band, computed as the Fourier transform of the constrainer (see example later), enables us to write down the N ×N bilinear interaction matrix in momentum space explicitly featuring N −1 bottom flat bands at ω = 0 and a top band with a dispersion relation ω Let us analyze the behavior of this eigenvector for small wave-vectors k = q − q 0 .The gap-closing condition implies that T(q 0 ) = 0. Hence we can express its components, denoted as T a (q 0 + k), as Taylor expansions in k x and k y without the zeroth order constants.The leading terms in this polynomial expansion are where m a is the degree of the leading order terms in T a (q) and the constants c aj are given by the Taylor expansion.Ground states are associated with the (N−1)-dimensional space of (complex-valued) eigenvectors S orthogonal to T(q 0 + k), which satisfy Performing the inverse Fourier transformation into real space yields the generalized Gauss's law: and the ground state condition is ρ = 0.
Crucially, this analysis also yields the equal-time spin structure factor, the intensity distribution of which, 1 − | a T a | 2 /T 2 , exhibits singular patterns at q 0 , known as pinch points [8,9,[45][46][47][Fig.1(b)].Such Gauss's laws play a central role in describing the properties of the ground state manifold of algebraic CSLs.Specifically, the long-wavelength expansion results in an effective Hamiltonian given by H eff = ( N a=1 D (ma) a S a ) 2 , which properly captures the algebraic spin correlation of the system.The generalized Gauss's laws can also give rise to non-trivial physics, such as multipole conservations and fracton charges, which have garnered attention in various fields of physics [48][49][50][51][52]. Equipping this Hamiltonian with quantum dynamics provides the starting point for building the emergent (generalized) electrodynamics describing the corresponding QSL.In addition, the presence of multiple gap-closing points allows for the coexistence of different generalized Gauss's laws that describe the same ground state manifold, depending on the background wave-vectors q in the long-wavelength limit.
We can now distinguish different algebraic CSLs by comparing their gap-closing points.Concretely, two algebraic CSLs can be viewed as in the same class if one can adiabatically deform the constrainer Hamiltonian and turn the Gauss's law of one CSL into that of the other, without going through singular processes of merging/splitting/lifting some of these points.On the other hand, two algebraic CSLs are categorically different, if they have a different number of gap-closing points, or the associated Gauss's laws involve a different number of effective electric field degrees of freedom, or a different order of ∂ x , ∂ y .These gap-closing points cannot be made identical without going through certain singular transitions.Below, we provide an example to demonstrate how the Gauss's law is extracted from a concrete model, and how the merging/splitting of the gap-closing point happens.
Fragile Topological CSLs.We now turn to the second category of CSLs, FT-CSLs.They are characterized by bottom flat bands that are completely separated from other bands above them by a gap, resulting in exponentially decaying spin correlations [13].
The FT-CSLs can be further classified based on the homotopy class of the bottom band eigenvector configuration, which can only change when the system undergoes a gap-closing topological phase transition.Hence the phase boundaries of such FT-CSLs are inhabited by the algebraic CSLs.This is similar to the concept of topological band transitions, where the Chern number of a band cannot change without closing a gap.
The specific classification scheme for 2D FT-CSLs with N sub-lattice sites and one (or, analogously, N − 1, see below) bottom flat bands works as follows.Consider a normalized single bottom flat band eigenvector, denoted by B(q) (or the top band eigenvector T(q) if there is only one top band).The components of B(q) are generally complex, but for certain symmetries they can be also real.The flatness of the band ensures that it has a zero Chern number [53], which means the eigenvector B(q) is smoothly defined over the entire Brillouin zone (BZ).Since the BZ is a 2-torus, B(q) defines a map from the torus to the target space of CP N −1 or RP N −1 : B(q) : T 2 → CP N −1 (or RP N −1 ); q → B(q).( This eigenvector can still "wind" nontrivially over the BZ -captured by the homotopy class of the corresponding map [T 2 , CP N −1 ] or [T 2 , RP N −1 ].The homotopy class can only change at gap-closing points, where B(q) becomes ill-defined.
In the case where B(q) is complex, we can use the fact that the complex projective space is simply connected (π 1 (CP N −1 ) = 0) for any N − 1 ≥ 1 to obtain: For real-valued B(q), the homotopy class of [T 2 , RP N −1 ] does not have a simple formula.However, in a special case when one can consistently assign directions to the RP N −1 eigenvectors over the boundary condition of the BZ, the homotopy group simplifies: The only non-trivial case is when N = 3, which is the skyrmion number on the torus, Q sk , given by In the case of three-dimension models, we need to com- Note that the homotopy class is a fragile topological quantity in the following sense: if new spins are added to each unit cell and interact with the original spins, the non-trivial homotopy class may become trivial in the new model.By padding each unit cell with auxiliary spins, adiabatically tuning the CSL Hamiltonian and then decoupling them, one can change the homotopy class without closing the spectral gap.Example.We illustrate the classification scheme with the concrete example of a kagome lattice [Fig.2(a)] model.The constrainer Hamiltonian reads: Here, the sites 1, . . ., 6 and 1 ′ , . . ., 6 ′ in C KGM (R) are labeled in Fig. 2(b) for the hexagonal star located at R. The case of ξ 1 = ξ 2 = 0 was studied in Ref. [13].Our model has one constrainer per unit cell and three sub-lattice sites, leading to a spectrum that consists of two degenerate bottom flat bands and a top dispersive band.The eigenvector of the top band can be expressed as the Fourier transform of C KGM (R): and its dispersion is ω(q) = |T(q)| 2 .
We compute Q sk , Eq.( 9), to determine the homotopy class of the eigenvector T(q) of the top dispersive band (instead of the bottom band eigenvector B(q), which amounts to replacing B → T in the above formalism).Tuning the two parameters ξ 1,2 yields the diverse phases shown in Fig. 2(c), labeled by their skyrmion numbers.
Decreasing the value of ξ 2 moves along the phase boundary, but the spectrum of the Hamiltonian and the emergent Gauss's law changes, with the 4-fold pinch point (4FPP) splitting into two 2-fold pinch points [Fig.3(a,b)]: this transition between algebraic CSLs involves the merging and splitting of gap-closing points.
By contrast, upon increasing ξ 2 > −1, the gap-closing point is lifted, and a gap opens up between the flat and dispersive bands.This yields a FT-CSL, as shown in Fig. 3(e,f).Other phase boundaries of algebraic CSLs are also interesting, but we will not delve into them extensively.We do note that with our classification methodology, analysing these phase boundaries is now a straightforward, basic algebraic calculation.
Summary.In this work, we have presented a classification scheme for classical spin liquids, which we divide into two broad categories: algebraic and fragile topological CSLs.In our extended companion paper [35], we present a comprehensive analysis of the classification scheme, including aspects omitted here such as higher-dimensional CSLs, more complex band-closing structures, and connection to flat band theories.Additionally, we construct a variety of new models using the constrainer Hamiltonian formalism to illustrate the different aspects of the classification scheme, and the new physics arising from it.

FIG. 1 .
FIG. 1.(a) The landscape of CSLs consists of Fragile Topological CSLs (FT-CSLs) whose boundaries are Algebraic CSLs.(b) Algebraic CSLs feature gap closings between the bottom flat bands and higher dispersive bands of the spectrum of the exchange Hamiltonian.The band touching point determines the emergent Gauss's law.(c) FT-CSLs have no such gapclosing points, being classified by their eigenvector homopoty.

FIG. 2 .
FIG. 2. The kagome model of HKGM in Eq. (10) (a) The Kagome lattice with different sublattice sites colored differently.(b) The sites for the constrainer CKGM defined in Eq. (11).(c) The phase diagram of the model Eq.(10).Different colored regions are phases of fragile topological CSL with different skyrmion numbers Q.The three pink stars are the parameter sets shown and Fig.3.

FIG. 3 .
FIG. 3. Spin structure factor and spectrum of three parameter sets (pink stars in Fig. 2) of the model Eq.(10), highlighting the merging and lifting of the gap-closing point (blue circle), which indicates transition between algebraic CSLs and fragile topological CSLs.(a) Spin Structure for parameter ξ1 = −1, ξ2 = −1.3.(b) The two gap-closing points in its spectrum each hosts a 2-fold pinch point (2FPP).(c) Spin Structure for parameter ξ1 = −1, ξ2 = −1.(d) The previous two gap-closing points merge and form a single gapclosing point with 4FPP.(e) Spin Structure for parameter ξ1 = −1, ξ2 = −0.7.(f) The gap-closing point lifts up and opens up the gap.