Monte Carlo simulations in the canonical ensemble
Two types of move
- Move a randomly selected ion by a given distance
- Swap a pair of ions (charge switch)
These are accepted automatically if the energy decreases and with a Boltzmann weighted probability otherwise.