Simulation of confined primitive electrolytes: Application of a new method of summing the coulomb field

Recently, Lekner has presented a new method to sum the Coulomb forces between charged particles of a central system and its images extended periodically in 2 and 3-dimensions. In this paper we apply the new method in canonical ensemble Monte Carlo (CMC) simulations of the primitive electrolyte confined between two planar surfaces: one is charged and the other is neutral. The anions and cations have identical size with diameter d = 4.25 Angstroms and interact with a hard sphere repulsion and Coulomb interaction. In Lekner's method the long-range Coulomb potential is computed from a series of Bessel functions. We have demonstrated that the series converges after about 10 terms and so is computationally simpler that the Ewald sum method. In our simulation, we obtained the density distributions and mean electrostatic potentials of the confined system for the 1:1 electrolyte having concentrations equal to those of 1M and 2M bulk electrolyte and having different surface charge densities. For large separation of confining walls, the canonical ensemble Monte Carlo results agree with previously reported grand canonical Monte Carlo results.