Dualities and Topological Classification of the $S=1$ Pyrochlore Spin Ice

Abstract: We resolve the phase diagram of the $S=1$ pyrochlore spin ice, which exhibits trivial paramagnetic, U(1) Coulomb, and spin nematic phases. In the monopole-free limit, the system can be effectively mapped onto 3D $XY$ and Ising loop-gas models depending on the spin anisotropy, which provides theoretical estimates for the phase boundaries, while a macroscopic flux vector classifies the topological sectors via geometric parity rules. At finite temperatures, thermal monopoles act as a symmetry-breaking field in the continuous $XY$ wave picture and topologically sever defect strings in the loop-gas picture, rounding the phase transitions into continuous crossovers. These theoretical findings are corroborated by classical Monte Carlo simulations.

• The paper establishes duality mappings that identify 3D XY criticality at the low-fugacity transition and 3D Ising criticality at the high-fugacity transition, with predicted critical fugacities of approximately 0.394 and 1.884.
• The paper classifies monopole-free phases through the parity of a quantized global flux: zero flux in the trivial paramagnet, integer sectors in the U(1) Coulomb phase, and even-integer sectors in the spin nematic.
• The paper shows that finite-temperature monopoles act as dual symmetry-breaking fields and string endpoints, eliminating sharp transitions and producing crossovers confirmed by nonsingular specific-heat peaks and continuously leaking flux quantization.

The $S=1$ pyrochlore spin ice of Blume-Capel type extends the classical Ising spin ice by allowing spins to occupy a nonmagnetic $S^z=0$ state, whose statistical weight is controlled by a fugacity $w = e^{-\mu/T}$. Recent Monte Carlo work by Pandey and Damle identified three regimes in this model—a trivial paramagnet, a U(1) Coulomb phase, and a spin nematic with confined $\mathbb{Z}_2$ flux—and conjectured 3D $XY$ and 3D Ising criticality at the two boundaries (Pandey et al., 12 Dec 2025). The paper under review supplies the missing theoretical foundation: explicit duality mappings onto standard lattice models, an exact parity-based topological classification via a macroscopic flux vector, and a demonstration that thermal monopoles destroy both transitions at any finite temperature (2603.03852).

Model and monopole-free limit

The model lives on the pyrochlore lattice, viewed as links $\ell$ of the diamond lattice, with Hamiltonian

$$H = \frac{J}{2}\sum_{\bm{r}} Q_{\bm{r}}^2 + \mu \sum_\ell (S_\ell^z)^2,$$

where $Q_{\bm{r}} = \nabla \cdot \bm{S}_{\bm{r}}$ is the discrete divergence (monopole charge) at each diamond site. Writing fugacities $w = e^{-\mu/T}$ for $S^z=0$ states and $S^z=0$0 for monopoles, the partition function becomes a constrained sum over charge configurations. In the low-temperature limit $S^z=0$1, monopoles are suppressed ($S^z=0$2) and the ice rule $S^z=0$3 holds exactly; all analytic results below pertain to this monopole-free manifold.

Duality mapping to the 3D XY model

For small $S^z=0$4, enforcing Gauss's law with a continuous phase field $S^z=0$5 on the diamond lattice and performing a discrete integration by parts factorizes the partition function over links into $S^z=0$6. For small $S^z=0$7, expanding the logarithm yields precisely the classical 3D $S^z=0$8 model with ferromagnetic coupling $S^z=0$9. Because the neglected higher harmonics are irrelevant operators in the renormalization-group sense, this mapping rigorously establishes that the transition at $w = e^{-\mu/T}$0 belongs to the 3D $w = e^{-\mu/T}$1 universality class—converting the earlier numerical conjecture into a theoretical result. Using $w = e^{-\mu/T}$2, the authors estimate $w = e^{-\mu/T}$3, in reasonable agreement with the numerical value $w = e^{-\mu/T}$4; the residual discrepancy is attributed to higher-order terms in the generalized $w = e^{-\mu/T}$5 action.

Loop-gas representation and the 3D Ising transition

In the large-$w = e^{-\mu/T}$6 regime, $w = e^{-\mu/T}$7 states become minority defects that, by the ice rule, form closed strings. Pauling counting gives a loop weight $w = e^{-\mu/T}$8 per unit length, so that

$w = e^{-\mu/T}$9

which is isomorphic to the high-temperature expansion of the zero-field 3D Ising model with $\mathbb{Z}_2$0. Loop proliferation—the onset of the U(1) Coulomb phase—therefore maps onto the Ising ferromagnetic transition, establishing the 3D Ising nature of $\mathbb{Z}_2$1. Equating $\mathbb{Z}_2$2 with $\mathbb{Z}_2$3 predicts $\mathbb{Z}_2$4, remarkably close to the numerical value $\mathbb{Z}_2$5. The authors are careful to note that Pauling's approximation neglects short-range correlations such as loop self-intersections, but argue these constitute irrelevant local perturbations that renormalize the effective fugacity without changing the universality class—an assumption rather than a rigorous proof.

Topological classification by global flux parity

In the monopole-free limit, the macroscopic flux $\mathbb{Z}_2$6 is quantized: $\mathbb{Z}_2$7. A geometric congruence argument shows that the parity of each component equals the parity $\mathbb{Z}_2$8 of the number of $\mathbb{Z}_2$9 links crossing a horizontal plane. This yields a clean classification:

Phase	Flux sectors	Mechanism
Trivial paramagnet ($XY$0)	$XY$1	Directed loops cannot wrap
U(1) Coulomb ($XY$2)	$XY$3	Odd-crossing proliferating strings unlock parity
Spin nematic ($XY$4)	$XY$5	Even-crossing constraint confines odd flux

This parallels the exact macroscopic-loop-parity description of the 3D toric code, but with a crucial difference: here the classification relies entirely on the absence of monopoles.

Finite-temperature rounding of transitions

At any $XY$6, the monopole fugacity $XY$7 is finite, and the topological distinctions dissolve. Repeating the $XY$8 duality with monopoles included produces a sine-Gordon term $XY$9—a uniform field on the dual $\ell$0 model that explicitly breaks the emergent continuous symmetry. In the loop-gas picture, monopoles act as string endpoints, generalizing the sum from closed loops to arbitrary graphs with endpoint weight $\ell$1, which matches the high-temperature expansion of the 3D Ising model in a magnetic field. Since a uniform field eliminates the Ising transition, both phase boundaries round into continuous crossovers. The Supplemental Material derives these mappings term by term, including the local state sums $\ell$2 for sites of defect degree $\ell$3. This fragility contrasts sharply with 3D $\ell$4 topological phases such as the toric code and Kitaev spin liquids, where strict local constraints forbid string endpoints and preserve sharp finite-$\ell$5 transitions.

Monte Carlo verification

Classical Monte Carlo simulations (heat-bath updates supplemented by a short-loop algorithm in the frustrated regime, with up to $\ell$6 thermalization sweeps) confirm the theory. The specific heat exhibits broad maxima whose height saturates with system size—no thermodynamic divergence—and peak positions drift away from the monopole-free phase boundaries as temperature increases. Most strikingly, the measured flux components deviate continuously from their quantized values, with the deviation measure $\ell$7 tracking the finite monopole density $\ell$8: the fractional leakage of winding numbers directly visualizes the loss of topological quantization. One caveat noted by the authors is that $\ell$9 is size-dependent, approaching $$H = \frac{J}{2}\sum_{\bm{r}} Q_{\bm{r}}^2 + \mu \sum_\ell (S_\ell^z)^2,$$0 in the random limit, so it serves as a diagnostic rather than a scaling observable.

Limitations and open questions

Several assumptions qualify the results. Both duality estimates of the critical fugacities rest on controlled approximations (truncation of the generalized $$H = \frac{J}{2}\sum_{\bm{r}} Q_{\bm{r}}^2 + \mu \sum_\ell (S_\ell^z)^2,$$1 action; Pauling counting), which account for the systematic deviations from the numerically exact values of Pandey and Damle. The claim that geometric corrections to the loop weight are irrelevant perturbations is argued, not proven. All analytic statements hold strictly only in the monopole-free limit; at finite temperature the phases are adiabatically connected crossovers, so no sharp order parameter survives. The extension to quantum $$H = \frac{J}{2}\sum_{\bm{r}} Q_{\bm{r}}^2 + \mu \sum_\ell (S_\ell^z)^2,$$2 spin ice—with transverse fluctuations, quadrupolar degrees of freedom, and quantum loop dynamics—is identified as an open direction rather than treated analytically.

Conclusion

This work provides a complete dual-model description of the $$H = \frac{J}{2}\sum_{\bm{r}} Q_{\bm{r}}^2 + \mu \sum_\ell (S_\ell^z)^2,$$3 pyrochlore spin ice: a 3D $$H = \frac{J}{2}\sum_{\bm{r}} Q_{\bm{r}}^2 + \mu \sum_\ell (S_\ell^z)^2,$$4 mapping explains the deconfinement transition into the U(1) Coulomb phase, a loop-gas/Ising high-temperature expansion explains the confinement transition into the nematic, and a flux-parity argument classifies the topological sectors. Its central result is that thermal monopoles—acting simultaneously as a symmetry-breaking field and as string endpoints—round both transitions into crossovers at any nonzero temperature, a conclusion verified by simulations showing saturated specific-heat peaks and fractionalized flux. The same mechanism is applied to argue that ice-VII and ice-X are adiabatically continuous, underscoring the generality of the gauge-theoretic framework across frustrated magnets and hydrogen-bonded networks.