A version of configuration bias Monte Carlo suitable for the insertion
and deletion of ions in grand canonical Monte Carlo simulations of io
nic solutions has been developed. The method has been optimized and ap
plied to a system where the straightforward insertion and deletion of
ions is prohibitively impractical. The present approach is likely to b
e useful in other systems as well.