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 a **time dependent impurity** in the system, which contains a charge of Q=Q_{0}*cos(?*t), where t is the number of iterations performed and Q_{0} and ? are tunable constants. Once added, the screening charge (Q_{sc}) and the charge (Q) are plotted on a chart as time develops (Q vs -Q_{sc}).

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.