Simulating spin glasses at scale has always been a pain — Python is too slow, and C++ is too painful. Peapods is my attempt at a middle ground: a Monte Carlo package for Ising spin systems where the hot loops are written in Rust and everything else stays in Python, glued together with PyO3. You just pip install peapods and go. ...
Unlike ferromagnetic spin models, a spin glass has a highly non-trivial phase transition profile. For mean field models (the SK spin glass), the glass experiences a replica-symmetry breaking phenomenon that is only fully understood recently. However, for short-range models, almost nothing is known rigorously, including basic problems such as the nature of the glass transition. A potential route towards settling such problems may rely on a suitable generalization of the random cluster model as a universal statistical tool. ...