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

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.