N-state Potts ices as generalizations of classical and quantum spin ice

Classical and quantum spin ice models are amongst the most popular settings for the study of spin liquid physics. $N-$state Potts ice models have been constructed that generalize spin ice, hosting multiple emergent $\text{U}(1)$ gauge fields and excitations charged under non-trivial combinations of these fields. We present a general treatment of classical $N-$state Potts ices relating their properties to the $\mathfrak{su}(N)$ Lie algebras, and demonstrate how the properties of charged excitations in the classical model can be related to this symmetry group. We also introduce quantum generalizations of the Potts Ice models, and demonstrate how charge flavor changing interactions unique to $N>2$ models dominate their low energy physics. We further show how symmetries inherited from the $\mathfrak{su}(N)$ can lead to flux vacuum frustration, greatly modifying the dynamical properties of charged excitations.

💡 Research Summary

This paper presents a comprehensive study of N‑state Potts ice models as a broad generalization of both classical and quantum spin ice. The authors first construct a classical Potts ice on a bipartite lattice (or its line‑graph) by assigning one of N_c colors to each edge (or vertex) and imposing an “equal‑color‑number” rule: at every vertex exactly z copies of each color must meet, where the coordination number is N = z N_c. This rule is the natural extension of the two‑in‑two‑out ice rule of spin ice. By introducing N_F‑component electric fields E_i^c for each color, the authors rewrite the Hamiltonian as a sum of squared divergences Q_i^C = η_C ∑_{r∈C} E_i^{c_r}. The ground‑state manifold corresponds to Q_i^C = 0, i.e., a divergence‑free condition, while violations create charged excitations.

A key insight is that the set of diagonal matrices built from the vectors E_i^c spans the Cartan sub‑algebra of the Lie algebra su(N_c). The non‑diagonal generators R_α, labelled by the root vectors α, act as color‑swap operators. Acting with R_α on a ground‑state bond flips a color and creates a pair of defects carrying charges ±α. Thus every elementary excitation (“roton”) carries a charge equal to a root of su(N_c). The emergent gauge structure is therefore (N_c − 1) independent U(1) gauge fields generated by the Cartan sub‑algebra, i.e., an abelian projection of an SU(N_c) gauge theory.

The entropic interaction between two rotons is proportional to the Killing‑form inner product (α,β). Because for su(N_c) this product can only take the values 2, 1, or 0 (after appropriate normalization), the Coulomb‑like interaction is strongest for opposite charges of the same color pair, weaker when the two rotons share only one color, and vanishes when they share none. Monte‑Carlo simulations on several lattices—square (N_c = 4, 2D), pyrochlore (N_c = 4, 3D), and the “octachlore” (line‑graph of the cubic lattice, N_c = 6, 3D)—confirm the presence of a Coulomb phase. Worm‑length statistics display power‑law distributions identical to those of conventional spin ice, with exponents ≈ −2 in two dimensions and a crossover from ≈ −2.5 to ≈ −1.0 in three dimensions, indicating random‑walk behavior of the emergent electric field lines.

The authors then turn to quantum extensions. By adding kinetic terms for the electric fields they generate (N_c − 1) gapless photon modes and gapped vison excitations, mirroring the quantum spin‑ice picture. However, for N_c ≥ 3 a new class of three‑field “flavor‑changing” interactions appears, originating from the off‑diagonal su(N_c) generators. In a gauge‑mean‑field treatment these terms dominate the low‑energy physics, allowing rotons to change their charge flavor without emitting photons. This mechanism has no analogue in the N_c = 2 (Ising) case.

Furthermore, the symmetry of the Cartan sub‑algebra imposes constraints on the allowed flux configurations in the ground state. For N_c ≥ 3 the system cannot simultaneously satisfy all flux‑free conditions, leading to “flux‑vacuum frustration”. This frustration modifies the dynamics of both photons and rotons, producing anisotropic photon dispersion and restricting roton motion to sub‑manifolds of the lattice.

To illustrate phenomena absent in ordinary quantum spin ice, the authors construct exactly solvable toric‑code‑like models that incorporate the flavor‑changing terms and flux frustration. These models display topologically protected degeneracies and non‑trivial braiding statistics of rotons, highlighting the richness of the N‑state Potts ice landscape.

In summary, the work establishes a deep connection between Potts ice models and the algebraic structure of su(N). Classical Potts ice realizes a multi‑U(1) Coulomb phase with root‑charged excitations, while its quantum counterpart introduces novel flavor‑changing dynamics and frustrated flux vacua. The results open pathways for realizing higher‑symmetry gauge theories in experimental platforms such as ultracold atoms, Rydberg arrays, and superconducting qubit lattices, and suggest many future directions including non‑equilibrium dynamics, disorder effects, and extensions to higher dimensions.