| Abstract: | The charge line (CL) method had been used in the past to represent the periodic charges in Monte Carlo simulations of ions in a cylindrical pore. In this method, there exists a possible singularity when the edge of the image line overlaps with an ion in the central cylinder. This singularity is more problematic for molecular dynamics when the force is evaluated. Molecular dynamics simulations with the CL method have not been reported in the literature. By replacing the first section of the image charge line with an image point, we show that the CL method can be improved and be applicable in the molecular dynamics simulation of electrolytes in a cylindrical geometry. The modified method is demonstrated to be effective by simulations of a high packing primitive model electrolyte, representing the state of a molten salt. |
