You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
<pclass="caption">Drag the temperature slider across 2.269 — the exact critical point — and watch the system tip between ordered ferromagnetic domains and disordered paramagnetic noise. The magnetization readout tracks the transition in real time.</p>
80
80
81
-
<p>Wilhelm Lenz introduced the model in 1920 in a paper titled "Beitrag zum Verständnis der magnetischen Erscheinungen in festen Körpern" — a contribution to understanding magnetic phenomena in solid bodies. His setup was simple almost to the point of absurdity: a grid of sites, each holding a spin that could point in only one of two directions, each interacting only with its four nearest neighbors, each either aligned or anti-aligned, energy lowered by agreement and raised by disagreement. He handed the problem to his doctoral student, Ernst Ising, who solved the one-dimensional version in his 1925 thesis and promptly announced that the model didn't work. A chain of spins with only nearest-neighbor coupling exhibits no phase transition. There is no temperature at which alignment suddenly takes over. Ising concluded, reasonably, that the model was too simple to capture ferromagnetism in real materials.</p>
81
+
<p>Wilhelm Lenz introduced the model in 1920 in a paper titled "Beitrag zum Verständnis der magnetischen Erscheinungen in festen Körpern" — a contribution to understanding magnetic phenomena in solid bodies. The setup was deliberately minimal: a grid of sites, each holding a spin pointing in one of two directions, each interacting only with its four nearest neighbors, energy lowered by agreement and raised by disagreement. He handed the problem to his doctoral student, Ernst Ising, who solved the one-dimensional version in his 1925 thesis and concluded the model couldn't capture ferromagnetism in real materials. A chain of spins with only nearest-neighbor coupling exhibits no phase transition — no temperature at which alignment takes hold and persists.</p>
82
82
83
-
<p>He was right about one dimension and wrong about two. The error was understandable: the jump from a chain to a grid is not trivial, and the mathematics needed to solve the 2D case didn't exist yet. It took nineteen years. In 1944, Lars Onsager published the exact solution to the two-dimensional Ising model on a square lattice, deriving the partition function using transfer matrices and the representation theory of a Lie algebra. The result was exact, not approximate, and it contained a specific number: <code>T<sub>c</sub> = 2 / ln(1 + √2) ≈ 2.269</code> (in units where the coupling constant and Boltzmann constant are both 1). Below that temperature, the system spontaneously magnetizes — spins align into coherent domains that persist indefinitely. Above it, thermal fluctuations overwhelm the tendency toward alignment and the net magnetization vanishes. At exactly 2.269, the correlation length diverges and the system displays scale-invariant fractal structure: patches of aligned spins appear at every size simultaneously, with no characteristic length.</p>
83
+
<p>Ising's conclusion held for one dimension. The jump to two dimensions proved harder: the mathematics needed to solve the 2D case didn't exist yet, and it took nineteen years. In 1944, Lars Onsager published the exact solution to the two-dimensional Ising model on a square lattice, deriving the partition function using transfer matrices and the representation theory of a Lie algebra. The derivation produced a specific number: <code>T<sub>c</sub> = 2 / ln(1 + √2) ≈ 2.269</code> (in units where the coupling constant and Boltzmann constant are both 1). Below that temperature, the system spontaneously magnetizes — spins align into coherent domains that persist indefinitely. Above it, thermal fluctuations overwhelm the tendency toward alignment and the net magnetization vanishes. At exactly 2.269, the correlation length diverges and the system displays scale-invariant fractal structure: patches of aligned spins appear at every size simultaneously, with no characteristic length.</p>
84
84
85
-
<p>What made Onsager's result significant beyond the model itself was that it provided the first exact description of a continuous phase transition in a many-body system. The specific heat has a logarithmic singularity at <em>T<sub>c</sub></em> — not a jump, not a spike, but a divergence of the form <code>C ~ −ln|T − T<sub>c</sub>|</code> that approaches infinity as the critical temperature is approached from either side. This kind of critical behavior — diverging susceptibilities, power-law correlations, scale-free fluctuations — turned out to be a universal feature of second-order phase transitions, not a peculiarity of the Ising model. The universality class that the 2D Ising model defines encompasses real physical systems: thin magnetic films, adsorbed gases on surfaces, certain liquid-vapor transitions.</p>
85
+
<p>What made Onsager's result significant beyond the model itself was that it provided the first exact description of a continuous phase transition in a many-body system. The specific heat has a logarithmic singularity at <em>T<sub>c</sub></em> — a divergence of the form <code>C ~ −ln|T − T<sub>c</sub>|</code> that approaches infinity as the critical temperature is approached from either side. This kind of critical behavior — diverging susceptibilities, power-law correlations, scale-free fluctuations — proved to be a general feature of second-order phase transitions, extending well beyond the Ising model itself. The universality class it defines encompasses real physical systems: thin magnetic films, adsorbed gases on surfaces, certain liquid-vapor transitions.</p>
86
86
87
87
<p>The simulation runs a Metropolis-Hastings sweep, introduced in the 1953 paper by Nicholas Metropolis, Arianna Rosenbluth, Marshall Rosenbluth, Augusta Teller, and Edward Teller. Each sweep visits N² randomly chosen sites; at each site, the energy cost of flipping the spin is computed from its four neighbors, and the flip is accepted with probability determined by the Boltzmann factor:</p>
<p>If flipping lowers the energy (dE ≤ 0), the flip is always accepted. If it raises the energy, it is accepted with probability <code>exp(−ΔE / kT)</code>. At high temperature, <code>beta = 1/kT</code> is small and the exponential stays close to 1, so energy-raising flips are accepted often — the system is disordered. At low temperature, the exponential falls steeply and only energetically favorable alignments survive — the system orders. The Metropolis paper was the first application of what is now called Markov chain Monte Carlo, and the original problem it was designed to solve was, in fact, a system closely related to the Ising model: computing equations of state for a hard-sphere gas using random sampling on the MANIAC computer at Los Alamos.</p>
98
98
99
-
<p>What you see in the artifact at the critical temperature is the signature of a universality class rather than a coincidence. The domains at T = 2.269 do not have a preferred scale — there are large clusters and small clusters and clusters within clusters — because at the critical point the system's correlation length has formally diverged. No finite simulation can capture an infinite correlation length, but the qualitative character is visible: the boundary between order and disorder is not a smooth surface but a fractal one. The 2D Ising model at criticality has been mapped exactly to a conformal field theory, and the boundary curves between domains are described by the Schramm–Loewner evolution SLE<sub>3</sub>, one of the random curve families that Oded Schramm introduced in 1999. A lattice of binary spins, solved first in one dimension in 1925 and exactly in two dimensions in 1944, is still generating new mathematics.</p>
99
+
<p>At the critical temperature, the domains do not have a preferred scale — there are large clusters and small clusters and clusters within clusters — because the system's correlation length has formally diverged. No finite simulation can capture an infinite correlation length, but the qualitative character is visible: the boundary between order and disorder is a fractal one. The 2D Ising model at criticality has been mapped exactly to a conformal field theory, and the boundary curves between domains are described by the Schramm–Loewner evolution SLE<sub>3</sub>, one of the random curve families that Oded Schramm introduced in 1999. A lattice of binary spins, solved first in one dimension in 1925 and exactly in two dimensions in 1944, is still generating new mathematics.</p>
100
100
101
101
<divclass="viz-frame">
102
102
<iframesrc="../gallery/grayscott.html" loading="lazy" title="Gray–Scott reaction-diffusion — another system defined entirely by local interaction rules"></iframe>
103
103
</div>
104
104
<pclass="caption">The Gray–Scott reaction-diffusion system belongs to the same family of models: local rules, global pattern. Feed and kill rates play the role that temperature plays in the Ising model — they are the parameters that determine which phase the system occupies.</p>
105
105
106
-
<p>What connects the Ising model to the other artifacts on this site is the argument implicit in its construction: that global, large-scale behavior can be derived entirely from local rules applied uniformly across a lattice. There is no long-range interaction in the Ising model, no global coordinator, no information about the system's overall state available to any individual spin. The phase transition emerges from the statistics of neighborhood interaction repeated across many sites at many temperatures. This is not a metaphor for social dynamics or biological pattern formation — though both have been modeled with Ising-like frameworks — but a mathematical fact about a precisely defined system. The model earns its generality by being the simplest possible system that can exhibit the behavior, stripped of everything that isn't structurally necessary. That economy is what has made it one of the most studied models in all of theoretical physics.</p>
106
+
<p>What connects the Ising model to the other artifacts on this site is the argument embedded in its construction: that global, large-scale behavior can be derived entirely from local rules applied uniformly across a lattice. There is no long-range interaction, no global coordinator, no information about the system's overall state available to any individual spin. The phase transition emerges from the statistics of neighborhood interaction repeated across many sites and many temperatures. Both social dynamics and biological pattern formation have been modeled with Ising-like frameworks, and in each case the explanatory leverage comes from the same place: a precisely defined system in which nothing is assumed beyond the local exchange. The model's reach follows from its constraint — everything that isn't structurally necessaryhas been stripped away, and what remains is enough to produce a phase transition.</p>
0 commit comments