GCEMC Simulations of “Swiss Cheese” Ion Exchange Membranes in Dilute Solutions of Primitive Model 1:1 Electrolytes with Different Sizes of Ions or 2:1 Electrolytes with Different Bjerrum Parameters

Abstract

Monte Carlo simulations using a Markov process corresponding to a (generalized) Grand Canonical Ensemble have been performed for a number of spherical micropores in equilibrium with dilute external bulk solutions of primitive model electrolytes. Dilute solutions of 1:1 electrolytes with a Bjerrum parameter B = 1.546 with cations three times larger than the anions have been simulated. Also, dilute solutions of 2:1 electrolytes with ions of equal size and reduced Bjerrum parameters B_r = 1.546 and 3 have been simulated. The pores are primitive pores with hard walls and the same dielectric permittivity in the wall and in the pore solution. They range from a pore radius = 5 times the mean ionic diameter to 35 times this diameter, and they carry a fixed charge equal to + 5,0 and ?5 elementary charges. The fixed charge is modelled as smoothly distributed on the pore-wall interface. In addition to the electric potential of the interfacial charge and the electric potential of the spherical double layer, a potential Δ between the pore solution and the bulk solution may be deliberately added. For single pores we may take Δ = 0, but then the pore is generally not electroneutral. In a “Swiss cheese” membrane with a lot of (equally sized) pores, the membrane phase has to approach electroneutrality for growing size of the phase. This is approximated by means of a membrane generated potential Δ in each pore (from the electrostatic interactions with the other pores). The potential A so chosen to obtain electroneutrality is the GCEMC Donnan potential. These non-ideal Donnan potentials are compared to the ideal values (with activity coefficients equal to zero). From the mean occupation numbers of cations and anions in the pores, the average pore values of the mean ionic and the single ionic activity coefficients of the ions are calculated. These are very dependent on pore sizes and on the potential in the pore. The excess energy and the electrostatic Helmholtz free energy of the ions in the pores are also simulated directly. The electrostatic entropy is found as the difference.