Simulation of a classical liquid with long-range Coulomb interactions V(R)=1/R on a partially-filled square lattice, using Monte Carlo methods. By default, the system shows a lattice at 5% filling.

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 create a random disorder in the system, by adding M frozen impurities. These impurities have randomly distributed charges between -W/2 and W/2, where W is the "Strength Parameter."

We can add an external electric field of strength S, which can be chosen to be positive or negative.

We can enable particle tracking, either for all particles or for individual particles, 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.