Simulation of a classical liquid with long-range Coulomb interactions V(R)=1/R on a half-filled square lattice, using Monte Carlo methods. At low temperatures the liquid freezes into a "checkerboard pattern" (e.g. Wigner crystal). As the temperature is raised above Tc ~ 0.1, the crystal melts. The melting temperature is dramatically reduced by fluctuation effects induced by strong frustration characterizing the long-range Coulomb interaction.

We can add frozen impurities with tunable charge 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 add an external electric field of strength S, which can be chosen to be positive or negative.

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.