Direct Calculation of the Electric Potential Distributions around Ions from High Precision Canonical Ensemble Monte Carlo Simulations of Some Primitive Model Electrolyte Systems

Abstract

Abstract The conditional electric potential around the ions are found by Poisson integration of the conditional charge densities calculated from the radial distribution functions found by the high precision canonical ensemble Monte Carlo simulations of some Z:Z,2:1 and 1:1:1 primitive model electrolyte systems described in earlier studies. In most cases, the correspondance with the solution given by solving the linear Poisson-Boltzmann equation is very nearly perfect – at least not too far from contact. This is so in spite of the fact, that the linearisation condition of Debye and Huckel is violated near to contact in dilute solutions and in spite of the anticipation that the Debye-Huckel approximation should be of little relevance in more concentrated systems, where the Debye length is of the same order of magnitude as the ionic diameters. The linear Poisson-Boltzmann equation even continues to be of relevance in zones around ions of different sizes, where only part of the ions appearing in the solution ...