Topological Phase Transitions and Their Thermodynamic Fate in Arbitrary-$S$ Pyrochlore Spin Ice

TL;DR

This paper develops a theoretical framework to classify topological phases and phase transitions in classical pyrochlore magnets with arbitrary spin, revealing spin-dependent critical behaviors and the effects of monopoles and thermal fluctuations.

Contribution

It introduces a comprehensive classification scheme for topological phases in pyrochlore spin ice with arbitrary spin, including new insights into the nature of phase transitions and monopole effects.

Findings

Integer spins exhibit a 3D XY deconfinement transition.

Half-integer spins remain in a Coulomb liquid without transition.

For S=3/2, the system maps to the 3-state Potts model with a first-order transition.

Abstract

We develop a self-contained theoretical framework that classifies the topological phases and critical phenomena of classical pyrochlore magnets with arbitrary spin $S$, subject to competing exchange and single-ion anisotropies. In the small-$w$ regime, where the single-ion term favors low spin amplitudes, exact dualities reveal a dichotomy: integer spins exhibit a continuous 3D $XY$ deconfinement transition, whereas half-integer spins remain in a $U(1)$ Coulomb liquid without any transition. In the large-$w$ regime, where the local spin amplitudes are maximized ($|S^z| = S$), the macroscopic flux is quantized to multiples of $2S$. By mapping the defect structure to topological loop gases, we prove that the compatibility between the physical ice rule and the emergent $\mathbb{Z}_{2S}$ flux conservation holds if and only if $S \le 3/2$. For $S=3/2$, this maps the system to the 3-state Potts model, whose symmetry-allowed cubic invariant drives a first-order transition. For $S \ge 2$, monopole contamination breaks the discrete clock mapping. Using an exact decomposition of the partition function, we show that the hierarchical string fusion cascade exponentially suppresses the discrete perturbations, which act as a dangerously irrelevant operator at the 3D $XY$ fixed point, protecting 3D $XY$ criticality. Finally, incorporating thermal monopoles, we show that they act as a symmetry-breaking effective magnetic field that severs defect strings. Consequently, the continuous transitions are rounded into crossovers, whereas the first-order $S=3/2$ transition is predicted to survive at finite temperatures, terminating at a critical endpoint. Classical Monte Carlo simulations for $S$ up to $7/2$ corroborate these analytical predictions.