Discrete aqueous solvent effects and possible attractive forces

We study discrete solvent effects on the interaction of two parallel charged surfaces in ionic aqueous solution. These effects are taken into account by adding a bilinear non-local term to the free energy of Poisson-Boltzmann theory. We study numerically the density profile of ions between the two plates, and the resulting inter-plate pressure. At large plate separations the two plates are decoupled and the ion distribution can be characterized by an effective Poisson-Boltzmann charge that is smaller than the nominal charge. The pressure is thus reduced relative to Poisson-Boltzmann predictions. At plate separations below ~2 nm the pressure is modified considerably, due to the solvent mediated short-range attraction between ions in the the system. For high surface charges this contribution can overcome the mean-field repulsion giving rise to a net attraction between the plates.