The-Powder-Toy VS lammps

The Power Toy - Physics Sandbox. C++, SDL2

I love these kinds of sandbox games, Orb Farm was super nice to have it running directly in the browser for a little procrastination break. Although my favorite is probably The Powder Toy [0], I've been following this project for years.

[0] https://powdertoy.co.uk/

I play The Powder Toy https://powdertoy.co.uk/

Dwarf Fortress N Cave Story Powder Toy Knytt Stories

yes, if you look in the sidebar on Reddit you can see a link pointing you to powdertoy.co.uk

Powder Toy is not only my favorite but available on Android

Look into verlet lists or cell lists for your task, which are typically used with a discrete element method such as LAMMPS for what you seem to be describing.

Your choices are an n-body simulation (e.g., LAMMPS) with Coulomb interactions or, if your electrons are sufficiently sparse, a particle-in-cell (e.g., Starfish). Your best bets for visualization are going to be matplotlib or something more user-friendly like Visit or Paraview. Without a neutralizing background, however, your electrons are just going to repel each other, hit the walls, and disappear - there's not going to be much interesting to visualize. What are you actually trying to simulate? With more information maybe you could receive some more targeted advice.

Is this just a personal /learning project? I see "chemical simulations" and think quantum-espresso, LAMMPS, NWChem, etc. If this is a personal project, I apologize for my bluntness above.

Actually, the LAMMPS website has a similar example: Just hover over the "LAMMPS" logo in yellow and blue.

Modify the file to add what you need. In your case, modify the post_force functions in the files, since those will allow you to define the integration steps for the Vicsek model and its variations.