This scene visualizes a 3D Ising model on a cubic lattice with periodic boundary conditions.

Each cube is a spin with value $s = \pm 1$. The Metropolis algorithm updates spins according to the coupling $K$ and chosen neighbor range.

K controls how strongly neighboring spins prefer to align ($K = J/(k_B T)$). At $K \approx 0.221654$ the system is near its critical point.

But note that our system is quite small so finite-size effects are significant. For systems with finite $L$, the critical point shifts to lower values of $K$ and the phase transition is smoothed out.

Use the controls to change lattice size, visualization options, and simulation speed. XYZ coloring maps lattice coordinates to RGB.