The notebook simulates a 2D Ising Model with nearest neighbor interactions. A system of L x L square lattice with periodic boundary conditions is quenched using the Metrolpolis Algorithm. The following are snapshots of the lattice configurations taken at different temperatures, of a 40 x 40 lattice with nearest neighbor ferromagnetic interactions.
Observables such as the magnetization, magnetic susceptibility, energy and specific heat are also calculated at each temperature step. The Magnetization is zero at high temperatures (disordered phase), wheras it gradually increases (or decreases) to +1 (-1) below the critical temperature Tc ≈ 2.4. The specific heat and heat capacity exhibits a peak at the critical temperature.
