Simulation of a classical liquid with repulsive nearest neighbor interactions V=1 on a square lattice using Monte Carlo methods. At low temperatures the liquid freezes into a "checkerboard pattern" (e.g. Wigner crystal), since the lattice is chosen to be half-filled. As the temperature is raised above Tc ~ 2, the crystal melts.

We can also add frozen impurities which interact with the particles on neighboring sites with an interaction Q, which can be chosen to be positive or negative. This can be used to investigate how the frozen impurities affect the freezing/melting.

We can enable particle tracking in order to observe trends in particle travel.

Created by Eric Pelz with the help of Graduate Student Yohanes Pramudya and Professors Vladimir Dobrosavljevic and Stratos Manousakis at the National Magnetic Field Laboratory at the Florida State University.