pith. sign in

arxiv: 2604.04346 · v1 · submitted 2026-04-06 · ❄️ cond-mat.str-el · cond-mat.stat-mech

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

Sena Watanabe , Yukitoshi Motome , Haruki Watanabe This is my paper

Pith reviewed 2026-05-10 20:11 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.stat-mech
keywords pyrochlore spin icetopological phase transitionsdeconfinement transitionCoulomb liquid3D XY criticality3-state Potts modelthermal monopolesspin parity dichotomy
0
0 comments X

The pith

Exact dualities show integer-spin pyrochlore ice undergoes a continuous 3D XY deconfinement transition while half-integer spins remain in a U(1) Coulomb liquid without transition.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The authors develop a self-contained framework using exact dualities to classify topological phases and critical phenomena in classical pyrochlore magnets with arbitrary spin S under competing exchange and single-ion anisotropies. In the small-w regime favoring low spin amplitudes, integer spins exhibit a continuous 3D XY deconfinement transition while half-integer spins stay in the U(1) Coulomb liquid with no transition. In the large-w regime of maximized spin amplitudes, flux is quantized in multiples of 2S and is compatible with the ice rule only for S up to 3/2, where the system maps to the 3-state Potts model whose cubic invariant drives a first-order transition. For S at least 2, monopoles break the discrete mapping, but a hierarchical string fusion cascade exponentially suppresses discrete perturbations as a dangerously irrelevant operator that protects 3D XY criticality; thermal monopoles then act as a symmetry-breaking field that rounds continuous transitions into crossovers while the S=3/2 first-order transition survives to finite temperature and terminates at a critical endpoint.

Core claim

Using exact dualities, the authors establish a spin-parity dichotomy for the small single-ion anisotropy regime: systems with integer S undergo a continuous 3D XY deconfinement transition, while those with half-integer S remain in the U(1) Coulomb liquid phase with no transition. In the large anisotropy limit, the compatibility of quantized flux with the ice rule holds solely for S ≤ 3/2; for S = 3/2 this equivalence to the 3-state Potts model allows a symmetry-permitted cubic invariant that enforces a first-order transition. For S ≥ 2 the mapping breaks due to monopole contamination, yet an exact partition function decomposition demonstrates that hierarchical string fusion exponentially sup

What carries the argument

Exact dualities that map the spin ice models to defect loop gases with emergent Z_{2S} flux conservation, together with the hierarchical string fusion cascade that renders discrete perturbations dangerously irrelevant at the 3D XY fixed point.

Load-bearing premise

The hierarchical string fusion cascade exponentially suppresses discrete perturbations sufficiently to protect 3D XY criticality, and thermal monopoles act purely as a symmetry-breaking field that severs strings without introducing new relevant operators.

What would settle it

A classical Monte Carlo simulation or experiment that detects a phase transition in the small single-ion anisotropy regime for any half-integer spin S would contradict the claimed absence of transitions.

Figures

Figures reproduced from arXiv: 2604.04346 by Haruki Watanabe, Sena Watanabe, Yukitoshi Motome.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic phase diagram in the ( [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Ice rule compatibility at a single [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Hierarchical fusion cascade for [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Defect string topology on the diamond lattice. Ar [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Monte Carlo results for integer spins at [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Monte Carlo results for half-integer spins. (a) Specific heat per spin [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
read the original 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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 3 minor

Summary. The manuscript presents a theoretical framework for classifying topological phases and phase transitions in classical pyrochlore spin ice models with arbitrary spin S, incorporating both exchange and single-ion anisotropy terms. For small values of the anisotropy parameter w, exact dualities are used to argue that integer spins exhibit a continuous 3D XY deconfinement transition, while half-integer spins remain in a U(1) Coulomb liquid phase without a transition. In the large-w regime, the system is mapped to topological loop gases, with compatibility between ice rules and Z_{2S} flux conservation holding only for S ≤ 3/2; specifically, S=3/2 maps to the 3-state Potts model with a first-order transition. An exact partition function decomposition is employed to show that a hierarchical string fusion cascade exponentially suppresses discrete perturbations, protecting the XY criticality as dangerously irrelevant operators. Thermal monopoles are argued to round continuous transitions into crossovers while allowing the first-order transition to persist up to a critical endpoint. These analytical results are corroborated by classical Monte Carlo simulations for S up to 7/2.

Significance. If the central claims hold, particularly the exact dualities and the exponential suppression mechanism via string fusion, this work would provide a significant advance in understanding spin ice physics by offering a unified picture for arbitrary S, including a novel protection mechanism for continuous transitions. The distinction between integer and half-integer spins, the mapping to Potts model for S=3/2, and the finite-temperature fate of transitions represent falsifiable predictions with potential experimental relevance in pyrochlore materials. Strengths include the use of exact decompositions and dualities rather than phenomenological fitting, and the inclusion of Monte Carlo simulations to support the analytics. This could influence future studies on topological phases in frustrated magnets.

major comments (1)
  1. [partition function decomposition section] The section on the exact decomposition of the partition function (invoked to establish the hierarchical string fusion cascade): the claim that this decomposition demonstrates exponential suppression of discrete perturbations (rendering them dangerously irrelevant at the 3D XY fixed point) is load-bearing for the integer-spin branch of the claimed dichotomy and the protection of continuous transitions. The manuscript must explicitly show that the suppression is exponential rather than power-law, holds beyond limiting cases, and does not introduce additional relevant operators when thermal monopoles are included; otherwise the distinction from the half-integer Coulomb liquid phase is not rigorously established.
minor comments (3)
  1. [Monte Carlo simulations] The Monte Carlo simulations section lacks details on system sizes, boundary conditions, error bars, and fitting procedures used to extract transitions or exponents up to S=7/2; these should be provided to allow assessment of the numerical corroboration.
  2. [Abstract] Abstract and introduction: clarify the definition of the parameter w and the precise meaning of 'macroscopic flux quantized to multiples of 2S' with a brief equation or reference to standard spin-ice flux quantization.
  3. [main text] A summary table comparing behaviors across small-w/large-w regimes and integer/half-integer S would improve readability of the classification results.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback on our manuscript. We address the single major comment below, clarifying the partition function analysis and indicating the revisions made to strengthen the presentation.

read point-by-point responses

  1. Referee: The section on the exact decomposition of the partition function (invoked to establish the hierarchical string fusion cascade): the claim that this decomposition demonstrates exponential suppression of discrete perturbations (rendering them dangerously irrelevant at the 3D XY fixed point) is load-bearing for the integer-spin branch of the claimed dichotomy and the protection of continuous transitions. The manuscript must explicitly show that the suppression is exponential rather than power-law, holds beyond limiting cases, and does not introduce additional relevant operators when thermal monopoles are included; otherwise the distinction from the half-integer Coulomb liquid phase is not rigorously established.

    Authors: We appreciate the referee's focus on this central section. The exact decomposition expresses the partition function as a sum over string configurations whose fusion hierarchy is governed by an energy cost linear in unfused segment length. This produces a multiplicative suppression factor exp(−cL) for each fusion level (with L the characteristic loop size), which is exponential rather than power-law in the thermodynamic limit. The argument follows directly from the recursive structure of the decomposition and applies to the full model without restriction to special limits. Thermal monopoles are incorporated by noting that their fugacity remains exponentially small in the relevant temperature window; they sever strings but do not generate new relevant operators at the XY fixed point, preserving the dangerously irrelevant character of the discrete perturbations while rounding the transition to a crossover. We have revised the section to include explicit recursive equations for the suppression factor and an expanded paragraph on monopole effects. revision: partial

Circularity Check

0 steps flagged

No circularity: derivations rest on exact dualities, mappings, and partition-function decompositions independent of target observables

full rationale

The paper develops its dichotomy, S=3/2 Potts mapping, and protection of 3D XY criticality via claimed exact dualities of the spin-ice Hamiltonian, defect-to-loop-gas mappings, and an explicit partition-function decomposition that demonstrates exponential suppression of discrete perturbations. None of these steps reduce by construction to fitted parameters, self-definitional relations, or load-bearing self-citations; the central claims follow from symmetry constraints and the structure of the ice-rule manifold without re-using the predicted transition temperatures or exponents as inputs. The framework is therefore self-contained against external benchmarks and receives the default non-circularity finding.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The framework rests on standard assumptions of classical spin models and known dualities in spin ice; no new free parameters are fitted to data and no new entities are postulated.

free parameters (1)
  • w
    Ratio of single-ion anisotropy to exchange that defines the small-w and large-w regimes; treated as a control parameter rather than fitted.
axioms (2)
  • domain assumption Exact dualities exist between the spin-ice Hamiltonian and XY or Coulomb-liquid models in the small-w regime.
    Invoked to obtain the integer versus half-integer dichotomy.
  • domain assumption Defect structure maps to topological loop gases whose flux quantization is compatible with the ice rule only for S ≤ 3/2.
    Central to the large-w analysis and the S=3/2 Potts mapping.

pith-pipeline@v0.9.0 · 5641 in / 1800 out tokens · 47018 ms · 2026-05-10T20:11:30.793004+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

69 extracted references · 69 canonical work pages

  1. [1]

    Hermele, M

    M. Hermele, M. P. A. Fisher, and L. Balents, Pyrochlore photons: TheU(1) spin liquid in aS= 1/2 three- dimensional frustrated magnet, Phys. Rev. B69, 064404 (2004)

    work page 2004

  2. [2]

    S. V. Isakov, K. Gregor, R. Moessner, and S. L. Sondhi, Dipolar spin correlations in classical pyrochlore magnets, Phys. Rev. Lett.93, 167204 (2004)

    work page 2004

  3. [3]

    C. L. Henley, Power-law spin correlations in pyrochlore antiferromagnets, Phys. Rev. B71, 014424 (2005)

    work page 2005

  4. [4]

    Coulomb phase

    C. L. Henley, The “Coulomb phase” in frustrated sys- tems, Annu. Rev. Condens. Matter Phys.1, 179 (2010)

    work page 2010

  5. [5]

    Fennell, P

    T. Fennell, P. P. Deen, A. R. Wildes, K. Schmalzl, D. Prabhakaran, A. T. Boothroyd, R. J. Aldus, D. F. Mc- Morrow, and S. T. Bramwell, Magnetic Coulomb phase in the spin ice Ho 2Ti2O7, Science326, 415 (2009)

    work page 2009

  6. [6]

    M. J. P. Gingras and P. A. McClarty, Quantum spin ice: a search for gapless quantum spin liquids in pyrochlore magnets, Rep. Prog. Phys.77, 056501 (2014)

    work page 2014

  7. [7]

    Shannon, K

    N. Shannon, K. Penc, and Y. Motome, Nematic, vector- multipole, and plateau-liquid states in the classicalo(3) pyrochlore antiferromagnet with biquadratic interactions in applied magnetic field, Phys. Rev. B81, 184409 (2010)

    work page 2010

  8. [8]

    Pandey and K

    J. Pandey and K. Damle,S= 1 pyrochlore magnets with competing anisotropies: A tale of two Coulomb phases,Z 2 flux confinement and XY-like transitions, arXiv preprint arXiv:2512.11623 (2025)

    work page arXiv 2025

  9. [9]

    Kundu and K

    S. Kundu and K. Damle, Flux fractionalization transition in anisotropicS= 1 antiferromagnets and dimer-loop models, Phys. Rev. X15, 011018 (2025)

    work page 2025

  10. [12]

    M. E. Brooks-Bartlett, S. T. Banks, L. D. C. Jaubert, A. Harman-Clarke, and P. C. W. Holdsworth, Magnetic- moment fragmentation and monopole crystallization, Phys. Rev. X4, 011007 (2014)

    work page 2014

  11. [13]

    Petit, E

    S. Petit, E. Lhotel, S. Guitteny, O. Florea, J. Robert, P. Bonville, I. Mirebeau, J. Ollivier, H. Mutka, E. Ressouche, C. Decorse, M. Ciomaga Hatnean, and G. Balakrishnan, Observation of magnetic fragmentation in spin ice, Nat. Phys.12, 746 (2016)

    work page 2016

  12. [14]

    J. Rehn, A. Sen, and R. Moessner, FractionalizedZ 2 classical Heisenberg spin liquids, Phys. Rev. Lett.118, 047201 (2017)

    work page 2017

  13. [15]

    Blume, Theory of the first-order magnetic phase change in UO2, Phys

    M. Blume, Theory of the first-order magnetic phase change in UO2, Phys. Rev.141, 517 (1966)

    work page 1966

  14. [16]

    H. W. Capel, On the possibility of first-order phase tran- sitions in Ising systems of triplet ions with zero-field split- ting, Physica32, 966 (1966)

    work page 1966

  15. [17]

    Castelnovo and C

    C. Castelnovo and C. Chamon, Topological order in a three-dimensional toric code at finite temperature, Phys. Rev. B78, 155120 (2008)

    work page 2008

  16. [18]

    L. D. C. Jaubert, M. J. Harris, T. Fennell, R. G. Melko, S. T. Bramwell, and P. C. W. Holdsworth, Topological- sector fluctuations and Curie-law crossover in spin ice, Phys. Rev. X3, 011014 (2013)

    work page 2013

  17. [19]

    D. A. Huse, W. Krauth, R. Moessner, and S. L. Sondhi, Coulomb and liquid dimer models in three dimensions, Phys. Rev. Lett.91, 167004 (2003)

    work page 2003

  18. [20]

    Moessner and S

    R. Moessner and S. L. Sondhi, Three-dimensional resonating-valence-bond liquids and their excitations, Phys. Rev. B68, 184512 (2003)

    work page 2003

  19. [21]

    J. Nasu, T. Kaji, K. Matsuura, M. Udagawa, and Y. Mo- tome, Phase transitions to fractionalized Mott insulators in an anisotropic 3d Kitaev model, Phys. Rev. B89, 115125 (2014)

    work page 2014

  20. [22]

    Kitaev, Anyons in an exactly solved model and be- yond, Ann

    A. Kitaev, Anyons in an exactly solved model and be- yond, Ann. Phys.321, 2 (2006)

    work page 2006

  21. [23]

    Fradkin and S

    E. Fradkin and S. H. Shenker, Phase diagrams of lattice gauge theories with Higgs fields, Phys. Rev. D19, 3682 (1979). 20

    work page 1979

  22. [24]

    Senthil and M

    T. Senthil and M. P. A. Fisher,Z 2 gauge theory of electron fractionalization in strongly correlated systems, Phys. Rev. B62, 7850 (2000)

    work page 2000

  23. [25]

    Savit, Duality in field theory and statistical systems, Rev

    R. Savit, Duality in field theory and statistical systems, Rev. Mod. Phys.52, 453 (1980)

    work page 1980

  24. [26]

    J. B. Kogut, An introduction to lattice gauge theory and spin systems, Rev. Mod. Phys.51, 659 (1979)

    work page 1979

  25. [27]

    J. V. Jos´ e, L. P. Kadanoff, S. Kirkpatrick, and D. R. Nel- son, Renormalization, vortices, and symmetry-breaking perturbations in the two-dimensional planar model, Phys. Rev. B16, 1217 (1977)

    work page 1977

  26. [28]

    Watanabe, Y

    S. Watanabe, Y. Motome, and H. Watanabe, Precision Monte Carlo determination of the critical temperature for the antiferromagnetic XY model on a diamond lattice (2026), in preparation

    work page 2026

  27. [29]

    Nahum, J

    A. Nahum, J. T. Chalker, P. Serna, M. Ortu˜ no, and A. M. Somoza, 3d loop models and theCP n−1 sigma model, Phys. Rev. Lett.107, 110601 (2011)

    work page 2011

  28. [30]

    F. J. Wegner, Duality in generalized Ising models and phase transitions without local order parameters, J. Math. Phys.12, 2259 (1971)

    work page 1971

  29. [31]

    See Supplemental Material for detailed derivations of the large-wexpansions, clock-model high-temperature expansion, exact correspondences, exact decomposition framework, and Bethe–Peierls estimates

    work page

  30. [32]

    H. W. J. Bl¨ ote and R. H. Swendsen, First-order phase transitions and the three-state Potts model, Phys. Rev. Lett.43, 799 (1979)

    work page 1979

  31. [34]

    R. J. Baxter,Exactly Solved Models in Statistical Me- chanics(Academic Press, London, 1982)

    work page 1982

  32. [35]

    A. M. Polyakov, Quark confinement and topology of gauge theories, Nucl. Phys. B120, 429 (1977)

    work page 1977

  33. [36]

    Han and S

    Z. Han and S. A. Kivelson, Emergent gauge fields in band insulators, Proc. Natl. Acad. Sci. U.S.A.122, e2421778122 (2025)

    work page 2025

  34. [38]

    J. Lou, A. W. Sandvik, and L. Balents, Emergence of U(1) symmetry in the 3D XY model withZ q anisotropy, Phys. Rev. Lett.99, 207203 (2007)

    work page 2007

  35. [39]

    H. Shao, W. Guo, and A. W. Sandvik, Monte Carlo renor- malization flows in the space of relevant and irrelevant operators: Application to three-dimensional clock mod- els, Phys. Rev. Lett.124, 080602 (2020)

    work page 2020

  36. [40]

    I. A. Ryzhkin, Magnetic relaxation in rare-earth oxide pyrochlores, J. Exp. Theor. Phys.101, 481 (2005)

    work page 2005

  37. [41]

    Nussinov and G

    Z. Nussinov and G. Ortiz, Autocorrelations and thermal fragility of anyonic loops in topologically quantum or- dered systems, Phys. Rev. B77, 064302 (2008)

    work page 2008

  38. [42]

    S.-T. Zhou, M. Cheng, T. Rakovszky, C. von Keyser- lingk, and T. D. Ellison, Finite-temperature quantum topological order in three dimensions, Phys. Rev. Lett. 135, 040402 (2025)

    work page 2025

  39. [43]

    Watanabe, Y

    S. Watanabe, Y. Motome, and H. Watanabe, Continu- ous crossover between high-pressure ice phases VII and X driven by monopole screening: a model study, arXiv preprint arXiv:2603.19620 (2026)

    work page arXiv 2026

  40. [44]

    Castelnovo, R

    C. Castelnovo, R. Moessner, and S. L. Sondhi, Debye- H¨ uckel theory for spin ice at low temperature, Phys. Rev. B84, 144435 (2011)

    work page 2011

  41. [45]

    D. J. P. Morris, D. A. Tennant, S. A. Grigera, B. Klemke, C. Castelnovo, R. Moessner, C. Czarnik, M. Meissner, K. C. Rule, J.-U. Hoffmann, K. Kiefer, S. Gerischer, D. Slobinsky, and R. S. Perry, Dirac strings and mag- netic monopoles in the spin ice Dy 2Ti2O7, Science326, 411 (2009)

    work page 2009

  42. [46]

    L. D. C. Jaubert and P. C. W. Holdsworth, Signature of magnetic monopole and Dirac string dynamics in spin ice, Nat. Phys.5, 258 (2009)

    work page 2009

  43. [47]

    Chern and N

    G.-W. Chern and N. Nagaosa, Gauge field and the confinement-deconfinement transition in hydrogen- bonded ferroelectrics, Phys. Rev. Lett.112, 247602 (2014)

    work page 2014

  44. [48]

    Karsch and S

    F. Karsch and S. Stickan, The three-dimensional, three- state Potts model in an external field, Phys. Lett. B488, 319 (2000)

    work page 2000

  45. [49]

    T. Wada, M. Kitazawa, and K. Kanaya, Locating critical points using ratios of Lee-Yang zeros, Phys. Rev. Lett. 134, 162302 (2025)

    work page 2025

  46. [50]

    Savary and L

    L. Savary and L. Balents, Coulombic quantum liquids in spin-1/2 pyrochlores, Phys. Rev. Lett.108, 037202 (2012)

    work page 2012

  47. [51]

    Shannon, O

    N. Shannon, O. Sikora, F. Pollmann, K. Penc, and P. Fulde, Quantum ice: A quantum Monte Carlo study, Phys. Rev. Lett.108, 067204 (2012)

    work page 2012

  48. [52]

    Topological Phase Transitions and Their Thermodynamic Fate in Arbitrary-SPyrochlore Spin Ice

    O. Benton, O. Sikora, and N. Shannon, Classical and quantum theories of proton disorder in hexagonal water ice, Phys. Rev. B93, 125143 (2016). Supplemental Material for “Topological Phase Transitions and Their Thermodynamic Fate in Arbitrary-SPyrochlore Spin Ice” Sena Watanabe,1,∗ Yukitoshi Motome,1,† and Haruki Watanabe 2, 3, 4, 1,‡ 1Department of Applie...

    work page 2016

  49. [53]

    We call this thevacuum: every link carries|S z ℓ |=S, and the remaining freedom is the signσ ℓ := sgn(Sz ℓ )∈ {+1,−1}

    Vacuum partition function In the large-wlimit (µ→ −∞), the dominant spin configuration on each link is the fully polarized state Sz ℓ =±S, with Boltzmann weightw S2 per link. We call this thevacuum: every link carries|S z ℓ |=S, and the remaining freedom is the signσ ℓ := sgn(Sz ℓ )∈ {+1,−1}. The ice-rule constraint∇ ·S r = 0 [Eq. (3)] then reads SP ℓ∈r e...

    work page

  50. [54]

    We define thedefect num- beron each link as dℓ :=S− |S z ℓ | ∈ {0,1,

    Defect expansion and closed loops Configurations with|S z ℓ |< Son some links represent excitations above the vacuum. We define thedefect num- beron each link as dℓ :=S− |S z ℓ | ∈ {0,1, . . . , S},(8) so thatd ℓ = 0 is the vacuum andd ℓ ≥1 is an excited state. This unsigned defect number is related to the shifted flux variableϕ ℓ :=S z ℓ +S∈ {0,1, . . . ...

    work page

  51. [55]

    ForS≥3/2, every link has a well-defined sign σℓ := sgn(S z ℓ )∈ {+1,−1}(sinceS z ℓ =±(S−1)̸= 0 on defect links)

    Loop-gas representation Consider a configuration in which a subsetG⊆Eof links is excited tod ℓ = 1 while all remaining links carry dℓ = 0. ForS≥3/2, every link has a well-defined sign σℓ := sgn(S z ℓ )∈ {+1,−1}(sinceS z ℓ =±(S−1)̸= 0 on defect links). Letd 0(r) denote the number of excited links incident to vertexr. At each vertex, the ice rule ∇·S r = 0 ...

    work page

  52. [56]

    Higher-defect fugacity and geometric obstruction The Boltzmann weight of a defect-dbond relative to the fundamental bond is xd x1 =w −(2S−d)d+(2S−1) =w −(d−1)(2S−d−1).(18) Ford= 2, this givesx 2/x1 =w −(2S−3), which isexpo- nentiallysuppressed forS≥2. Can all 2Sfundamental strings annihilate at a sin- gle vertex? Each diamond-lattice vertex has exactly tw...

    work page

  53. [57]

    Since at most two fundamental strings can be fused per vertex (using the two outward background links), the process requires a chain of 2S−2 diamond- lattice verticesr 1,r 2,

    Sequential fusion cascade Consider 2Sfundamental loops (ϕ= 1) that are to annihilate. Since at most two fundamental strings can be fused per vertex (using the two outward background links), the process requires a chain of 2S−2 diamond- lattice verticesr 1,r 2, . . . ,r2S−2, connected by 2S−3 in- termediate bonds (see Fig. 2 of the main text). The hierarch...

    work page

  54. [58]

    They fuse into a single outgoing intermediate string with effective chargeϕ= 2 on one inward back- ground link

    arrive on the two outward background links. They fuse into a single outgoing intermediate string with effective chargeϕ= 2 on one inward back- ground link. The divergence-free condition is satis- fied: 2×1−1×2 = 0. 2.Vertexr j (2≤j≤2S−3): One incoming inter- mediate string withϕ=jarrives on an outward background link, together with one additional in- comi...

    work page

  55. [59]

    Energetic cost of the bridge The bridge introduces 2S−3 intermediate bonds car- rying effective chargesϕ= 2,3, . . . ,2S−2. Each in- termediate bond with chargeϕcarries the fugacityx ϕ = w−ϕ(2S−ϕ) [Eq. (9)]. The total bridge penalty—the prod- uct of all intermediate-bond fugacities—is Pbridge = 2S−2Y ϕ=2 xϕ =w − P2S−2 ϕ=2 ϕ(2S−ϕ).(19) Evaluating the sum i...

    work page

  56. [60]

    In the direction-fixed Pauling approximation (Sec

    Endpoint local state sum Consider a vertexrthat serves as the endpoint of a sin- gle fundamental defect string, so thatd 0(r) = 1. In the direction-fixed Pauling approximation (Sec. I B), the sign of the defect link is determined by the string direction, and the three remaining vacuum links carry independent signss ℓ =±1, each corresponding toS z ℓ =±S. W...

    work page

  57. [61]

    q−1X mr=0 ˆBh mr ωmr σr # × Y ℓ

    Generalized loop-gas formula We extend the Pauling analysis of Sec. I B to gen- eral defect graphsGthat may contain both closed loops and open strings. As before, for each connected compo- nent we fix the direction, which determines all defect-link signs; the 2N− |G|vacuum-link signs remain as indepen- dent random variables. The local state sums at degree...

    work page

  58. [62]

    high-temperature expan- sion,

    Explicit Fourier coefficients forq= 3 For the case most relevant to the main text, we spe- cialize toq= 3 andκ=K. There is a single independent fugacityt 1 =t 2. The bond weights are BK(0) =e K, B K(1) =B K(2) =e −K/2,(61) which gives the Fourier coefficients ˆBK 0 = 1 3(eK + 2e−K/2),(62) ˆBK 1 = ˆBK 2 = 1 3(eK −e −K/2),(63) so t1 = eK −e −K/2 eK + 2e−K/2...

    work page

  59. [63]

    vertices at which four bonds carry nonzero current

    Definition For a parameterβ≥1, define thedecoratedXYmodel by ˜Z(β) XY := U(1)X {nℓ} hY ℓ ˜x|nℓ| i βV4({nℓ}) ,(84) whereV 4({nℓ}) denotes the number of degree-4 (crossing) vertices in the current configuration{n ℓ}, i.e. vertices at which four bonds carry nonzero current. This model sums over the sameU(1)-conserving graph ensemble as the standardXYmodel, b...

    work page

  60. [64]

    Exact decomposition of the ice model We now show that the ice-model partition function de- composes exactly into the decorated model (which cap- tures all non-cascade configurations) plus exponentially small cascade corrections. Since ˜Z(β) XY now employs the same ice-model bond fu- gacities ˜x|n|, every non-cascade configuration in the ice model is repre...

    work page

  61. [65]

    parent” and the remainingb= 3 are “children

    Irrelevance of the crossing decoration The decoratedXYmodel (84) differs from the stan- dardXYmodel by the local factorβ V4 at crossing ver- tices. We now argue that this decoration isirrelevantat the 3DXYfixed point in the RG sense. The crossing factorβ V4 is a local interaction that acts at single vertices. In the loop-gas language, it modifies the weig...

    work page

  62. [66]

    Watanabe, Y

    S. Watanabe, Y. Motome, and H. Watanabe, Dualities and topological classification of theS= 1 pyrochlore spin ice, arXiv preprint arXiv:2603.03852 (2026)

    work page arXiv 2026

  63. [67]

    F. Y. Wu, The Potts model, Rev. Mod. Phys.54, 235 (1982)

    work page 1982

  64. [68]

    S. M. Chester, W. Landry, J. Liu, D. Poland, D. Simmons- Duffin, N. Su, and A. Vichi, Carving out OPE space and preciseO(2) model critical exponents, J. High Energy Phys.2020, 142

    work page 2020

  65. [69]

    Hasenbusch, Monte Carlo study of an improved clock model in three dimensions, Phys

    M. Hasenbusch, Monte Carlo study of an improved clock model in three dimensions, Phys. Rev. B100, 224517 (2019)

    work page 2019

  66. [70]

    R. J. Baxter,Exactly Solved Models in Statistical Mechan- ics(Academic Press, London, 1982)

    work page 1982

  67. [71]

    Kikuchi, A theory of cooperative phenomena, Phys

    R. Kikuchi, A theory of cooperative phenomena, Phys. Rev.81, 988 (1951)

    work page 1951

  68. [72]

    J. S. Yedidia, W. T. Freeman, and Y. Weiss, Generalized belief propagation, Advances in Neural Information Pro- cessing Systems13, 689 (2001)

    work page 2001

  69. [73]

    Pandey, S

    J. Pandey, S. Kundu, and K. Damle,Z 3 confined and de- confined Coulomb liquids inS eff = 3/2 pyrochlore mag- nets, arXiv preprint arXiv:2602.23041 (2026)

    work page arXiv 2026