The grand canonical ensemble Monte Carlo method applied to the electrical double layer

The grand canonical ensemble Monte Carlo method is applied to the electrical double layer, using a discrete surface charge distribution and the restricted primitive model for the electrolyte. A systematic examination of the effects of the dimensions of the basic Monte Carlo cell shows that the usual minimum image technique, for determining the potential energy of the system, must be supplemented by long range corrections. Alternatively, very large dimensions or numbers of ions must be employed. The calculated ion distribution functions are compared with the Gouy–Chapman theory for bulk electrolyte concentrations up to 0.2 C/m2. The Gouy–Chapman theory compares very well with the Monte Carlo results for electrolyte concentrations up to 0.1 mole/dm3, but quantitatvive differences appear at 1 mole/dm3. These differences are more pronounced at 2 mole/dm3 for which the Monte Carlo results indicate an inversion of the positive and negative charge distribution functions.

[1]  I. Snook,et al.  The grand canonical ensemble Monte Carlo method applied to electrolyte solutions , 1980 .

[2]  I. Snook,et al.  Solvation forces in simple dense fluids. I , 1980 .

[3]  J. Valleau,et al.  A Monte Carlo study of an electrical double layer , 1979 .

[4]  D. Henderson,et al.  Application of the hypernetted chain approximation to the electric double layer at a charged planar interface , 1979 .

[5]  D. Henderson,et al.  Some exact results and the application of the mean spherical approximation to charged hard spheres near a charged hard wall , 1978 .

[6]  R. James,et al.  Preparation and characterization of amphoteric polystyrene latices , 1977 .

[7]  Christopher W. Outhwaite,et al.  A treatment of solvent effects in the potential theory of electrolyte solutions , 1976 .

[8]  J. Valleau,et al.  A Monte Carlo method for obtaining the interionic potential of mean force in ionic solution , 1975 .

[9]  I. R. Mcdonald,et al.  Monte Carlo simulation of the average force between two ions in a Stockmayer solvent , 1975 .

[10]  Harold L. Friedman,et al.  Study of a Refined Model for Aqueous 1‐1 Electrolytes , 1971 .

[11]  D. N. Card,et al.  Monte Carlo Study of the Thermodynamics of Electrolyte Solutions , 1970 .

[12]  E. Teller,et al.  Monte Carlo Study of a One‐Component Plasma. I , 1966 .

[13]  F. Stillinger,et al.  Statistical Mechanical Theory of Double-Layer Structure and Properties , 1963 .

[14]  S. Levine,et al.  THEORY OF A MODIFIED POISSON-BOLTZMANN EQUATION. I. THE VOLUME EFFECT OF HYDRATED IONS , 1960 .

[15]  W. G. McMillan,et al.  The Statistical Thermodynamics of Multicomponent Systems , 1945 .

[16]  D. Chapman,et al.  LI. A contribution to the theory of electrocapillarity , 1913 .

[17]  C. W. Outhwaite.  Modified Poisson–Boltzmann equation in electric double layer theory based on the Bogoliubov–Born–Green–Yvon integral equations , 1978 .

[18]  S. Levine,et al.  Comparison of theories of the aqueous electric double layer at a charged plane interface , 1978 .

[19]  J. Rasaiah,et al.  Monte Carlo study of two ions in a Stockmayer solvent , 1977 .

[20]  L. Blum.  Theory of electrified interfaces , 1977 .

[21]  C. Radke,et al.  Thermodynamics of the effects of adsorption on interparticle forces , 1973 .

[22]  D. N. Card,et al.  Calculations on the ``Restricted Primitive Model'' for 1–1 Electrolyte Solutions , 1972 .

[23]  H. Friedman,et al.  Ionic Solution Theory , 1962 .